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ABSTRACT. The present position of our knowledge of antiferromagnetism, including 
ferrimagnetism, is reviewed, and some very interesting phenomena concerning the magnetic 
behaviour of certain ferrites and of pyrrhotite are described and explained. 


I—ANTIFERROMAGNETISM 


$1. THE NEGATIVE MOLECULAR FIELD AND 
ANTIPARALLEL SYSTEMS 


E know that in his celebrated theory of ferromagnetism (Weiss 1907) 
represented the interactions between the magnetic moments of neigh- 
bouring atoms by means of a molecular field 


HH, =nJ Ayre a) 


proportional to the magnetization J, m being aconstant. ‘The magnetic suscepti- 
bility s of the substance is therefore given by the Curie—Weiss equation 


PO ee (2) 


The reciprocal of the susceptibility 1/s decreases linearly with absolute temperature 
and, when x is positive, becomes zero at a temperature 0=nC, called the Curie 
temperature, or better the Weiss temperature, below which the substance is 
spontaneously magnetized under the action of its own molecular field: in this 
way one can account for the principal properties of ferromagnetic substances like 
iron, nickel, etc. When 7 is negative s takes a finite value at absolute zero equal 
to —I/n. 

About 1932 I suggested (Néel 1932 a) that eqn. (2), which is based on the 
hypothesis of a random orientation of the atomic magnetic moments in the 
absence of an applied field, ought to break down at low temperatures because 
antiparallel arrangements of atomic moments should be formed under these 
conditions. In fact, the original crystalline lattice ought to break down into two 
sub-lattices A and B, with the moments of the atoms in A directed in one sense and 
those in B directed in the opposite sense, as shown in fig. 1. Energy considerations 
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show that these two sub-lattices should be such that an atom in the sub-lattice A 
should have the maximum possible number of neighbours in sub-lattice B and 
vice versa: in fig. 1 all the neighbours of an atom A are of type B. Under the 
action of a field / this antiparallel arrangement is deformed (fig. 2). 

An elementary calculation on classical lines, applied to elementary moments 
which are capable of being freely rotated, shows that when the temperature 1s 
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Fig. 1. Resolution of a plane lattice into two Fig. 2. Action of a field h on an antiparallel 
sub-lattices. arrangement; the direction of antiferro- 
magnetism is given by D. 


Fig. 3. Variation with temperature of the susceptibility of an antiparallel arrangement of moments; 
the asymptote corresponds to the Curie-Weiss law. * 


markedly below —nC the magnetic susceptibility of the system is constant, 
independent of the temperature, and given by 


Sse Rit a eae . (a9 
where k is a numerical coefficient equal, for example, to } for a body-centred-cubic 
lattice or for the two-dimensional case represented in fig. 1. At high temperatures 
with respectto — nC the antiparallel arrangement is destroyed and eqn. (2) is again 
obeyed. In conjunction with the properties of a diatomic arrangement one 
deduces that 1/s varies with T in the manner indicated by the curve of fig. 3: 
the asymptote of this curve corresponds to the Curie-Weiss law given by eqn. (2). 


* In diagrams for which blocks were already available y and S are sometimes used to indicate 
susceptibility. 
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§2. THE INTERPRETATION OF THE MAGNETIC PROPERTIES 
OF MANGANESE AND CHROMIUM 


The mechanism which has just been described gives at low temperatures 


- something which may be conveniently termed constant paramagnetism. I proposed 


at the time (Néel 1932b) to explain on these lines the magnetic properties of 
manganese and chromium which, in effect, exhibit a constant paramagnetism too. 
large to be explained by the Pauli theory of paramagnetism arising from free or 
partially free electrons. According to the point of view proposed here, one 
therefore attributes to the atoms of manganese or of chromium a definite magnetic 
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Fig. 4. Variation with temperature of the reciprocal of the susceptibility of solid 
solutions of manganese in silver. 


moment. I estimated the value of nC to be 1720°k for manganese and 4150° k 
for chromium, so that the region of temperature at which measurements can be 
made is limited to that where the susceptibility is practically independent of 
temperature. Yet, on decreasing the interactions between the magnetic atoms, 
it should be possible to reach the region of the curve in fig. 3 which is close to the 
asymptote and to obtain direct evidence of the magnetic moments of the substances. 
under consideration. By dissolving manganese or chromium in a non-magnetic 
metal like gold, silver or copper, I showed that one obtains solid solutions obeying 
the Curie-Weiss law: fig. 4 shows, for example, the experimental results for 
manganese dissolved in silver. With the aid of the Curie constant C thus obtained 
one finds for manganese or chromium an atomic moment a little less than 4 Bohr 
magnetons. 
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§3. THE EFFECT OF MAGNETO-CRYSTALLINE COUPLING 


A little later, in 1936, I examined (Néel 1936 a) the effect of magneto-crystalline 
coupling between the lattice and the atomic moments. Let us call D the common 
direction along which the atomic moments of the two sub-lattices are aligned at 
absolute zero, in one sense or the other; this may also be called the direction of 
antiferromagnetism. When the moments are not coupled with the lattice, that is 
to say where there exist only exchange forces of an isotropic kind, the direction D 
is orientated perpendicular to the applied field. When coupling exists D occupies 
a privileged direction and is fixed with respect to the lattice in the absence of an 
applied field, and remains orientated in this direction when the applied field is 
weak. If the applied field is already normal to the privileged direction the results 
are the same as those obtained when coupling is absent; but if the field is applied 
parallel to the privileged direction and if the temperature is close to absolute zero 
the susceptibility is at first zero in weak fields, then for a certain critical field Hy 
decoupling takes place and the direction D turns suddenly to set perpendicular to 
the field: at the same time the susceptibility suffers a discontinuous change and 
takes a finite value, identical with that of the preceding case. Figure 5 gives the 


Fig. 5. Susceptibility as a function of field for different orientations of the field with respect 
to the direction of antiferromagnetism. 


values of the susceptibilities for different directions of the field H with respect to 
the orientation of the privileged direction as a function of the reduced field H/H,. 
This phenomenon and its discontinuity have recently been observed by 
C. J. Gorter and his collaborators (Poulis et al. 1951, Poulis 1952) in the course of 
a study of the magnetic properties of a crystal of copper chloride (CuCl, .2H,O) 
at liquid helium temperatures. : 
When the directions D are randomly distributed, as is probably the case in 
a polycrystalline substance, one recognizes that the susceptibility must vary 
with the field: curve A of fig. 6 gives an example of the thermal variation of the 
susceptibility in weak fields and curve B that in strong fields. The theory of the 
phenomena of D coupling has recently been taken up in Japan by Nagamiya (1951) 
and his collaborators (Yosida 1952). - It is possible that the hitherto unexplained 
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variations of the magnetic susceptibility of platinum as a function of field and 
temperature, studied in detail by Collet and Foéx (1931) and Foéx (1939), have 
their origin in a similar mechanism. 


§4. THE MOLECULAR FIELD APPLIED TO ANTIFERROMAGNETISM 


The method proposed above only gives results applicable to the two extreme 
cases of very low or very high temperatures. In order to obtain a more precise 
idea of the phenomena in the intermediate region I proposed in 1936 (Néel 1936 b) 
to treat negative interactions by the method of the molecular field, considering the 
two sub-lattices A and B to be spontaneously magnetized, each in an inverse sense, 
under the action of two molecular fields 7,,J, OF %pJp, proportional to the magneti- 
zation J, or J,,, of its own sub-lattice and the other 7,,,J, or %p,J proportional to the 
magnetization of the other sub-lattice. The molecular field acting on sub-lattice 
A will be, for example, 

Hi=nt rt dy) | ee en (4) 


Here, on account of symmetry, ”,,=Mpp and m,)=Mpa, but the two coefficients 
n,, and n,,, are not equal because they represent short-distance interactions which 
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Fig. 6. Action of the field on the susceptibility Fig. 7. Variation of the reciprocal of the 
of an antiferromagnetic arrangement : susceptibility of an antiferromagnetic 

A, weak fields; B, strong fields. with temperature. A, the direction of 


antiferromagnetism is perpendicular to 
the applied field; B, parallel to the 
applied field; C, in any direction. 


are different. In addition, the coefficient n,, is essentially negative, otherwise one 
would have ordinary ferromagnetism. 

The analysis of the problem thus shows that the two sub-lattices acquire 
spontaneous magnetizations which are equal and opposite in sign and decrease 
with temperature to become zero at a definite temperature of transition T,. In 
fact, at temperatures below 7, the magnetic moments are ranged in a regular 
manner, while above this temperature complete statistical disorder obtains in the 
substance. From an energy point of view this ordered phase behaves like an 
ordinary ferromagnetic, and is characterized bya specific heat anomaly of which the 
order of magnitude and type of thermal variation are the same as for a ferromagnetic 


substance. 
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As far as the magnetic properties are concerned, since the resultant spontaneous 
magnetization of the ordered phase is zero because the spontaneous magnetizations 
of the sub-lattices exactly compensate one another, the substance always remains 
paramagnetic whatever the temperature. The application of an external field 
perturbs the system and gives rise to a magnetic susceptibility which, in the absence 
of magnetocrystalline coupling with the lattice, remains constant up to the transition 
temperature 7,. Above this temperature the susceptibility obeys the Curie-Weiss 
law represented by eqn. (2) with a negative Curie point 9=nC. On the whole 
the thermal variation of the susceptibility is represented by Curve A of fig. 7. 
One notes that the transition point (6, in the figure) is marked by a sharp change, 

- that is to say by a discontinuity of the rate of change of susceptibility with respect 
to temperature. Such are the essential features of what it is now convenient to 
call antiferromagnetism. 

When coupling exists between the atomic moments and the crystalline lattice 
the susceptibility always follows curve A when the applied field H is perpendicular 
to the direction of antiferromagnetism D, but, when the field H is parallel to this 
direction, Bitter has shown (1938) that one obtains a curve like B (fig. 7) with a 
susceptibility which disappears at absolute zero. Finally, with the directions D 
randomly distributed, one ought to observe a mean curve C, as Van Vleck (1941) 
has pointed out: the susceptibility at absolute zero is thus equal to two-thirds of 
its value at the transition point. 

At the outset even manganese hardly seemed to satisfy the predictions of the 
theory: its specific heat anomaly was rather controversial, but Bizette, Squire and 
Tsai were soon able to prove (1938) that the oxide MnO possessed the properties 
characteristic of an antiferromagnetic body, with a sharp break in the susceptibility 
curve at 122°K accompanied by a specific heat anomaly and an anomaly in its 
thermal expansion. A little later Foéx and Graaf (1939) discovered that Cr,O, 
showed analogous properties. The list of antiferromagnetic substances then grew 
rapidly, mainly due to Bizette and Tsai (1951), with the addition of MnS, MnTe, 
FeO, FeS, CoO, NiO, etc. The properties of these substances agree quite well 
with the curve C given by Van Vleck (fig. 7). In addition, one observes an 
increase of susceptibility with field, below the transition temperature, which 
represents a progressive decoupling for the direction D, predicted by the original 
theory. 

Yet the most satisfactory confirmation of the fundamental ideas was really given 
by Shull. The theory of antiferromagnetism is based in effect on the breakdown 
of the crystalline lattice into two sub-lattices magnetized in opposite senses. 
If atoms possessing magnetic moments directed in opposite senses may be 
considered different from one another, then we have to deal with a kind of super- 
lattice ; such a superlattice cannot be revealed by means of x-rays whose diffraction 
is not affected by the orientation of the moments, but one can show it to exist by 
using a beam of neutrons, when extra diffraction lines are obtained which are 
characteristic of the antiferromagnetic state and which disappear when the 
temperature of the substance is raised above the transition temperature. Shull 
et al. (1951 b) should be given the credit for the experiments which permit us to 
fill in details of the antiferromagnetic structure. Figvre 8, due to these authors, 
allows us to compare the diffraction diagrams of MnO above and below the 
transition point. One observes from the (111) plane the extra diffraction lines due 
to the magnetic superlattice. These experiments have shown that the atoms of 


. 
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_ Manganese are arranged in successive planes perpendicular to a ternary axis and 
that these planes are magnetized alternately in opposite senses: now this is 
exactly the structure which I was led to propose some time ago (Néel 1948) from 
purely magnetic considerations. Neutron diffraction experiments also permit 
us to fix D the direction of antiferromagnetism. This is parallel to a quaternary 
axis for the cases of MnO, MnS, MnTe, CoO and NiO, while it is parallel to a 
ternary axis in FeO. 


$5. PROPERTIES OF SOME ANTIFERROMAGNETIC CRYSTALS 


In the case of a single crystal the relations between the axis of symmetry of the 
crystalline lattice and the direction D can be studied by purely magnetic methods. 
Let us take for example manganese difluoride MnF,, which is a tetragonal uniaxial 
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Fig. 8. Diffraction of neutrons by MnO (C. G. Shull). One notes the extra peak due to the 
antiferromagnetic resolution into two sub-lattices. 


crystal of the rutile type. From the measurements of Bizette and Tsai (1939) on 
powdered specimens, and of Stout and Griffel (1949) on the single crystal, it is 
seen that in the antiferromagnetic region the susceptibility remains practically 
independent of temperature when it is measured in a direction perpendicular to 
the quaternary axis, which corresponds to curve A of fig. 7, whilst it tends to become 
zero at 0°K after the manner of curve C along the direction of the axis: these are 
the properties which one would expect when the antiferromagnetic direction 
coincides with the quaternary axis. The rhombic sesquioxide of iron «Fe,Os is 
equally interesting: it exhibits a fundamental antiferromagnetism superimposed 
upon a very weak ferromagnetism whose origin is not very clear (Néel 1949) but 
which we may neglect here. The study of its magnetic properties suggested to me 
(Néel 1948, 1949) that at ordinary temperatures the successive layers of iron atoms, 
perpendicular to the ternary axis, were magnetized consecutively in opposite 
directions, with a direction of antiferromagnetism perpendicular to the ternary 
axis. ‘These two predictions were confirmed by Shull et al. (1951 b) by means of 
neutron diffraction. 

As has been shown by Morin (1950) and by Guillaud (1951) by magnetic 
studies, the sesquioxide of iron exhibits an anomaly at a temperature of —78°c, 
M. Pauthenet has undertaken the study of this question in my laboratory using a 
single crystal of haematite. The results (Néel and Pauthenet 1952) were roughly 
as follows: in a direction perpendicular to the ternary axis the susceptibility 
remains practically constant from liquid air temperatures up to the transition 
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temperature 675°c.. On the other hand, in the direction of the axis the suscepti- 
bility is at first very weak at low temperatures, taking a value in the region of 
— 78° c almost equal to the value along the perpendicular direction and maintaining 
that value up to the transition point. These facts may be explained very easily if 
one supposes that the direction of antiferromagnetism coincides with the ternary 
axis below —78°c and moves into the perpendicular plane, with freedom of orienta- 
tion, at higher temperatures. This change in orientation may be likened to those 
undergone by the directions of easy magnetization in uniaxial ferromagnetic 
bodies like cobalt or MnBi. It is to be noted that in the case of the latter substance 
Guillaud (1943) has shown that there exists a plane of easy magnetization perpen- 
dicular to the hexagonal axis below 84°, whilst the hexagonal axis becomes an 
axis of easy magnetization above that temperature. 


§6. THE APPROXIMATE NATURE OF THE THEORY 

As we have just seen, the theory of antiferromagnetism, based on the reduction 
to two sub-lattices and extension thereto of the hypothesis of the molecular field, 
has permitted us to predict and to coordinate a large number of experimental 
facts: one could hardly expect more. On one hand, its classical nature does not 
allow it to be rigorously applied except to atoms possessing a relatively large 
magnetic moment; on the other hand, the representation of interactions by a 
molecular field is the better the larger the number of neighbours which act magneti- 
cally upon a given atom. It is, therefore, a theory of limited application veg 
badly adapted, for example, to the case of a linear chain of atoms, of spin s=4, 
as Kramers (1952) has noted. One might also raise the same objections to 
the molecular field theory of ferromagnetism, whose utility one nevertheless 
appreciates: it is only necessary to bear in mind the limited and approximate 
nature of the theory which one uses. 


II_FERRIMAGNETISM 


§1. IMPERFECT ANTIFERROMAGNETISM 

In an antiferromagnetic substance the two sub-lattices are physically and 
essentially identical, so that their spontaneous magnetizations compensate one 
another exactly, and the resultant spontaneous magnetization is consequently 
exactly zero. But one could very well imagine, for some reason or other, that 
the two sub-lattices are not exactly equivalent; the spontaneous magnetization 
of one of the sub-lattices will then be greater than that of the other, and the resultant 
spontaneous magnetization will be different from zero, as in a ferromagnetic 
substance. ‘This is the reason why the iron sulphide FeS is a normal antiferro- 
magnetic, where the number of atoms of sulphur is equal to the number of atoms 
of iron, for the two sub-lattices are then exactly the same. But when there is a 
net excess of sulphur, as in the case of pyrrhotite FeS,.,,, the two sub-lattices are 
different and ferromagnetism is exhibited: we will return to this example later. 
It is probable that the small superimposed ferromagnetism of «Fe,Oy, arises in an 
analogous manner. 

However, the outstanding example of these phenomena is that of the ferrites 
Fe,0,;MO, where M is a bivalent metal, such as copper, nickel, cobalt, iron, 
manganese, etc. The crystalline lattice is then of the spinel type Al,O,MgO. 
The metallic ions may occupy two different categories of sites from the crystallo- 
graphic point of view: A sites surrounded by four atoms of oxygen and B sites 
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surrounded by six atoms of oxygen. For each molecule there are one A site and 
two B sites. Verwey (Verwey and Heilmann 1947, Verwey and de Boer 1936) 
showed by x-ray studies that the ferromagnetic ferrites which are of interest to us 
possess a structure which is called inverse, in which the A site is occupied by a 
Fe’* ion and the two B sites by the M?* ion and the other Fe*+ ion. Experiment 
shows that magnetically these substances are characterized by a feeble saturation 
magnetization with respect to the total magnetic moment of the ions which are 
present, and above the Curie point by a thermal variation of the susceptibility 
which is distinctly different from that of ordinary ferromagnetics; by representing 
the reciprocal of the susceptibility 1/y as a function of the temperature 7’ one 
obtains a hyperbola instead of astraight line, as is evident in the early measurements 
of Mlle Serres (1932) and the more recent results of Fallot and Maroni (1951) 
(cf. fig. 12). ‘This was an irritating puzzle for the theoreticians. 

In order to interpret these properties I supposed (Néel 1948) that in these 
substances the principal magnetic interactions were those of the ions situated in 
A sites with the ions situated in B sites, and that the interactions AB were negative. 
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Fig. 9. Comparison of molecular moments predicted by theory and experimental 
values for certain ferrites. 


One may easily show that if the interactions AB are sufficiently large with 
respect to the interactions AA between atoms situated in A sites and the inter- 
actions BB between atoms situated in B sites, two magnetic sub-lattices are formed, 
composed of the ions situated in sites A and of the ions situated in sites B 
respectively: these two sub-lattices are both magnetized to saturation at absolute 
zero, but in opposite senses. As these two sub-lattices are essentially different, 
there exists a resultant spontaneous magnetization which is not zero and a ferro- 
magnetism which arises, paradoxically enough, from negative interactions. In 
view of this fact and of the rather special properties which follow from it, I thought 
it appropriate to give the name ferrimagnetism to this mechanism. 


§2. PROOFS OF FERRIMAGNETISM 
It follows from the above hypothesis that at low temperatures the magnetic 
moments of the two ferric ions in the molecule compensate one another because 
they belong to different sub-lattices: the molecular moment is therefore merely 
equal to that of the ion M?+. Figure 9 gives in Bohr magnetons the saturation 
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molecular moments of anumber of ferrites (black circles) from determinations made 
in different laboratories. ‘The hatched region corresponds to the predicted 
theoretical moments, taking account of probable orbital moments ; the continuous 
line refers to moments due to spin only. ‘The agreement is very satisfactory 
except in the case of copper: we shall see later a reason for this apparent 
exception. 

Neutron diffraction has brought an even more direct confirmation of the 
ferrimagnetic structure: Shull et al. (1951) has shown that in the case of 
magnetite the ions situated in A sites possess magnetic moments oppositely 
directed to the resultant total magnetization. 

Another striking confirmation of the theory is given by the variation with 
concentration of the molecular moment of certain mixed ferrites in which a 
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Fig. 10. Comparison of molecular moments predicted by theory and the experimental values 
for the series of mixed ferrites of nickel and zinc: G, Guillaud’s values; P, Gorter’s values. 


portion of the ions M is replaced by non-magnetic Zn?*+ ions. For zinc ferrite is a 
normal ferrite, i.e. a ferrite in which the two Fe®+ ions are in the two B sites, 
whilst the Zn?* ion occupies the A site, for which it has a special affinity. When 
in an inverse ferrite one replaces an M ion in a B site by a Zn ion, that Zn ion 
goes into an A site causing the movement of an Fe*+ ion from a B to an A site. 
This results in the reversal of the magnetic moment of this Fe®*+ ion, which is 
equal to tive Bohr magnetons. The net effect is that the replacement of an 
M ion by a Zn ion produces an increase of moment equal to 10—m Bohr 
magnetons, m denoting the moment of the M ion. ‘These predictions were 
verified in a remarkable manner by the experiments of E. W. Gorter (1950) at 
Eindhoven and by Guillaud (1949) at Bellevue. Figure 10 shows as an example 
the results obtained with mixed ferrites of nickel and zinc: on the NiO side the 
experimental moments fall remarkably well upon the theoretical line AB. 
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The case of copper ferrite, to which we have just referred, is particularly 
interesting: here we have an inverse ferrite in which the affinity of the M ion, 
namely the Cu?* ion, for the B sites is particularly weak. It therefore happens 
that at high temperatures these ions are randomly distributed over the A and 
B sites: onthe average one-third are on the A sites and two-thirds on the B sites. 
When the temperature is lowered they all progressively reoccupy the B sites. 
Quenching permits us to keep in unstable equilibrium the different intermediate 
pel | 
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Fig. 11. Saturation magnetization of copper and magnesium ferrites as a function of quenching 
temperature; theoretical curves shown by full lines. Upper curve Fe,0,;CuO; lower 
curve Fe,O;MgO. 


states and, therefore, makes manifest the variations of molecular moment as a 
function of the quenching temperature. Pauthenet and Bochirol (1951) have 
studied this phenomenon in my laboratory and have shown that the observed 
moments agree well with the moments predicted by the statistical theory. In 
fig. 11 the observed values are compared with the theoretical curves (full lines) 
in the case of copper ferrite and the closely analogous magnesium ferrite. ‘These 
atomic migrations within the crystalline lattice, thus demonstrated by magnetic 
methods, have been confirmed by x-ray analysis by one of my collaborators, 
M. Bertaut (1951). 

In conclusion, we see that it is possible to predict the molecular moment of 
ferrimagnetic bodies with much greater exactness than that of the ordinary 
ferromagnetics, iron, nickel or cobalt. ‘This was a priori very unexpected. 


§3. APPLICATION OF THE MOLECULAR FIELD THEORY 
TO FERRIMAGNETISM 


We have thus been able to extend the theory of ferrimagnetism by generalizing 
the hypothesis of the molecular field and by representing the different interactions 
present, AA, AB and BB, by three molecular field coefficients. We thus show that 
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the susceptibility s above the Curie point is represented by a hyperbolic type of 


expression of the form 
is=TiCe lyse T= 8) 90h COU ae RARE (5) 


where C is the Curie constant corresponding to the given ions in the absence of 
interactions, and x, o and @ are functions of the three molecular field coefficients 
2, 4,N,p and Ny_. Fallot and Maroni have shown that the susceptibility of ferrites 
can be effectively represented by an expression of the form (5). We merely note, 
as is obvious from fig. 12, that in the region of the Curie point there occur some 
rather marked discrepancies between experiment and theory whose origin is 
analogous to that of the discrepancies observed in the case of ordinary ferro- 
magnetics between the paramagnetic Curie point and the ferromagnetic Curie 
point: they are simply due to the effect of the fluctuations of the molecular field 
which I studied previously (Néel 1932 a). In fact, when one represents magnetic 
interactions by a molecular field it is necessary to suppose that this field fluctuates 
about the mean value in space and time with an amplitude which is the greater 


Fig. 12. Experimental values (full lines) of the reciprocal of the susceptibility of the ferrites of 
nickel, cobalt and iron as functions of temperature. The theoretical hyperbolae are shown 
by the broken lines. 


the smaller the number of neighbouring atoms reacting with one another. Con- 
sequently, the disappearance of the ordered state takes place at a temperature a 
little lower (20° to 50°c) than the theoretical temperature corresponding to a 
molecular field which does not fluctuate. 

The study of paramagnetic properties above the Curie point thus provides 
values of the three coefficients of the molecular field which, with the values of the 
moments of the ions, are sufficient to specify completely, to a first approximation, 
a ferromagnetic ferrite; one can thus calculate the curve of the thermal variation 
of the spontaneous magnetization without recourse to any other data. The com- 
parison of these curves with the experimental curves provides a proof of the value 
of the method of the molecular field applied to ferrimagnetism. Figure 13, taken 
from a work by Pauthenet and Bochirol (1951), allows us to compare the experi- 
mental values of the spontaneous magnetization with the theoretical curves for the 
three ferrites of iron, cobalt and nickel. The agreement is very satisfactory. 
We notice that these curves appear very different from the corresponding curves. 
for ferromagnetics, being in general more flattened than these. 
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§4. FEATURES OF FERRIMAGNETISM 


One must emphasize that ferrimagnetism is essentially associated with a 
paramagnetic law of the hyperbolic form (eqn. (5)); the hyperbola may, of course, 
degenerate into its asymptotes. In order to conclude from the examination of the 
experimental data that ferrimagnetism rather than ferromagnetism exists, it is 
therefore necessary that the variation of the susceptibility with temperature above 


12) 0,2 0,4 O6 0,8 11/0 
FeONiO: Courbe th€orique et points P 
<3 ex perimentaux 


O6 
7 


0,2 0,4 O6 oy) 1 ve 
FeQFeO - FeOC 00 : Courbes théoriq ues et P 
oi points expe@rimentoux: 


Top diagram: Fe,O,NiO. 
Bottom diagram: (1) Fe,0;FeO; (2) Fe,.O,;CoO. 


Fig. 13. Spontaneous magnetization of the ferrites of nickel, cobalt and iron as a function of 
the reduced temperature; experimental points and theoretical curves calculated from 
paramagnetic data. ; 


the Curie point should be of the form (5). .But quite recently Zener (1952) 
proposed a breakdown of the iron lattice into two simple cubic sub-lattices, one 
composed of Fe*+ ions of five Bohr magnetons, the other of Fe~ ions of one Bohr 
magneton, the latter being orientated in the opposite sense to the former: the iron 
would then be a ferrimagnetic. Independently of the objections based on energy 
considerations which such a conception would raise, if it were true one should find 
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above the Curie point a hyperbolic variation of the reciprocal of the susceptibility, 
1/s. In practice one actually finds the ordinary Curie-Weiss law. It therefore 
does not seem possible to accept Zener’s point of view. 


§5. OTHER FERRIMAGNETIC SUBSTANCES 


Many other ferrimagnetic substances are recognized today in addition to the 
ferrites, for example, the antimonide of manganese, Mn, Sb, studied at low temper- 
atures by Guillaud (1943). In the crystalline lattice of this substance there are in 
effect two kinds of atoms of manganese possessing different electron structures: 
some are in a 3d’ state with three Bohr magnetons, the others are in a 3d° state with 
five Bohr magnetons pointing in the opposite sense to the former: the theoretical 
molecular moment would thus be equal to two Bohr magnetons, whilst experiment 
gives 1:87. The form of the thermal variation of 1/s confirms this interpretation, 
as shown in fig. 14, which reproduces the experimental points of Mlle Serres 
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Fig. 14. Variation of the reciprocal of the susceptibility of Mn,Sb with temperature (experimental 
points with fitted theoretical curve). 


(1947), while the theoretical hyperbola is given by the full curve. One notices the 
difference between the theoretical Curie point obtained by extrapolation from 
the hyperbola and the true Curie point: again this is a result of the effect 
of fluctuations of the molecular field already mentioned. The curves of the 
thermal variation of the spontaneous magnetization are likewise different from 
those of ordinary ferromagnetics: but the experimental values on the contrary 
approximate very well to the values calculated from the theory of ferrimagnetism if 
we use the values of the coefficients deduced from the study of paramagnetism 
(Néel 1948). 

Another good example, of more interest, is that of the sulphides of iron. The 
sulphide of iron FeS in stoichiometric proportions is an antiferromagnetic sub- 
stance with a transition point at 325°c. It belongs to the Ewald type B8, like 
NiAs. Let us number consecutively the planes of iron atoms perpendicular to 
the ternary axis ; we may then imagine that one of the sub-lattices is formed by the 
planes of even order and that the other is formed by the planes of odd order. 

For the non-stoichiometric sulphides the magnetic properties change gradually 


as the proportion of sulphur increases, and for the compound Fed, osanprox 
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ferromagnetic properties appear. In particular, the saturation magnetization 
seems to pass through a maximum in the region of the composition FeS,.,,, which 
is approximately Fe,S,: one is reminded of pyrrhotite. The saturation magneti- 
zation of pyrrhotite is small, and of the order of 0-36 Bohr magnetons per atom of 
iron: it is reasonable to suppose that we have here a case of ferrimagnetism due to 
lack of equality between the two sub-lattices, which were originally the same in 
FeS. ‘This inequality should probably be associated with the unequal proportion 
of atoms of iron and sulphur in pyrrhotite. Moreover, as sulphur must always 
remain bivalent, the composition FeS,,, must necessarily contain 1 —2x of Fe?+ 
ions and 2x of Fe*+ ions. Also, one may suppose, on one hand with Hirone 
and T’suya (1951) and on the other hand with Yosida (1951), that the ferric ions 


Fig. 15. Lattice structure of pyrrhotite. 


collect on the iron planes of even order perpendicular to the ternary axis, and 
that the planes of odd order contain only ferrous ions. In this way one obtains 
the necessary lack of symmetry and the correct order of magnitude for the 
saturation moment. 

Unfortunately, this explanation is not very satisfactory, for it neglects one 
essential fact: the crystalline lattice of FeS,,., must be an incomplete lattice with 
x unoccupied places which are normally occupied by Fe?*+ionsin FeS. I therefore 
think that these vacancies, or holes, must arrange themselves regularly on iron 
planes of even order and give rise to the required lack of symmetry. If this idea 
is correct, this lack of symmetry must correspond to a kind of superstructure 
discernible by x-rays, because the scattering powers of the holes and of the iron 
ions are different. I therefore asked M. Bertaut, the x-ray specialist in my 
laboratory, to find if pyrrhotite actually had the structure of NiAs and if the 
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diffraction pattern did not contain by accident some supplementary lines. By 
making exposures which were sufficiently prolonged, Bertaut (1952) has actually 
found 240 diffraction spots which may be attributed to this superstructure and 
which indicate an ordered distribution of the holes: consequently it appears that 
this is the first known example of a structure containing holes in ordered arrange- 
ment. Figure 15 shows the structure obtained, corresponding to the formula - 
Fe,S, JT, where T'isahole. ‘This figure shows the arrangement of the iron atoms 
and the holes in the planes of odd order. For the sake of clarity, planes of iron 
of even order, which contain no holes, have net been depicted, neither have the 
intermediary planes of sulphur. One of the sub-lattices is therefore probably 
formed by the planes of iron with holes, the other by planes without holes. One 
would like to study the position of the ferric ions, but this is not possible by x-ray 
methods. In any case, as the two sub-lattices are unsymmetrical, ferromagnetism 
necessarily follows. ‘ 

Inthe paramagnetic region the susceptibility of pyrrhotite shows the character- 
istic hyperbolic behaviour of ferrimagnetism, but as Benoit (1952) showed at 
Grenoble, a new phenomenon appears at 600°c: the susceptibility undergoes a 
sharp discontinuity and becomes much more weak. It possesses, above this 
temperature, a value and a thermal variation which very much resemble those of 
FeS. It therefore appears most reasonable to attribute this discontinuity to the 
disappearance of the superstructure. Above 600°c the holes distribute them- 
selves randomly on all the iron planes, and the two sub-lattices should become 
identical; the substance should then show properties very close to those of FeS, 
as is found by experiment. 

§6. CONCLUSION 

It follows from this account that the idea of ferrimagnetism and its treatment 
by the method of the molecular field are on the whole fairly fruitful in spite of their 
somewhat undue simplicity. Undoubtedly it would be desirable to apply more 
highly developed methods to these problems, but when we consider that there 
do not yet exist moderately rigorous and satisfactory solutions for the case of a 
simple cubic lattice occupied by ions with spin s =4, and with one single exchange 
integral, it seems very unlikely that we shall get an immediate solution of the 
problem of spinels, with two different kinds of atoms with different spins, greater 
than 4, and nine exchange integrals, all different! For the time being it is not 
unwise to use a simple method even if it is a little rough and ready. 

Finally, let us note that the studies of magnetic properties of antiferromagnetic 
and ferrimagnetic substances have enabled us, from the experimental data, to 
prove the necessary existence between magnetic atoms of indirect exchange 
interactions in which an intermediary, non-magnetic atom plays an essential part; 
there exists in this way magnetic coupling by the intermediary action of atoms of 
sulphur or of oxygen of the type Fe~-O—Fe, Mn—O--Mn, Fe-S--Fe, etc. This is 
the superexchange whose possibility Kramers (1934) previously envisaged on 
theoretical grounds. 
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ABSTRACT. ‘The theory of ferromagnetism is often developed on the assumption that 
the exchange energy responsible for the magnetization is adequately represented as a sum 
of exchange integrals between pairs of electrons with mutually parallel spins. ‘This is 
equivalent to using for the total wave function a determinant built up of one-electron 
functions. For the magnetized states of interest in ferromagnetism a single determinant 
wave function must be regarded as inadequate on account of the so-called spin or exchange 
degeneracy, which requires us to use instead appropriate linear combinations of the degen- 
erate determinants. In the present paper a model is presented for which the results of 
neglecting or taking account of this spin degeneracy are identical. ‘The model can be handled 
exactly and leads to equations for the free energy, magnetization, etc., which are a general- 
ization of Stoner’s theory of ferromagnetism to a system with two energy bands; these 
could be the overlapping 3d and 4s bands in the iron transition series. ‘The calculations 
are not necessarily restricted to ferromagnetic metals and would apply equally well to those 
that are paramagnetic. 
* Now at Department of Physics, University of Illinois, Urbama, U.S.A. 
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§1. INTRODUCTION 

T is possible to distinguish two main stages in the calculation of the quantum 
| mechanical properties of the valence electrons in metals. In the first we 

calculate the energy of an electron moving in the ion-core field and in the average 
field of the remaining valence electrons; the energy of all the electrons together i is® 
then taken to be simply the sum of the resulting one-electron energies. This 
stage is represented by the elementary band theory of solids (Seitz 1940, especially 
Chapter 8). Thermal properties can be handled by using the simple Fermi—Dirac 
statistics. Inthe next stage we take the one-electron wave functions, which may be 
those found by some stage-one calculation, build up a product wave function ‘Y 
with the correct symmetry properties and then calculate the energy of the selected 
configuration by using the true N-electron Hamiltonian. This more elaborate 
procedure is less often used in discussing electronic thermal properties, though 
one is forced in principle to use it in the theory of ferromagnetism. Generally 
the simple Fermi—Dirac statistics are not applicable in this approximation 
(Wohlfarth 1950). 

Now inthissecond stage it is commonly assumed that the spin dependent energy 
of a state is asum of exchange integrals for all those pairs of electrons with mutually 
parallelspin. Thisis equivalent to describing the whole system by a wave function 
which is a determinant of one-electron functions, each one-electron function being 
simply the product of an orbital space function and a spin function. The orbitals 
may be assumed localized about particular atoms as in the Heisenberg (1928) 
theory, or non-localized as, for example, in Slater’s calculations (1936). In both 
cases the electronic interaction term 2; —;e?/r,; in the Hamiltonian is treated as a 
perturbation, the eigenfunctions of the unperturbed Hamiltonian each being a 
determinant of one-electron functions. When however a given state (or con- 
figuration) contains only one electron in some orbitals, there will be a number of 
degenerate unperturbed determinant functions corresponding to the different ways 
of assigning spin functions (either « or f) to these singly occupied orbitals. This 
degeneracy is known as spin or exchange degeneracy. It follows that the correct 
wave function for calculating the perturbation energy due to the electron interactions 
is not a single determinant, as is frequently assumed in ferromagnetism, but some 
linear combination of the degenerate determinantal functions. ‘The method of 
calculating the coefficients and the perturbation energy is familiar from the theory 
of atoms and molecules (see for example Slater 1929, or Condon and Shortley 
1935) and leads to a classification of states by the eigenvalues of the square of the 
total spin vector operator and of the z component of the spin vector operator. 
Furthermore, Dirac (1929) showed that this problem was equivalent to a similar 
degeneracy problem, in which the perturbation part of the Hamiltonian was a sum 
of terms, each proportional to the scalar product of the spin vector operators of two 
electrons in singly occupied orbitals. ‘This leads to a very convenient way of 
solving some exchange degeneracy problems, and one which we shall employ in 
the present paper. A very clear account of this method has recently been given by 
Corson (1951). 

It is possible in principle to proceed from here to a calculation of magnetic and 
thermodynamic properties by constructing the partition function or sum-over- 
states, although for the large systems arising in the theory of ferromagnetism it is. 
generally only possible to evaluate both the energy levels and the partition sum 
approximately. Bloch (1932) showed that to neglect spin degeneracy is justifiable 
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at sufficiently high temperatures, where the partition function goes over into the 
function that would be obtained if we used single determinant wave functions 
rather than the correct linear combinations. Many accounts of the theory of 
ferromagnetism start either explicitly or implicitly from this high temperature 
approximation. It does not seem to have been noticed however that there is a 
model, probably nearly realized in practice, which can be handled exactly, and for 
which the results of including or neglecting the spin degeneracy are the same. 
We specify this model by the requirement that the available orbitals divide into 
two groups, A and B, such that the exchange integral between any two orbitals. 
depends only on which group(s) they belong to. This includes the special case 
of a model with only one group or one common exchange integral. On this basis 
we shall derive exact equations for magnetization, free energy, etc., which prove 
to be a generalization of those used by Stoner (1938, 1939) in his theory of collective 
electron ferromagnetism. ‘These equations are derived in the next section, and 
in the final §3 we list several observations on them and on their relation to previous 
work in this field. 


§2. EXPRESSION FOR THE FREE ENERGY 

We suppose that the orbitals ,(r) of the ‘ferromagnetic’ electrons form a set of 
localized orthonormal functions extending throughout the magnetic ‘domain’, 
and that the inner shell electrons may be replaced by a suitable potential field. 
The sum of this field and that due to the nuclei we denote by U(r). For definiteness 
the y,(r) may be supposed to be Bloch tight-binding wave functions (see Mott and 
Jones 1936), though in what follows we shall not need such an assumption explicitly. 
With each orbital ;, may be associated an re e,, thus: 


ae | i 01 Sm V24 Up bi(r) dr, eke (2.1) 


This may be regarded as the sum of the kinetic energy and the energy of interaction 
of the orbital ¥, with the lattice. We also form the Coulomb and exchange integrals, 
denoted by J,; and K,, respectively in Slater’s (1929) notation : 


ies | i '(ry) vi (ry) — = DO aduarsardre i 0 a. (2.2) 


Dee [|e ob; (ra) - BAK, (hei Gri ats.6. | ta eee (2.3) 


Dirac’s result for the exchange problem may now be applied in the following form 
(Corson 1951, § 81) ; insolving the secular problem arising from the spin degeneracy 
the perturbation term in the Hamiltonian may be replaced by the following 
operator : 


v= % Diy —4(1+4s,.s,)K,]+ PACE: ee ad eK.) + Sil BERS 
t< le i 


eaov 


in which doubly-filled orbitals are denoted by a Greek letter and singly-filled orbitals 
by a Roman letter. s, is the spin vector operator for the electron in orbital 2, 
and s,, is the g component of this operator. ‘The Bohr magneton has been denoted 
by 8. The last term in eqn. (2.4), being the operator for the energy in the external 
magnetic field H, is not really part of the Dirac result, which concerns only the 
electronic interactions, but it may conveniently be added into the perturbation 
58-2 
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operator V. ‘The possible energies for the various spin states of a particular 
configuration are the eigenvalues of the matrix of V. 

We now apply this result to a model for which the whole set of available orbitals 
may be divided into two groups, A and B, such that all pairs of A orbitals have a 
common exchange integral K,, all pairs of B orbitals have a common exchange 
integral Ky and all pairs formed from one A and one B orbital have a common 
exchange integral K,,. The spin dependent part of V then becomes 


—22 s,.$,Ky 
i<j 
=—2K, & 84;.84;-2Kg 2% $;-Spj—2K ap % Saz- Sp; 
Ai<Aj Bi<Bj Ai, B 
=—2K, & $y;.8,;—2Kg Sp; -Spj—2KapSa - Sp 
Ai<Aj Bi<Bj 
= —K,(S,?— 2 8,7)—Kp(Sp?— USp*)—2KypSa-Sp,  -+s--- (2.5) 
Ai Bi 
where S,=2s,, -and “Sy=23p,0 9 eee ee (2.6) 
Ai Bi 


are the total spin vector operators for those electrons in A and B orbitals respectively. 
Thus $,?and $,? commute with the right-hand side of eqn. (2.5) and, of course, with 
S,,+Sp,=S, and $?=(S,+S$,)?; therefore, since the unperturbed Hamiltonian 
is spin-independent, S$, and S,?, as well as S? and S,, are constants of the motion 
with eigenvalues S,(S,+1) and Sp(Sg+1) respectively. For a configuration 
with m A orbitals and m B orbitals singly occupied (z and m are assumed even for 
convenience) the numbers S, and Sg may be any of the integers from 0 to 2/2 
and 0 to m/2 respectively. ‘The eigenvalues of S$? are S(S+1) with the allowed 
values of S being the integers from S, + Sp to |S, —S,|, in the usual way for the 
sum of two angular momentum vectors (Corson 1951, §21). Lastly, the eigen- 
values of S, are.S, S—1,... —S+1, —S. Since 2S, .S,=S?—S,?—S,?, we 
see that the secular matrix of V has eigenvalues 


2 2 
=O RY Pap WH i ee NE Yip a0 515 joc a ok —— w =e 
i<j di, << ue 
+5 (Kap— Kp) ta Kaew 4) > a ope 
=V,-2eHS. (2.7) 


where N, and Ng are the total numbers of electrons in A and B orbitals, respectively. 
Since N,, Np, Sy, Sp and S are all very large, of the order of the number of atoms 
in a magnetic domain, the distinction between S(S+1) and S?, N,(N,—1) 
and N,”, etc., has been ignored in obtaining eqn. (2.7). Fora given configuration 
which we denote by C(, Ny ; m, Ny), the number of (spin) states with given values 
of S,, Sz, S and S, is obtained by a simple generalization of the Branching Rule 
(Corson 1951, §82). ‘This number, which we shall denote by g(n, S, : m, Sz), is 


alts Sam Ss) =| ( Es )- is be 
oN BF oe sn + Sy n+ S,+1 3m+ Sp eee) | 


The total energy EF of each of the states of given S,, Sp, S and S, is obtained 
by adding to V the sum of the orbital energies for the configuration 
C(n, Nx; m, Nz). This can be written as 


ZU artar t UUpeepes eel l6 eee (2.9) 
Ar Br 
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where the occupation numbers u,, and up, are 0, 1, 2, according as the orbitals 
Ar and Br are unoccupied, singly or doubly occupied. (In eqn. (2.9) and the 
following we use a Roman letter to denote any orbital, not necessarily one singly 
occupied as formerly.) In this notation to specify C(n, N,; m, Ng) is to specify 
the two sets of numbers (u,,) and (up,). Obviously 
Sis sw, eandh lites Net on) Lob. eee (2.10) 
Ar Br 
for all sets (u,,) and (up,). 

This completes the enumeration of the states and corresponding energies. 
To derive now the thermodynamic properties of the system the most direct way is 
to evaluate the partition function Q, i.e. we must sum g(n, S,; m, Sy) e-“/** over 
all spin states, over all configurations and over all allowed distributions of the N 
electrons between the two groups. Thus 

N N-WN, Npl2 Na2 SatSp S Np Ma 


rae AS Snes: Ll Deore re es 


N,=0 Ng=0 Sg=0 Sa=0 S=|S,4—Spl Sz=0 m=28R n=28,4 C(n, Na; m, Ng) 
° —E/kT 
PUES SIRT CaN ee «gil 2 ale geet Mee tier (2.11) 


As all the thermodynamic quantities depend directly not on Q but on In Q, it is 
justifiable to replace O by 
Np Na 
O =5 2 ~ ns US Sys (Oto) Cabrales Gonone (2.12) 
m=2Sp n=2S,4 C(n, Nag; m, Ng) 
together with the conditions that this QO’ be made a maximum with respect to 
S,, 5, Sy, Sgand N,. ‘The five equations which result give the (thermodynamic- 
ally) mean values of these quantities. We denote such mean values by a bar. 
Substituting for E we find that 
SS SOR Fp eee Eaten. (2.13) 

In the following we replace S, and S in QO’ by S,+ Sg, maximizing the resulting 
expression with respect to S, and Sg. This is legitimate since there are no 
cross terms such as SiS, or S,S,. ‘The remaining problem, the evaluation of 
Q’, is now possible accurately only if we make additional assumptions about the 
coulomb integrals.J,,, which so far we have left entirely arbitrary. Such restrictions. 
on the J,,are not necessary for the equivalence of the two free energies calculated 
neglecting and also including spin degeneracy; this always holds under the condi- 
tions we have imposed on the exchange integrals K,,. But to complete the calcu- 
lation accurately we assume there to be only three distinct coulomb integrals 
J,, Jp and J, analogous to the three exchange integrals K,, Kp and Kay. ‘There 
seems in fact little point in working with a more general set of coulomb integrals 
unless we also take a general set of exchange integrals. ‘The function Q’ becomes. 


Q’ =exp a7 (2BH(Ss + Sp) Vo) 


Na Nz 1 
ely x x g(n, Sa; n, Sp) ex - gr (2 Uartart 3 tnren,) 
n=2S, m=2S, O(n, N4; m, Ng) 
oor (2.14) 
where Vo =4N 4? (254 —K,)+ 1N,? (27, — Kp) +4NiNp (2J,3—Kap) 
AE NSIS OT eee (2.15) 


is independent of the configuration of the system so long as N, electrons are in A. 
orbitals and N, in B orbitals. 
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Consider the contribution to Q’ from the term fae e lee eu in 
g(n,S,; m,Sp). The first factor is the number of ways in which }n+S, | 
upward spins and 4n—S, downward spins can be distributed over the » singly 
occupied A orbitals. Similarly the second factor is the number of ways m+ Sz 
upward spins and 4m—S, downward spins can be distributed over the singly 
occupied B orbitals. For all these arrangements © 4,u4,<€4,+ Up,Up,r€pr has the 


same value. We see therefore that 
ee eee n m 1 

> 2: 3 )( Jex {- (Su r€ pete oe )} | 

nN=28, m=28R O(n, Ny: m, Np) o + Sy] \gm+ Sp P kT ai ae Ries, | 

L Soalets)=n(- enone) 

= ex INE | 

= Ca(n, Na) eet c te kT 5 2 } | 


z py af ae > ) 

zs = ge, UBr€ Br ’ 
| 2s, Cp(m, Np) Ce =o 5 P { kT (= Br©B \ 
may be ieee as 


2, a a kT 7(2Pa car) 2, van ie kT 7(2 an car) 
1 1 
s =, exp{ — kT (Zeon) Le, ees { RE (= ce NH pee 


where the new occupation numbers p,,, 9,,, Pp, and gp, areeitherOor1. Thep’s 
may be regarded as occupation ‘numbers for parallel spins and the q’s for anti- 
parallel spins. ‘The outer summations in the four factors of eqn. (2.16) are over 
all sets of occupation numbers satisfying 


2 Pick SNA +-Sas) Bap =a was 
rn 


Epa BN ptSp, 24pr=3Np—Sp- 
B 


Each factor in (2.16) is in fact a simple Fermi—Dirac partition function ; if we 
denote these by Z4(v) and Z,(v) for v particles distributed over the A orbitals only 
and B orbitals only, then (2.16) is just 
+ Za(2Nq + Sa)4a(2Na — Sa)Zn(2Np + Sp)Zp(3-Np — Sp) - 

Dealing with the contributions to Q’ of the other three terms in g(7, S,; m, Sg) in 
the same way we get 

—Z4(3Na + Siy)Za(3Na — Sa) Zp(3Np + Spt 1) Zp (3-3 — Sp —-1) 

—Zy(3Na+Sg +1) Za(3Nq — Sy — 1) Z3(2Np + Sp) Zp(3.Np — ts) 


ZA(ENq + Sy +1) Z,(3-Na — Sg — 1) Zp(3.Np + Sp t+ 1) Z—(3-Np — Sp 1). 
The sum of these four terms can be simplified by using the eee i 
Z(v) , 
Zo+1) = Aly) ON eee (2.17) 


where A, the usual parameter arising in the Fermi—Dirac problem, is determined 
by the number of particles v: 


= 3 [\hexp (ee el (2.18) 
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In this way we find that 


. Nal Abe 
InQ =2(S, + Sp)En = ET 
+n {Z4($Na + Sy) Z,($Ng — Sa) Z3(3.Np + Sp) Zp(3-Np — Sp)} 
Ap (4Np—-S if A, 4N,-—S i 
Ind 1 — —B42— 8 = B/ Se Bet NS eal 2, a eG, Ot . 
et aR eRe aN ane 


The last term is entirely negligible except when S, and Sp tend to zero, as at the 
Curie point. Even then, however, it is only of the order of In N and so may be 
neglected since InQ must be proportional to N. The remainder of eqn. (2.19) 
is just the partition function we would obtain if we neglected spin degeneracy ; 
2S, would then be the difference between the numbers of parallel and antiparallel 
spins in A orbitals, with a similar interpretation for 2S,. 

It only remains to substitute for Z, and Z, to obtain an expression for the free 
energy F. This may be simplified by using the condition that F be a minimum 
with respect to S, and Sy. We also replace the discrete sums over orbitals by 
integrals by introducing density of states (orbitals) factors m,(e) and mp(e). Then 


F=Nana ae Nate “t iN,2 (27, —Ky)+ N,? (2J;, — Kp) 
+3NiNp3 (2J43—Kap)—Sa?Ky —Sp?Kp—2S,Sp Kap 


4 { ms(oin{t bexpgplna—<+ Ka +KanSn + BH) de 

+] paalain {1 + exp apna LANTOS BH} de 

- J male) In {1 + exp apis eK Soe KyySy+BH)} de 
a i np(€) In +exp arte —<— KS ~ KanSs—PH)} Ae. 


with the following five equations connecting 7,4, 7p; Simon Ncand Neecre 


na — 1p = #2 p— Kp) Np — (24 — Ka) Na +4(2Iap—Kap)(Na— Np), 
Peas (2.22) 


These equations can be applied to a system with two energy bands which may 
or may not overlap one another for a certain range of «; they can be worked out 
fully for simple forms of m,(e) and np(c) by using the numerical tables of McDougall 
and Stoner (1938) and of Rhodes (1950). We shall not attempt here the detailed 
application of these equations, but append several general comments prompted by 
the above calculation. 
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§3. GENERAL REMARKS 


(i) On the Density of States Functions 


The functions m,(¢) and mp(e) which we have introduced are not identical with 
the density of states functions usually constructed by energy band calculations, 
The energy « above is the sum of kinetic energy and energy of interaction of the 
orbital function :4(r) with the lattice; the energies obtained in a band calculation 
(for example by the cellular method) are approximations to the energy parameters 
in the Hartree-Fock equations, and these, as is well known, include a contribution 
from electronic coulomb and exchange interactions. This is probably no 
practical disadvantage since, as already mentioned, eqns. (2.20) etc. can only be 

worked out for certain assumed forms of m,(e) and mp(e), and it seems likely that 
we can obtain in this way an equally good representation of the density of states in 
our energy « as of the density of states in the < of the energy band calculations. 


(ii) Application to Iron, Cobalt and Nickel 


The A and B groups may be tentatively identified with the 3d and 4s bands in 
iron, cobalt and nickel, where it is generally assumed some of the 4s orbitals have 
energies « in common with some of the 3d orbitals. (Equations (2.20) to (2.23) do 
not necessarily lead to ferromagnetism, however, and we could equally well identify 
A and B with the 4d and 5s bands of the paramagnetics platinum and palladium.) 
If we can neglect the B orbitals as of lesser importance than the A, then the equations 
for the magnetization are an obvious generalization of those used by Stoner 
(1938, eqns. (2.10), (2.11)). If M bethe magnetic moment then 


M = 1 = =i 
FH =25s= | male) | exp gple—na—KSy—PH) +1 | de 
1 = = 

-| nale)| exp gle—na+KxSa+8H)+1 | dé, Sheer By 
A 

= at 

with n={ nale)| exp gz(e—70—KaSs—BH) +1 | de 
S a di 


1 = =f 
+{ n ,(e) | exp gp (e—70+ KS +8H)+1 | Heo gon ee (3.2) 
A 


These equations have been worked out by Stoner for ,4(e) proportional to «4? and 
by Wohlfarth (1951) for m4(e) independent of «. In both cases N,K,/2 is to be 
identified with their parameter k0’; no other interpretation of R6’ has so far been 
proposed. 

Equations (3.1) and (3.2) for a system with only one distinct exchange integral 
K and one coulomb integral J should also have been obtained by Hirone (1938), for 
this was precisely the model he assumed. In his calculations of the partition 
function, however, Hirone omitted the sum over all configurations of a given number 
of singly occupied levels, and thus obtained an incorrect result. * 


* Dr. H. Watanabe in a recent private communication to the author has said that Hirone and: 
Miyahara subsequently corrected this error in a paper entitled “‘ Ferromagnetism in Semi- 
Conductors ’’’ which appeared in Nippon Sugaku-Butsuri_ Gakkai-shi (1943), 7, 92. Since this 
paper is in Japanese and not generally accessible in the West we have considered it worth while 
to refer here to Hirone’s earlier omission. 


Ferromagnetism 893 


(iit) More General Models 


The particular model that we have examined is the most general for which the 
calculated properties are the same whether we include or neglect spin degeneracy. 
Neglecting spin degeneracy we would say that the energy is entirely determined by 
the two quantum numbers (S,), and (S,), ; including spin degeneracy the energy 
is specified by the two total spin quantum numbers S, and S,. However, for a 
more complicated model, for example one with three groups of orbitals, this is no 
longer true. Neglecting spin degeneracy, the energy would be completely 
determined by (S,),4, (S,)p, (S,)o etc., although including spin degeneracy the 
energy is not now determined uniquely by S,, Sp, Sg ete. 

For these more general systems the equivalence of the results obtained by 
neglecting and including spin degeneracy is restored only if we make the mean 
value approximation familiar in the simple Heisenberg theory. (We note that 
such an approximation in a general model with non-localized orbitals does not 
lead directly to Stoner’s theory in the way that it should do if this were a ‘zeroth- 
order approximation’ as has been suggested by Van Vleck (1945).) In developing 
more accurate theories it seems likely that the correct expressions for the eigen- 
energies are required. ‘This will apply not only in the theory of ferromagnetism; 
taking account of the spin degeneracy will modify the results whenever the average 
number of singly occupied levels 7 is different from zero, and for this it is not 
necessary that S be different from zero; i.e. the equilibrium spin state could be a 
singlet. For example, in the calculation of the influence of the exchange energy on 
the specific heat of free electrons (Wohlfarth 1950, Lidiard 1951) 7 is not zero, 
although so far no account has been taken of spin degeneracy. It is possible that 
inclusion of spin degeneracy would remove some of the notable disagreement 
between this calculation and experiment. Work to examine this possibility is in 
progress. 
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ABSTRACT. Proton magnetic resonance widths and relaxation times have been measured 
in ammonium chloride, bromide and iodide, methyl alcohol and dimethyl sulphate from 
90° kK to 290° x. Line width transitions were observed in NH,Cl and NH,Br at 135° k 
and 99° respectively. There were no sharp changes in either the relaxation time or 
line width at the temperatures of the A-specific heat anomalies of the ammonium salts. 
Evidently, there is no marked increase in the rate of rotation of the ammonium ion at these 
transitions. The results for NH,Br suggest that rotation of the NH,* ion is hindered by a 
potential barrier of about 3 500 cal/mole. The N-H distance is 1:03+0-01 A. The line 
width and relaxation time show pronounced changes at the A-point of methyl alcohol, 
indicating a sudden change in the motion of the CH;OH molecule. In solid dimethyl 
sulphate, a narrow line of small amplitude appears at about 170° K. 


§1. INTRODUCTION: NUCLEAR RESONANCE AND MOLECULAR 
MOTION 

N atomic nucleus of spin s and magnetic moment yp in a magnetic field H 

possesses (2s+1) energy states separated by energy differences uAZ/s. 

In molecular beam experiments, Rabi, Zacharias, Millman and Kusch 

(1938) demonstrated that transitions between adjacent states may be induced 

by an oscillatory magnetic field perpendicular to the main field H and having a 

frequency v given by hy =H /s or 2nv = yH where y =1/(sh/27) the gyromagnetic 
ratio for the nucleus. 

When Purcell; Torrey and Pound (1946) and Bloch, Hansen and Packard 
(1946) succeeded in detecting these transitions in bulk matter electromagnetically, 
it was found that while the resonance field H for a given frequency-is determined 
by the gyromagnetic ratio of the particular nucleus investigated, the shape and 
intensity of the resonance line are determined by the molecular conditions in 
the substance chosen. Since each nucleus is subject to a magnetic field due 
to its neighbours, resonance at a given frequency is observed as a more or less 
broad line spread over a range of magnetic fields corresponding to the distribution 
of the internal fields. In diamagnetic solids the only important contribution to 
the nuclear magnetic resonance line width comes from magnetic interaction 
between neighbouring nuclei. When the nuclei are stationary a formula for the 
second moment or mean square width of the line has been given by Van Vleck 
(1948). Gutowsky, Kistiakowsky, Pake and Purcell (1949) have shown that in 
many cases widths and shapes of nuclear resonance lines may be used to deduce 
internuclear distances in the solid state. 

Molecular motion causes the internuclear fields to fluctuate, and if many 
fluctuations occur during the lifetime of the spin states, it is the average value 
of the internal fields that determines the line width. As these lifetimes are often 
greater than 10-° second, even comparatively slow molecular motion may 
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considerably reduce the line widths (Bloembergen, Purcell and Pound 1948). 
Bitter suggested that observations on the line width of nuclear resonance in 
solids would give information on molecular motion (Bitter, Alpert, Poss, Lehr 
and Lin 1947). In particular, it was thought that definite information could be 
obtained with regard to certain A-anomalies in the variation of the specific heat 
with temperature in substances such as methane, methyl alcohol, the hydrogen 
halides and ammonium salts. If the A-point marked the onset of molecular 
rotation (Pauling 1930) we should expect a very marked narrowing of the proton 
magnetic resonance as the temperature is raised through the A-point. If on the 
other hand, it is due to an order—disorder transition (Frenkel 1935) the magnetic 
resonance would be little affected. 

In the absence of electronic paramagnetism, the internal fields acting on a 
given nucleus may be classified under two heads: 

(a) The field due to the magnetic moments of nuclei in the same molecule 
as the nucleus considered, and sometimes to the molecular magnetic moment 
of this molecule; this field is generally the major part of the internal field, since 
the separation of nuclei in the same molecule is less than that of nuclei in adjacent 
molecules. If molecular rotation, either free or hindered, or rotational oscillation 
causes this field to fluctuate sufficiently rapidly, it is the average value over the 
motion which is effective. Often, as in the ammonium halides, this average is 
zero. 

(6) The field due to the magnetic moments of nuclei in neighbouring molecules, 
and sometimes to their molecular magnetic moments. When relative motion of 
molecules occurs, this field fluctuates with a period of the order of the time for 
the radius vector joining two molecules to turn through a large angle, say 90°. 
In all liquids with normal viscosities, molecular translational motion is sufficiently 
rapid for this class of internal field to average to zero as far as radio-frequency 
resonance is concerned. In solids molecular rotation will reduce this field by 
an amount dependent on the ratio of the intermolecular distance to the radius 
of rotation. It is calculated that in ammonium bromide, considered below, in 
which the separation of adjacent ions is 4-0 A, rotation of the protons of the NH, ion 
about the nitrogen with a radius of 1-0A will reduce this intermolecular field 
from an r.m.s. value of 2:3 gauss to 1-9 gauss. Andrew (1950) has calculated 
that in the normal paraffins molecular rotation will reduce the r.m.s. intermolecular 
field from 2-8 gauss to 1-6 gauss. ‘These calculations are in good agreement 
with experiment, and may be taken as typical. In some solids, however, the 
nuclear resonance lines are narrower than the limit of measurement set by 
inhomogeneities of the applied field (say 0-1 gauss).. This is not compatible with 
the supposition that the only molecular motion is independent rotational motion 
about fixed centres. Gutowsky (1951) and Norberg and Schlichter (1951) have 
shown that self-diffusion accounts satisfactorily for the very narrow nuclear 
resonance lines observed in sodium metal. As the temperature is raised the 
diffusion frequency increases until it produces fluctuations of the internuclear 
fields during the lifetime of the spin states, reducing the effective interaction field 
to zero. It may be necessary to consider such a molecular motion in the other 
cases in which very narrow lines are observed in solids. 

A second way in which molecular motion affects nuclear resonance is through 
the spin-lattice relaxation time. ‘The nuclear magnetic moment of the specimen 
in a given magnetic field is inversely proportional to the temperature of the 
nuclear spin system (Curie’s law of paramagnetism). Application of an 
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oscillatory field near the resonance frequency results in the absorption of energy 
from the field into the spin system. If this energy cannot be transferred 
sufficiently rapidly from the spin system to the molecular lattice, the temperature 
of the spin system will rise and the nuclear magnetization will be reduced. ‘This 
effect shows itself as a saturation of the absorption in high intensities of 
oscillatory field. From measurements of this saturation effect the time constant 
for the transfer of energy between the spins and the lattice—the spin-lattice 
relaxation time—can be determined. When molecular motion takes place the 
important mechanism for the transfer of energy from the spins to the lattice 1s 
by means of transitions induced by components of the fluctuating internuclear 
fields at frequencies near the resonance value. Bloembergen, Purcell and 
Pound (1948) have given a theory relating the relaxation times in liquids to the 
Debye time constant for molecular motion and hence to the viscosity of the 
liquid. In solids in which molecular motion occurs, a somewhat similar 
relation can be expected to hold between the spin-lattice relaxation time and 
the period of molecular motion. In the experiments below an ammonium 
bromide, relaxation time measurements have been applied to the estimation of 
the frequency of molecular rotation at temperatures between 90°K and 290°k. 


§2. APPARATUS 


Except for one or two experiments, all measurements have been made at a 
frequency of 8-155 Mc/s provided by a crystal-controlled oscillator. For protons 
y =2:675 x 10% gauss"! sec-!, giving a resonance field of 1915 gauss. The 
radio-frequency circuit used resembles that of Bloch, Hansen and Packard (1946) 
in that resonance is observed by the voltage induced in a secondary coil 
surrounding the specimen instead of by the change in voltage across the primary 
coil providing the oscillatory field. By disposing the two coils at right angles 
the coupling between them may be made so small that no special balancing 
circuits are necessary to prevent saturation of the detecting system by direct 
pick-up from the oscillator. ‘To distinguish the resonance signal from the 
background, an a.c. modulation method is used (Bloembergen, Purcell and 
Pound 1948). An audio-frequency field whose amplitude H, is a small fraction 
of the width of the resonance line is superimposed on the main steady field H_ 
by means of auxiliary sweep coils round the magnet pole-pieces. In a region 
where the resonance response varies with the magnetic field, for example on the 
flanks of a resonance absorption curve, this small variation of magnetic field 
causes the output of the detector to vary in the same way. This variation of 
detector output is amplified by a tuned audio-frequency amplifier, and fed to a. 
phase-sensitive ‘lock-in’ amplifier (Dicke 1946) in which it is mixed at a second 
detector with a voltage of the same frequency obtained directly from the supply 
energizing the sweep coils. ‘The d.c. output from the second detector is read on 
a meter. Provided the resonance curve may be considered straight over the 
region covered by the sweep (and in the absence of saturation effects) a graph of 
the output meter reading against magnetic field is the differential of the resonance 
curve. ‘I'he method is highly sensitive and is not affected by the frequency — 
response of the amplifiers. 

The main magnetic field was produced by an electromagnet energized by a 
battery of accumulators. Its poles were shaped to give a uniform field at the: 
specimen, which had a volume of about 1 cm’. The small variations of magnet. 
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current required for the plotting of a resonance curve were measured by a 
potentiometer method. ‘To calibrate the potentiometer system in terms of 
magnetic field, use was made of the !"F resonance which at the frequency used 
is separated from the 1H resonance by 125 gauss. The !°F resonance was provided 
by the coil formers, which were made of P.T.F-.E. (polytetrafluorethylene)* to 
avoid the possibility of proton resonance from the formers. 

In order to make relaxation measurements, and also to allow time for thermal 
equilibrium to be established in the specimen in regions of high specific heat, 
it was desirable to keep the temperature of the specimen constant within 1°c 
over periods of about 1 hour. Adequate control of the temperature was secured 
by an apparatus operating on the ‘heat leak’ principle (Alpert 1949). The 
specimen and coils were surrounded by a copper can 2 mm thick which was 
screwed to a copper block connected by a heat conducting rod to a bath of liquid 
oxygen. A heating coil was wound non-inductively on the copper block, and the 
whole system was surrounded by a Dewar vessel. By varying the power supplied 
to the heating coil the temperature of the block and can could be held at any 
temperature between room temperature and liquid oxygen temperature. ‘The 
equilibrium time was about half an hour, and about 10 watts heating was 
required to maintain the highest temperature used, 290°k. The temperature of 
the specimen was measured by a copper-—constantan thermocouple having one 
junction mounted close to the specimen and the other in melting ice. For 
checking and control purposes other thermojunctions were fixed to the copper 
block and at the bottom of the copper can. 


§3. RESULTS 
(i) The Ammonium Halides 


The ammonium halides were investigated in view of the possibility of changes 
in the proton magnetic resonance associated with the A-anomalies in their specific 
heats, discovered by Simon, von Simson and Ruhemann (1927). Our results 
are in substantial agreement with those of Gutowsky and Pake (1948) and 
Sachs and Turner (see Purcell 1951). 

The structure of the ammonium halides is cubic, each cube of eight halogens 
having a nitrogen atom at its centre, with the four hydrogen atoms disposed 
about it in a tetrahedron. Pauling (1930) suggested that the anomaly in the 
specific heat of each salt is due to the onset of rotation of the tetrahedron of 
hydrogens about the nitrogen. Such a rotation would reduce the width of the 
proton resonance line. Frenkel (1935) on the other hand suggested that the 
tetrahedra are not capable of free rotation, but of oscillation about the nitrogen 
atom, the axes of oscillation being ordered in space below the A point, 
while above it their directions are disordered. Such a change would have no 
marked effect on the nuclear resonance. Frenkel’s theory was supported by the 
experiments of Lawson (1940) who from measurements of the elastic constants 
of ammonium chloride was able to derive the value of the specific heat at constant 
volume from the results of Simon, von Simson and Ruhemann. He found that 
the specific heat both above and below the transition is approximately 
18 cal/deg mole, consistent with the attribution of 12 cal/deg mole to lattice 
vibration and 6 cal/deg mole to three oscillational degrees of freedom of the 
hydrogen tetrahedron. Measurements were made on ammonium chloride, 
ammonium bromide and ammonium iodide. In no case was any discontinuity 

* This material was kindly supplied by Messrs I.C.I. Plastics. 
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found in either the width of the nuclear resonance line or the relaxation time at 
the A-point, but line width transitions were observed at very much lower 
temperatures. The line width in ammonium bromide fell from 25 gauss to 
5-0 gauss as the temperature rose above 99° k, and in the chloride from 25 gauss 
to 5-3 gauss above 128°. (These temperatures are somewhat below those given 
by Gutowsky and Pake (1948), viz. 108°k and 135°xK respectively.) In_ 
ammonium iodide the line showed some sign of broadening at the lowest 
temperature reached, 78°K; the transition is probable at about 73°x. A graph 
of line width, defined as the separation of the maximum and minimum of the 
differentiated absorption curve, for the bromide is given in fig. 1. 
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Fig. 1. Line width of proton magnetic resonance Fig. 2. Spin-—lattice relaxation time for proton 
in ammonium bromide plotted against magnetic resonance in ammonium bromide 
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To obtain quantitative information on the nuclear magnetic interaction it is 
necessary to determine the shape of the whole absorption curve, from which 
the mean square width of the resonance line can be calculated. ‘The mean 
square width of the proton resonance from ammonium bromide at 90°K is 
52 gauss”, consistent with the ammonium ions being stationary. Assuming a 
lattice spacing of 4-0 A between like ions (Bartlett and Langmuir 1921) we deduce 
a N-H distance of 1:03+0-01A. This may be compared with the values 
1-025+0-005A from nuclear resonance (Gutowsky, Kistiakowsky, Pake and 
Purcell 1949) and the values from neutron diffraction measurements on deuterated 
ammonium salts, 1-03 +0-02A for the chloride (Goldschmidt and Hurst 1951) 
and 0-99 + 0-02 for the bromide (Levy and Peterson 1951). Above the transition, 
at 137°, the mean square width in ammonium bromide is 3-6 gauss?, indicating 
that rotational motion of the ammonium ion is sufficiently rapid to render 
interactions between protons in the same ion ineffective. A mean square width 
of 3-2 gauss? is given by an approximate calculation of the mean square interaction 
field of protons in neighbouring ions, averaged over all possible positions, made 
by assuming that the protons move independently over spherical surfaces centred 
on the nitrogen atoms. 

A frequency of rotation sufficient to narrow the nuclear resonance may yet 
have no influence on the thermal properties of a substance. An estimate of the 
rate of rotation may be obtained from observations of the relaxation time 7). 
Figure 2 shows the variation of the relaxation time of ammonium bromide 
with temperature. Assuming that the relaxation time is governed by molecular 
motion, the following relation was obtained by Bloembergen (1948) between the 
relaxation time 7, and the characteristic time 7, of molecular motion: 


BS =iK ie 276 
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In this relation, w/27 is the frequency of measurement and the constant K, 
depends on internuclear distances and nuclear moments. We may determine 
K, from the fact that the minimum in the curve of 7, against temperature must 
correspond to the minimum in the relation between 7, and 7,. Values of 7, 
obtained between temperatures of 90°K and 290°x then fit approximately the 
relation 1/7,=2-5 x 10% exp(—1800/T). This type of formula fits the 
supposition that the rotation of the ammonium ion is hindered (Gutowsky and 
Pake 1950), an energy « being required to cross the barrier to rotation. The 
probability of the ion crossing the barrier is then proportional to exp (—«/RT). 
A characteristic temperature of 1 800° k corresponds to a barrier of 3 500 cal/mole. 

Measurements were also made of the variation of relaxation time with 
temperature for ammonium iodide. Our results on this salt are less complete, 
but the general behaviour is similar to that of the bromide. The minimum 
relaxation time of 0-05 sec occurs at approximately 100°x. The barrier to 
rotation could only be calculated very approximately at 1700 to 2200 cal/mole. 
Discontinuities in the line width and relaxation time were found at the first-order 
phase transition of this salt, which occurs at 255°xk, the line width increasing 
from 4-5 to 5-8 gauss and the relaxation time increasing from 1-3 to approximately 
20 seconds. Measurements of relaxation time in the chloride were made only 
at temperatures near the A-point (243°K). They showed that the relaxation 
time is increasing with rise of temperature in this region, and that there is no 
discontinuity at the A-point. All the evidence is then consistent with the view 
that in these salts a hindered rotation occurs at a rate increasing with temperature, 
but much too low to contribute appreciably to the specific heat at temperatures 
of the order of the A-points. 


(ii) Methyl Alcohol 

Methyl alcohol has a A-anomaly in the specific heat at 157-4°k (Kelley 1929) 
and a sharp change in the dielectric constant at about 160°k (Smyth and 
McNeight 1936) which might be connected with the onset of rotation of the 
CH,OH molecules or part of them. 

The results of measurements of the proton line width between 77°K and 
290° are shown in fig. 3. The line width remained constant at about 4-5 gauss 
from 77° K to just below the A-point. Above the A-point the line width was so 
small as to be limited by the inhomogeneities of the magnetic field (0-1 gauss). 
The line width transition took the form not of a steady decrease of the line width 
as the temperature was raised, but of a gradual replacement of the wide line by 
a narrow one. The narrow line first became detectable at 145°k, and became 
more prominent as the temperature was increased until at 156° the wide line 
could no longer be detected. ‘The results for the width of the wide component 
just below the A-point are rather scattered, and there is probably no significance 
in the apparent slight increase of line width. Line width measurements in this 
region of temperature were difficult owing to the long relaxation time, which 
necessitated the use of very small radio-frequency fields. 

To obtain more information on. the progress of the transition as the 
temperature was raised, some measurements were made of the variation with 
temperature of the amplitude of the narrow component. ‘The results shown 
in fig. 4 were obtained using a radio-frequency field of amplitude 2H, = 0-008 gauss, 
and a sweep field H,=0-48 gauss. Under these conditions the line-width of the 
narrow component was entirely determined by the sweep field, and H, was small 


goo A. H. Cooke and L. E. Drain 


enough to avoid saturation effects. The variation of spin—lattice relaxation time 
with temperature is shown in fig. 5. The results between 110° and 157° are 
rather inaccurate, as the saturation method did not yield well-defined values for 
the relaxation time. 

The line width of 4-5 gauss below the A-point is too small to be accounted for 
on the assumption that the methyl alcohol molecules are completely stationary. 
The CH, group at least must have some freedom of rotation, for the assumption 
of stationary CH, groups would yield a line of second moment at least 22 gauss’, 
whereas measurements of this quantity at 90°K gave values between 5 and 
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10 gauss*. (The inaccuracy in the experimental value arises from the shape of 
the resonance curves, which tail off very slowly at each side.) The relaxation 
process below the A-point is presumably largely due to this motion. 

The complete transformation of the character of the proton resonance at the 
A-point indicates a considerable change in the rate of molecular motion in solid 
methyl alcohol at this temperature. ‘The extremely narrow line present between 
the A-point and the melting point could not be accounted for by independent 
rotation of the molecules about fixed centres. It is remarkable that the 
(relaxation time-temperature) curve is continuous through the melting point, 
indicating that the molecular motion responsible for nuclear spin relaxation in 
the solid just below the melting point is little different from that in the liquid. 
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Similar conclusions were reached for solid hydrogen (Hatton and Rollin 1949), 
and for methane (Thomas, Alpert and Torrey 1950), 

It may be noted that our results do not agree with those of Gutowsky and 
Pake (1950) who obtained a line width of 8-5 gauss in the solid state with no 
change at the A-point. It is difficult to account for this discrepancy. The 
temperature control is believed to have been better in the present experiments 
than in those of Gutowsky and Pake, in which there may not have been time for 
the transition to take place, but this would not account for the difference in line 
width below the A-point. Gutowsky and Pake found that the line width is 
considerably affected by impurities. The chief difficulty is removing water. 
The specimen we used, which was kindly supplied to us by Dr. D. V. Clifton 
of the Dyson Perrins Laboratory, Oxford, had been dried by the method of 
Bjerrum and Zechmeister (1923). The melting point was 176-:0+0-3°K. 
Moreover, a subsidiary experiment made on a specimen containing a few per cent 
of water gave exactly the same line width at 90°x. Later measurements have 
shown that the measurement of line width in terms of the separation of the 
maximum and minimum of the differentiated absorption curves is not entirely 
satisfactory, the line width so measured increasing as the amplitude of the 
radio-frequency field is decreased. 


(iii) Dimethyl Sulphate 

Dimethyl sulphate shows an interesting transition in the dielectric constant 
at 203°xK (Smyth and Hitchcock 1932), suggesting some freedom of molecular 
movement above this temperature. The melting point is 242°K leaving a 
considerable temperature range for observations between the transition and 
melting points. 

The measured variation of line width with temperature is shown in fig. 6. 
It remains constant at 3-4 gauss up to about 170°k, when it decreases quite 
sharply to zero. Results in the transition region were extremely difficult to 
obtain owing to time or hysteresis effects. Above 175°K the resonance line 
consisted of two parts, a very narrow component and a component of the same 
width as below the transition. The variation with temperature of the amplitude 
of the narrow component is also shown in fig. 6, for the conditions H, =0-008 gauss, 
H,=0-5 gauss. There is no change at the dielectric transition point (203°K), 
the amplitude remaining constant from 190°K to the melting point, where it 
increases by a factor of ten. 

Relaxation measurements were not attempted except in the region below 
the transition. The relaxation time decreased as the temperature fell (0-6 sec at 
125°k, 0-14 sec at 77°). 

The line width below the transition is small enough to suggest that at least 
the methyl groups are rotating. Above the transition a more violent motion 
must occur, as in methyl alcohol. 


$4. CONCLUSION 

The variation of proton magnetic resonance line width has now been 

investigated in a number of materials. It appears that the rotation of the 

molecules in the solid state does not usually begin at a A-point or first-order 

transition, but that the rate of rotation generally increases steadily with 

temperature. Thus many substances are potentially interesting with regard to 
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their nuclear resonance behaviour, and not merely those which show specific 
heat anomalies. In some cases it has proved possible to relate the variation 
of the rate of rotation to a potential barrier. The correlation between evidence 
from nuclear line widths and from observations of dielectric constants is poor, 
The proton lines observed in many organic substances at the temperature of 
liquid oxygen are narrow enough to indicate molecular rotation at frequencies 
greater than 100 kc/s, whereas dielectric constant measurements do not. 

The very narrow lines observed in many substances present a difficulty of 
interpretation. They can occur only if the fluctuating magnetic field at one 
nucleus due to all its neighbours averages very nearly to zero. It is at least 
required that the field changes sign as a result of the motion. If neighbouring 
molecules are moving independently, the radii joining pairs of protons in 
neighbouring molecules must turn through angles of at least 90°, which is not 
in general possible as a result of molecular rotation about fixed centres. It may 
be necessary to consider some other relative motion of molecules, such as 
self-diffusion. 
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ABSTRAC T. For an isotropic two-band conductor in a transverse magnetic field, the 
various galvanomagnetic, thermomagnetic and thermoelectric coefficients are derived by 
elementary arguments in terms of those of the individual bands. The expressions are applic- 
able to quasi-free electrons of any degree of degeneracy : for fully degenerate electrons they 
reduce to those already given by Sondheimer and Wilson. When the square of the thermo- 
electric power is comparable with the Lorenz number, as may happen in semiconductors, 
the expressions for the adiabatic thermomagnetic and thermoelectric effects become 
extremely complicated, though those for the isothermal effects remain reasonably simple. 


§1. INTRODUCTION 


N extending the electron theory of metals and semiconductors beyond the 
] point represented by the simple free-electron model, a number of authors 

have studied the properties of a model consisting of two non-interacting 
bands of electrons. In particular, the behaviour of such a model in a transverse 
magnetic field has been studied in some detail, because it is qualitatively different 
from that of the simpler model, and. because it represents more closely the 
behaviour observed in real conductors. We may mention the work of Riecke 
(1898), Fowler (1936), Jones (1936), Kohler (1949 a—d), and Klahr (1951), and 
in particular of Sondheimer and Wilson (1947) and Sondheimer (1948), hereafter 
referred to as S-W and S respectively. ‘These last authors have considered a 
model consisting of two isotropic bands of degenerate quasi-free electrons and 
have given expressions which are approximately valid at all temperatures, including 
the region of the Debye temperature, for the electrical and thermal conductivities, 
and the various thermoelectric, thermomagnetic and galvanomagnetic effects, 
in a transverse magnetic field. ‘These results were obtained by first finding 
‘interpolation formulae’ for the properties of a single band in a transverse field, 
which reduced to the known correct expressions for temperatures well above and 
well below the Debye temperature, and then evaluating the properties of the 
two-band model in terms of these formulae. 

It is clear from the nature of the two-band model that it must be possible to 
express its properties quite generally in terms of those of the individual bands, 
since for given fields, the currents flowing in a two-band conductor are simply the 
vector sums of those produced by these fields in the individual bands. For 
degenerate quasi-free electrons, the results are implicit in the work of Sondheimer 
and Wilson, but their simple nature is somewhat obscured by the mathematical 
details of the ‘interpolation formulae’. In the present paper we give a brief 
derivation of the general expressions, applicable to any isotropic two-band 
conductor of any degree of degeneracy. In deriving these expressions it is 
convenient to make use of certain relations which are considered in the next 
section. 
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§2. DEFINITIONS. THE HEURLINGER RELATIONS 


Consider an isotropic conductor, in a uniform magnetic field H directed along 
the x axis, in which is flowing a longitudinal current, electrical or thermal, of 
density J, or W,. Then there are four main experimental arrangements to 
consider. In all of them, the transverse electrical current is zero: J,=0. In 
two of them the primary current is electrical, andin two thermal. In one electrical 
and in one thermal arrangement, the transverse temperature gradient is zero: 
oT/éy=T,' =0 (isothermal arrangements); in the other cases, the transverse 
heat flow is zero: W, =O (adiabatic arrangements). ‘These four arrangements, 
and the parameters defined by them, are classified below. In the vector diagrams, 
the various secondary fields and currents are shown having the signs they will 
have for an ideal free-electron metal, according to Kohler (1941 a), if the field H, 
is positive, i.e. directed out of the paper towards the reader. ‘Thus for type I, 
E,, is shown as negative, because the Hall coefficient Ay is negative for a free- 


electron metal. 


Isothermal Arrangements I,=T,' =0 


Type lame ten) 


v 
E, : 
L JE Sos E,|/HI,=A'y, W/L, =-'w- 


Lepealle al. W, E, 
i —=>W , t +E, 
T ,.<——_ 
= W,/ 1b =K', =— i el Oe = Bay, OWES =p'p. 
Adiabatic Arrangements I,=W,=0 
Ly poll aad =) ee Wage 
ee i Be 
E, 
LE = 07s E,/HL,=A*qg, —T,'/HI,=Apz, Wf]: = we 
Le 
Type Va <= (nel ee 
>W, E, 
—W,j/T aks SEATS SB TDS en ee 


‘These four arrangements thus suffice to define the isothermal and adiabatic 
electrical and thermal conductivities, o and «, Hall and Ettinghausen—Nernst 
coefficients Ay and Byy, and ‘'Thomson coefficients’ wy and py, and the Etting- 
hausen and Righi—Leduc coefficients Ay and Bp;. Now it is obvious from the 
vector diagrams that types II and IV differ only by a slight relative rotation of the 
x and y axes, and that type I can be constructed from type III by superposing a 
small field of type II in the y direction. It follows that among the numerous 
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coefhicients defined above, relatively few are independent, and in fact the following 
relations must hold (p =1/c): 


jap Unt Siete is els Cin i ee (1) 

ee at att 0) Vit = Sees Pend Bt (2) 

oe eet ee er Eh btm mg heute ae, (3) 

Bryn = B yy M'pBriy DOING (4) 

Peewee geal | uae (5) 

ace Ae SD eee BO? 5 a ee ee (6) 

also WS CACHIC ke Wika Se uaeet (7) 
W (1) ='B,,HT",. fe (8) 


These relations are necessarily valid for any conductor (two-band as well as 
one-band) which in a field H, remains isotropic in the xy plane. They are thus 
applicable for instance to a conductor of cubic symmetry, if z is a tetrad axis, as. 
well as to free-electron models. ‘They can obviously be generalized immediately 
to anisotropic conductors, provided that H is directed along a principal axis, but 
not so readily otherwise. They were first given by Heurlinger (1915, 1916); 
they have also been derived recently, in the course of an analysis of the problem 
by the methods of Onsager thermodynamics, by Mazur and Prigogine (1951). 
and by Callen (1952). We shall refer to them as the Heurlinger relations. 

The coefficients jz and ww which appear in eqns. (3) to (6) are of less interest 
in practice than the thermoelectric power ¢« and the Peltier coefficient II of the- 
material. For any conductor, these are related to x, and uy by 


[ep =e— = (d¢/aT), Squstaigte (9) 


Hog dd C1 Ci ete (10) 
where —e is the electronic charge and ¢ the chemical potential of an electron,. 
so that we may at once replace (5) and (6) by 


ef=e +B inBatt’, serene (5 a) 
Il?= Tl'+ Bp ApH”. iE kal a veilens (6 a) 


Now consider eqns. (3) and (4). It was pointed out by Heurlinger that A*y and 
B*,zx cannot be measured directly, but that the observations must be corrected 
for the thermoelectric e.m.f. in the leads (cf. also Kohler 1941 c). If, however, 
this e.m.f. is calculated in the usual way as Je dT rather than as Ju,dT, the 
quantities A®, and B*,y so obtained will not be those defined above, but will be 
given by 
Ang=A’n re Ag © 8 alae (3 a) 

Bay=Boon—€ Bey —~—sés tn nn (4a) 
The two alternative definitions of B*,y corresponding to (4) and (4a) have been 
used respectively by Kohler (1941 a,b), whose expressions satisfy eqn. (4), and 
Sondheimer (1948), whose equation (39) corresponds to our equation (4a). 
Since we shall be concerned with deriving generalizations of Sondheimer’s 
expressions, we shall henceforth adopt the definitions of A*y and B*gy corre— 
sponding to eqns. (3a) and (4a). 
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In addition to the above relations, the coefficients II and « are related by the 
second Thomson relation : 


I'= Te; 1 ik—-te Epon Ala NYY ote AS: (11) 
and finally, combining (5a), (6a), and (11), we obtain the Bridgman relation: 
Bowl =“ Ag. ccc ewe (12) 


Equations (11) and (12) have been obtained by Meixner (1939, 1941) and Kohler | 
(1941 c) using the electron theory of metals, and thermodynamically by Callen 
(1948) and Mazur and Prigogine (1951) using the Onsager reciprocal relations. 

In considering the two-band conductor, it is easiest to derive the isothermal __ 
coefficients; the adiabatic coefficients may then be found from them with the 
help of the above relations. 


§3. THE TWO-BAND CONDUCTOR 

From now on we shall usually be concerned with isothermal values of o and 
Ay, and the superscript i will be dropped; where the quantities o*, A*, occur 
they will be so written. 

Consider a conductor containing two isotropic, independent bands of conduc- 
tion electrons, whose individual properties we shall denote by subscripts 1 and 2. 
Then both bands will necessarily be subjected to the same electric fields and 
temperature gradients, but the resultant electrical and thermal currents in the 
individual bands will not in general have the same magnitude, and unless 
Ay 10;=Ay.20. and Bry 1=Baz,2 they will not have the same direction. 
The two-band effects arise simply because the total electrical and thermal currents 
in the conductor are then the vector sums of those in the individual bands, rather 
than the algebraic sums, so that, for instance, if the total transverse current is 
zero, the transverse currents in the individual bands will be not zero but equal and 
opposite. If Ay 4o,=Ay ec, and Bry ,=Bpz»2, so that the angles between 
field and current are the same in both bands, the generalizations introduced by 
‘the two-band model vanish. 

The properties of either band individually, or of the conductor as a whole 
{denoted by subscript 3) may be represented by equations of the form 


I, =F, +oph, tool, +apl,’ 7} 
yh => —Opk, tab, —apl, “aol, 
W= 82+ Pabya Pol atten. yal ie ee (13) 


W,= — Bale t+ Brky— Poly +Boly 
I, W,=f(E,, T,’) only, 
‘and if we write 
Ayo =i 7% Save (14), o/(1+f)=¢d ree ne (15) 


it is easily shown that the following relations hold between the coefficients «, 
and the quantities defined in §2: 


a, =¢, apy = >t, 
to= —¢(u'y— BiyyHt), tp = —P(u'gt + BiyyH), 
By =¢(u'p—*'AgHt), Bp=G(upttApH) == bonuses (16) 


Bos —K'+ bl (wyt +11 ApH) B yy — (uy —«'ApHt)u'g], 
Bp = —' Ba — Ol (u'pt + ApA)u'y t (wy —*' AgHt) B yy]. 
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The properties of the two-band conductor are now given at once, in, terms of 
those of the individual bands, by the set of equations 


Og = Ay + Xe, Pastit pr 6 es Aas (17) 
Moreover, using (9) and (10) in eqn. (17), it is easily seen that the terms in £, d¢/dT’ 
always cancel, so that if we define «’, 8’ from «, 8 by replacing ply and py every- 
where by «', Te’, eqns. (17) become 
Og’ =O,’ +0’, ured Oe eel ay re eres (17 a) 


[It may be noted that for a degenerate free-electron model, S—W, egns. (20) 
and (22), and S, eqn. (2), define quantities V, W, X, Y, P and Q related to our 
a”, Bp’ by 
P-t7 je W——Ca,/e", X= —Ci kT, Y=CPpp/R'T, 
P=Co' (/k?Te = — Cp’, /k?T?e, O = — Cu'p/k®? Te = Ch’,/k®T?e 
where C=372h4,] 


From (16) and (17 a) the parameters of the two-band model at once follow in 
terms of those of the individual bands. For instance, from «’, and «’, we have 


b3=$1 + do, Pstg= ity toto, = vane (19) 
which give, using (14) and (15), the expressions 
ee (0, +09)? +0;°0,°H*( Aq + Ane)? 


a6 o, + 0,+0,0,1H*(o, Aq 1? + 02.44 2”) er ea ere ek, (20) 
and 
a ae 0°? Ay 1+ 092Agq 2+ 0,70,"H? Ay 1 An o(Anit Ang) 21) 
— (9, + 0)? + 0,20,"A(Ayq 1+ Ane)” “fag ( 
which for small H reduces to 
= 0° Aqito2"Ap 2 
ie (0, +02)? 
Similarly, from «’, and «’p (or fh’, and fp) we obtain 
b5(1 + ts”)e's =P[(1 + tyts)e's + (t3— 41) Bay, 1 A) 
+ ol (1+ tots)en + (ts—to)BynoH],  —«.. + (22) 
and 
$3(1 + t;”)B gn 3 =¢,[(1 + t,t) B py 1H + (t,— ts)e', 
+ dol (1+ tots)Biun cH +(to—ts)e'o], see oe (23) 


for the isothermal thermoelectric power and Ettinghausen—Nernst coefhcient of 
the two-band conductor. 

The expressions for f’,, and fp are considerably complicated by the terms in d 
(cf. eqns. (16)): neglecting these terms for the moment, we obtain relatively 
simple expressions for x', and Bay, 3, and from these, using (3), (5a), and (12), 
expressions for the quantities x°;, <*; and Ag, may be found from those for 
ws, es and Biyys. To facilitate comparison with the expressions given by 
Sondheimer (1948) we have expressed the resulting equations (and also eqn. (23) 
for By, 3) in terms of the coefficients o, Ay, x*, Bay, c* and Bgy for the individual 
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bands, eliminating the coefficients ¢', B'yy and Ay by means of the Heurlinger 
relations. Writing for brevity 


=k",(1+ Bar," H*) C, =«*,(1 + Bp”) 
CsSe(t be Baro) Cy =«°(1 + Bar Brief”) 
D, =** (Bar — Bari)t Dz =«*,(Bari— Bar»2)4 
S, =0,2(1 + 4,2) + o,02(1 + t,¢2) S,=052(1 + t,7) +.0,0,(1 +t, 25) 
HV, =0405(t, — ty) HV, =0402(ty—t), 


the expressions for the two-band model are: 
(Cy + Cox) +49 (Co + Cia) 


Be eee btw 24 
ans CLG 2 
BariC, + Bere Ce 5 
Bar3= “< . Et aa ry pimeaas enarae ge (25) 
€*3(S) a SoC, a C2) = [e" (Cy ot Co) — B "an 11 D2|S; 
Fs ie 2 + B* gy, h(C, 4+ CL )IAV, + sae Nessie ce (26) 
Ag3(S1+S2)(C; +0,)=< * Barit Bran s)Si 
+ ae — B py Bai” VAI + ete o> Oe eae ete oe (27) 
B'py,3(Sy + Sa) = eo [(<*: Bait Bra) Si 
= (<*, a? B* py Bay") Val = nets re eee (28) 


where eqns. (26), (27) and (28) are completed by exactly similar terms with the 
suffixes 1 and 2 interchanged. For H=0 eqn. (26) reduces to 
e*, = (oe) + one*a)/ (Gi ey) —- 9 eens (28 a) 

(cf. Fowler 1936, equation (1241)) which can be derived very simply by considering 
the equal and opposite longitudinal currents in the individual bands caused by 
the net excess e.m.f’s (e*,;—¢*,) and (c*,—€",). 

Equations (24) to (27) have been obtained by neglecting the terms in ¢ in the 
eV for 8’, and B’y. If these terms are retained, the following expressions 
for x's and «',Bryz,, result: 


T 
iy =ael, +c + ae ee g, (Ae? — AB*H?)(o1+ 05) 
— ee + Ay») |, © ie), Seve (29) 
To,o. oe 
K3Bxy,3=K Bayi tk eBatot 75. [(Ac? — AB? H?)o,0,(Ay1 + Aga) 


+ZAcAB(o,+05)],,. (30) 


where AB= B'yy ;— B'yy and Ae=e',—e',. The expressions for *s, Brz,s, €*s 
and A» will correspondingly be a great deal more complicated than those given 
by eqns. (24) to (27), which are obtained on neglecting the terms in square brackets 
in (29) and (30). For small magnetic fields, these terms are smaller than the main 
terms in (29) and (30) by factors of order ¢?/L,, where L,, in the Lorenz number, 
equal to 2-24 x 10-4microvolt?deg-*. Thus for metals, where « is of order 
1-10 uv/degree, the correction terms are always negligible (for a degenerate 
free-electron model «?/L, is of order (RT/£)?), but for semimetals and semicon- 
ductors the correction terms may become comparable with the main terms 
(e.g. for bismuth, germanium and silicon, « is of order 100-400 uv/degree), and 
eqns. (24) to (27) must be replaced by the exact expressions derived from (29) 
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and (30). If thermal measurements could be made on semiconductors with 
£,=E,=0, rather than /,=J,=0, these complications would be avoided, © 
but this is probably impracticable. 


§4. DISCUSSION 

With the proviso just mentioned, eqns. (20) to (30) are formally applicable 
to any two-band conductor, of any degree of degeneracy, to which eqns. (13) are 
applicable, including not only models consisting of isotropic quasi-free bands of 
electrons, but also conductors of, for instance, cubic, tetragonal or hexagonal 
symmetry, if H, is parallel to the symmetry axis. The generalization to con- 
ductors of lower symmetry is straightforward, again provided that H, is parallel 
to a principal axis (cf. Jones 1936). If t,=t, and Bg, = Bayz the characteristic 
two-band behaviour vanishes, and eqns. (20) and (21), for instance, reduce to 
G3 =0, +02 and Aq303= Aq 10, = Ayqooy for all H. 

Applied to the model of two quasi-free degenerate bands used by Sondheimer 
and Wilson, and using for the parameters of the individual bands their interpolation 
formulae (S—W, eqn. (54); S, eqns. (37), (38), (41), (42)), our equations (20), 
(21), (24) to (28) become identical with theirs (S—-W, eqns. (57), (74); S, eqns. (22), 
(54) to (57)). Sondheimer (1950 a) has since given the exact evaluation of o for 
a single-band free-electron metal in the intermediate temperature region, and has 
also shown (unpublished) that the exact value of Ay in this region is inappreciably 
different from the value —1/nec. These results may be used in (20), (21) to 
find more accurate values for o and Ay for the two-band model in the intermediate 
temperature region. 

If the parameters of the individual bands vary with H, the values to be used 
in these equations are of course those appropriate to the applied field. For 
degenerate free electrons the variation of o and Ay with field is normally negligible, 
but Sondheimer (1950 b) has found that for a thin film in a transverse field normal 
to the film, curious oscillations of o occur with H. Unfortunately they occur 
only at such high values of H that for any real metal—even for sodium—the 
normal magneto-resistance effect due to anisotropy would tend to mask them. 
Using Sondheimer’s results for o(H) and Ay(H) for a single band in eqns. (20) and 
(21), it would at once be possible to evaluate the behaviour of a thin film of the 
two-band model, and by appropriate choice of the parameters of the model it 
would be possible to give it the same bulk magneto-resistance behaviour in a 
transverse field as, say, sodium; but such a model would be a considerable 
over-simplification (for instance, it would give no longitudinal magneto-resistance 
effect; the two-band effect, by its nature, cannot do this, whereas in sodium the 
longitudinal and transverse effects are comparable in size), and it seems doubtful 
whether it could safely be used for the quantitative evaluation of experimental 
results on thin films. 

For non-degenerate bands, i.e. semiconductors, the values of the parameters 
of the individual bands will in general vary with the applied field (cf. for example 
Harding 1933). In low fields, however, the zero-field values of the parameters 
may be used with, for example, eqn. (21a) to find the zero-field parameters for 
the two-band model (cf. Fowler 1936, eqn. (1247), Klahr 1951, eqn.(3)). Wright 
(1951) has made an extensive theoretical investigation of the zero-field parameters 
of a free-electron conductor of any degree of degeneracy, for various types of mean 
free path behaviour; his results may at once be extended, using the equations of 
§3, to two-band conductors. 
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It may be noted that the properties of an isotropic single band may be deduced 
by an extension of the methods of §3, by regarding it as a continuum of bands, 
each contributing do to the total conductivity. Equations (19), for instance, 
then become 43=| dd, d3t3=Jtd¢, and if we put t= —eH1/mwe for quasi-free 
electrons, and do = —(47e?l/3m)v?(0f,/dv) dv, we obtain the usual equations 
for o and Aj, first derived by Gans, and applied with various assumptions con- 
cerning /(v) and f, by Sommerfeld and Frank (1931), Harding (1933), Klahr (1951) 
and Wright (1951). Here, as before, the effects of a magnetic field vanish if 
t is the same for all electrons or if the field is longitudinal, since the effects again 
arise from the inequality of the Hall effect for the various groups of electrons. 

Lastly, we may consider the effect of recombination on the applicability of 

these equations. Fowler (1936) pointed out that in a transverse magnetic field 
there will be equal and opposite, non-zero transverse currents in the two bands 
when the total transverse current is zero, and that in equilibrium the consequent 
accumulations of charge in the individual bands must be balanced by either 
recombination or back-diffusion (cf. also Welker 1951). Experimentally it is 
found that for semiconductors recombination is rapid (Torrey and Whitmer 1948, 
p. 55) and occurs chiefly at the surface of the conductor (Shockley 1950). In 
dealing with degenerate two-band models it has generally been tacitly assumed 
that although the bands are non-interacting, no charge excesses are built up: 
these assumptions are compatible with the existence of rapid surface recombin- 
ation. For the results of the present paper to be applicable, it is necessary that 
recombination shall be rapid, so that the charge excesses in the individual bands 
remain small enough not to effect their properties: recombination may occur 
either at the surface or in the volume of the material, provided in the latter case 
that, in the presence of this interaction, it is still possible to define unambiguously 
the conduction parameters of the individual bands. 
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ABSTRACT. Samples containing elements of even Z have been examined for possible 
beta-activity using a photographic emulsion technique. For the elements Cr, Fe, Ni, Zn, 
Ge and Cd the results show that the effective half-value period for double beta-decay for 
the elements as a whole is in each case greater than 101” years if the disintegration energy is 
2 Mev or more (or greater than 101° years if the disintegration energy is 0-35 Mev or more). 
For Ca, Mo, Sr, Sn, Te, Ba, W, Os and Pt lower limits between 2 x 1014 and 6 x 101% years 
are similarly obtained on the assumption of 2 mev disintegration energy, the residual 
uncertainty in most cases being due to the slight contamination of the samples with natural 
radio-elements. A possible weak beta-activity shown by samples containing Mo requires 
further examination. 


§1. INTRODUCTION 

HE half-lives predicted for the process of double beta-decay depend 
strongly on the form of decay theory employed. On the basis of the 
Fermir theory the simultaneous emission of two electrons must be 
accompanied by the emission of two neutrinos. ‘Transition probabilities for 
such a disintegration have been investigated by Goeppert-Mayer (1935), who 
showed that they are so low that the process is almost certainly unobservable 
directly, unless the decay energy is improbably high. Furry (1939) calculated 
the disintegration probabilities on the basis of the Majorana theory of beta-decay. 
In this case a double beta-disintegration could occur without the emission of 
neutrinos and with a correspondingly greater probability. ‘The half-lives 
predicted in each of these cases for a nucleus for which Z=31, and for various 

values of the energy available, are quoted in table 1. 


Table 1 
Energy (Mev) 2 3 4 5 6 10 
Ty)2 (years) : Fermi LOMA 10 2 108* 07 108? 8X 1082 TO 
Majorana ee OS IOEY — ASEANOE SCIROCCO 


The experimental evidence available suggests that the probability of double 
beta-decay is considerably lower than that predicted by Majorana theory. 
Lower limits for the half-lives of 12°Te and ?8U of the order of 10!° years have 
been established by analysis for the reaction products (Inghram and Reynolds 
1949, Levine, Ghiorso and Seaborg 1950), and an examination of palladium 
by Winter (1952) has established lower limits for the half-lives of the isotopes 
of even mass number of the order of 101% years. A value for the half-life of 
130T'e of 1-4 x10?! years has also been determined by Inghram and Reynolds 
(1950) from an examination of the isotopic constitution of xenon from old 
tellurium ores. 
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The activity reported by Fireman (1949) for 4Sn, corresponding to a 
half-life for double beta-decay of (6:5 +2) x10!° years, has not been confirmed 
by the recent work of Kalkstein and Libby (1952), who established a minimum 
half-life of 1017 years for this isotope. Fireman has also recently failed to 
confirm any activity (Fireman and Schwarzer 1952). 

This paper reports an investigation of the possible activities of a number of 
elements of even atomic number from which minimum half-lives of the 
individual isotopes of the order of 1017 7 years have been established, where 7 
is the proportion of the isotope concerned in the natural element. 


§2>THE DETECTIONS OF (SMAI SACALVIahi> 


The electron-sensitive photographic plate has been employed for the 
detection of these very low activities, since it permits continuous observation 
of the suspected materials over long periods of time. Since such plates 
record the passage of every ionizing particle which enters the emulsion layer, 
large numbers of tracks, due to cosmic radiation and general radioactive 
contamination, are formed even in the absence of a localized source. However, 
by shielding the plates within steel and lead screens at a depth of 1860 ft. below 
ground, it has been found possible to reduce the background against which 
any activity must be observed from about 200 to 2-5 tracks per. mm? of a 
200 thick emulsion per day (Fremlin and Walters 1950). Also, with the 
co-operation of Messrs. Ilford Ltd. it has been possible to commence exposures 
with plates containing only some 50 tracks per mm? without using eradication 
processes. Such an initial track density is not a serious limitation when. 
exposures of one or two months duration are made. 

It was found convenient to expose thick sources of the substances examined. 
to the surface of the plates rather than to incorporate them in the emulsion. 
layer; this resulted in an increase in the ability to detect weak activities, since: 
it is usually impossible to incorporate a large percentage of foreign material. 
into the emulsion without loss of photographic sensitivity. 

The technique employed was to embed a few grammes of high purity material 
into recesses in graphite blocks, powdered samples being bound with a few 
drops of a dilute solution of collodion in acetone, previous measurements 
having shown these substances to be inactive. The graphite blocks then faced 
the emulsion layer during the exposure. This exposure technique avoided: 
chemical interaction by preventing actual contact between the sample and the 
emulsion, and also, by virtue of its simplicity, enabled a survey of a considerable 
number of elements to be made. Graphite was found convenient as a supporting 
material, since high purity samples could be obtained which showed no. 
detectable activity and since it does not react chemically with the emulsion. 
‘The graphite used in these experiments was obtained from the Atomic Energy 
Research Establishment, Harwell. During the exposures the plates were: 
stored in an atmosphere of nitrogen at a pressure slightly in excess of atmospheric,. 
relative humidity being controlled by the presence of phosphorus pentoxide. 
Under such conditions exposures of four months’ duration were made without 
latent-image fading or loss of photographic sensitivity. 


After development the plates were examined independently by several 


observers and a count made of the electron tracks crossing the upper emulsion 
surface both in the regions exposed to the samples and in those which had faced. 
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the graphite. This resulted in a greater sensitivity for the detection of a 
pure f-emission than a count of all the tracks present throughout the emulsion 
layer; the existence of K-capture or y-activity would not, however, be 
distinguished under these circumstances. Several exposures of different 
duration were made to each sample at intervals of a few months, so that the 
existence of any contamination by short-lived artificial radio-elements would 
be observed. All elements possessing stable isotopes isobaric with nuclei 
differing by two in atomic number were examined in this way, with the exception 
of some rare earth elements and the rare gases. ‘The results of the activity 
measurements obtained in the longest exposure made to each sample are 
quoted below for all those materials which were obtainable in an uncontaminated 
condition. It was considered highly unlikely that any contamination by the 
long-lived natural beta-emitters could be significant. Contamination by uranium, 
thorium or their products, on the other hand, was probable but in all cases 
would be indicated by the appearance of «-particle tracks, and a count of these 
would enable an estimate to be made of the order of magnitude of the /-particle 
output to be expected from the contamination. 

The activity of the graphite supports was examined in some detail and 
found to be of the order of 0-04 electron per mm? of its surface per day, which 
could, in general, be neglected in comparison with the statistical uncertainties 
involved in the measurement of the sample activities. ‘Thus counts in regions 
of the plate exposed to graphite could be regarded as the background against 
which similar counts obtained in areas facing the samples should be compared. 
An investigation into the origin of the background tracks observed in plates 
exposed under these conditions has been made and will be described in a future 
publication. 


§3. EXPERIMENTAL RESULTS 


Table 2 
Count in loaded area Paes Exposure Samp le Min. half- 
Sample Po hae ee ae ground dave) activity life Ty j./ 
; (B/mm*) (6/mm*/day) (years) 

CaCO; a 86:0+4:5 none 75°5+1°6 56 0:19+ 0-07 ZX1038 
Chromium a 66:0+4:-0 4:3+0°6 66:0+ 3:0 2G 0:00+0:15 2:2 10% 
Iron sponge a 76°5+4:3 1:0+ 0-4 1579 + 1:6 56 0:18+0:07 1:3x10" 
Nickel b 52:0+1:6 0:-5+0:1 52°5+1:6 Di 0:00+0°07 3°2<101? 
Zinc a 44°84+1°8 0:3+0-1 A] 5 at 1°38 Diy 0:00+ 0-12 2x10" 
Germanium a 50:0+1:-7 none 52:°5+1:6 De 0-00+0:08 2-8x10"7 
Molybdenum a 87:5+2:1 none 54:14+1:2 27 1:24+0:09 1:5x 1016 
Sree. a 79:14+3-4 none 63°38 115 56 0-27+0:07 3:-4x10% 

Cadmium 99:99 86:0+4-6 4:5+1-0 75°5+1°6 56 0:19+ 0-09 NOE? 
SnO a 94:04+4°8 1:5+0°5 OSeoise 125 56 0:56+0:09 1:4x10' 
124Sn oe ekeise of 20:6+ 1:4 46-44 2°3 15 6-10" 2>a11035 
Tellurium a7 154 540 none 63°83 1-5 56 0-24+0-11 1:3x 10% 
BaCO, By iMG ae Hac ll 47-2+2°-4 15 4:5 +0-4 1:8 x 10% 
‘Tungsten b 53:°5+1°-4 0:95+0:10 39-1+0-9 27 0:54+ 0-06 2x 1028 
Osmium ey Silva PAW 51+10 (5° o 1e6 56 40+3:°6 2°3°< 104 
Platinum b 67:2+3:°3 3:0+1-0 B25 it dic 27 0755450139 17 < 1028 


a: Johnson Matthey spectroscopically standardized; b: purity not known; 
* separated isotope, A.E.R.E. Half-lives for E=2 mev 
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Table 3 gives the expected disintegration energies for a few of the isotopes 
investigated, based both on the empirical mass tables drawn up by Metropolis 
and Reitwiesner (1950) and on experimental mass values where these are 
available. 


Table 3 
Energy difference between 
Isotope % abundance Isobars % abundance isobars (Mev) 
Theor. Exp. 
ON 99-57 Ca 96-92 —3:8 —0-3* 
6Ca 0-0032 wT; 7-95 ig 
18Ca 0-179 eT; 73-45 73 43% 
8Fe 0-33 Ni 67-76 ohm —1-35* 
sNj 1-16 817 n 48-9 —0-9 —1-85* 
20S 33-03 20Te 0-09 0:55 
1225p 4-78 12Te 2-43 2-8 
124Te 4-59 5-1 1540-44 
124 . HE 
28 ape 124Xe 0-094 —0-88 
sire 18-71 TADS 0-088 1-4 
EEM De 31°85 NONE 1-90 325 
13 OSES 4-07 5-6 
ere art site 0-101 —0-06 


* Collins, Nier and Johnson (1951, 1952). + Hogg and Duckworth (1952). 


§4. DISCUSSION OF MEASUREMENTS 


The estimation of the minimum half-life which can be ascribed to any 
isotope as a result of these activity measurements involves the calculation of 
the specific activity of the sample corresponding to the observed activity at 
its surface. Libby (1947) has shown on empirical grounds that an activity equal 
to one-tenth of the total number of disintegrations occurring in a layer equal 
in thickness to the beta-particle range will be observed at the surface of a thick 
source for the case of a normal energy spectrum. ‘The specific activity, and 
hence the minimum half-life ascribed to the element, therefore depends on the 
value assumed for the disintegration energy. ‘The values quoted in table 2 
correspond to a disintegration energy of 2:0 Mev. For any other assumed 
disintegration energy E the minimum half-life to be ascribed will be increased 
by a factor ¢ proportional to the f-particle range at energy E. Values of ¢ 
calculated from the range—energy relation given by Katz and Penfold (1952) are: 
E (mev) 0-25 0:5 1:0 ile 2 3 + 5 10 

0-06 0:17 0-44 0-62 1 LS 24 Dy 4:8 

From this it can be seen that few of the minimum half-lives ascribed are 
likely to be in error by more than a factor of three or four for reasonable values 
of the energy available. 

Considering now the detailed results given in table 2, it is clear that most 
of the elements examined show no evidence of any real activity. In several 
cases unmistakable f-activity was shown by the samples used, but the correlation 
between observed a- and f-activity is consistent with the assumption that 
practically all the activity is due to contamination by heavy radioactive elements. 
The only element for which clear evidence of activity exists is molybdenum and, 
even in this case, work on other samples with different chemical histories is 
desirable before the activity can be regarded as proved. 
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§5. CONCLUSION 


The simplicity of the method used has made possible the examination of a 
considerably larger number of elements than could have been done by previous 
methods. In special cases, such as those in which inert gases are produced in 
the reaction, far weaker activities could be detected by the ‘geological’ method 
used by Inghram and Reynolds (1949). Where the inert gases themselves are 
suspected of activity, very high pressures would be required if the present 
method were to be employed, and the sensitivity would not be comparable 
with that available by using a well-screened counter-controlled cloud chamber 
for example. For most substances, however, the method appears capable of 
competing in sensitivity with other existing techniques. 

Experimental values for the mass differences of the heavier pairs of isobars 
are greatly needed; it is believed, however, that the large number of elements 
which have been examined without yielding unmistakably positive activities 
has an appreciable cumulative value without these. At least a few pairs of 
isobars would be expected to possess sufficiently large mass differences to show 
marked activities if the no-neutrino hypothesis of double beta-decay were 
correct. Although, therefore, the minimum half-lives observable by the present 
method are not as large as those which have been recorded previously for one or 
two special cases, they represent useful evidence against this hypothesis in its 
simplest form. 
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ABSTRACT. Expressions for the angular distribution of tritons emitted from nuclei 
bombarded with deuterons have been obtained using Born’s approximation and the non- 
perturbation method of Butler. Discussion of the effect of different forms of the triton 
wave function on the angular distribution is given. The theory is shown to explain the 
experimental angular distributions satisfactorily. 


THE THEORY OF (d,t) REACTIONS 


XPERIMENTS on the emission of tritons from nuclei bombarded with 
Kk deuterons reveal angular distributions similar to those observed in (d, p) 
and (d,d’) reactions. This leads one to expect that this reaction also is 
a direct process, proceeding from the initial to the final state without the 
formation of a compound nucleus. The mechanism proposed here is that of an 
inverse stripping reaction. One imagines a neutron attached to the initial nucleus 
to be drawn out and captured by the deuteron in its flight without the deuteron 
actually penetrating the nucleus. 

There are two methods available for treating this type of reaction. One 
treatment is by Born’s approximation, which has been used with success in the 
interpretation of experiments on (d,p) reactions (Bhatia et al. 1952), and the 
inelastic scattering of deuterons (Huby and Newns 1951). ‘There is also a method 
which avoids the use of perturbation theory, and this has been used by Butler 
(1951) in his theory of the (d, p) reaction. Both these methods have been applied 
to the present problem.* Details of the calculation need not be given here since 
both methods have been published in full for the (d,p) reaction, and the 
modifications involved are straightforward. 

The differential cross section can be written 


o(6) = |G(p) PEAi| M4) P, 


where the scattering angle appears in the quantities k and p which are given by 
= ((Mg|M,)ha — hy)? + 4( My) Mi hgh, sin® 0/2, and p? = (hq/2 — hy/3)?-+ Shaky sin® 62, 
k, and k, being the wave numbers of the deuteron and triton given by 
Ra =(2m,Eq)'?/h and k, = (2m,E,)"?/h, where mg and m, are the reduced masses of 
the deuteron and triton and Fy and FE, the energies of the deuteron and triton 
in the centre-of-mass system of coordinates; M, and M, are the masses of the 
initial and final nuclei respectively. The values of J, which is the angular 
momentum with which the neutron leaves the nucleus, are restricted by the 
selection rules J;= J;+1+4, J; and J; being the spins of the initial and final 
nuclei respectively, and by the fact that / is even or odd according as the parities 
of the initial and final nuclei are the same or different. The A, contain all factors 


* I have not in fact used Butler’s technique, but have followed a more elegant presentation 
of this non-perturbation method due to Dr. R. Huby, which involves the use of S-matrix theory. 
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not affecting the angular distribution. ‘The factor M,(k) exhibits the maxima and 
minima which are characteristic of these reactions. This factor differs in the two 
methods of calculation and is, for Born’s approximation, 


aN 
Mih)= (see) Fesanla), seen (1a) 


and for the non-perturbation method, 


Mk) . yy hry) Kysgio(kto) — ( z yrs hr) K, 
=K1 | =— r, ig(K19) — RFo | r To)» 
I KT9 (oe 1+1/2( 0) 143/2\KTo 0 2kro 1+3/2( 0) 141/2(* 0) 


« being given by « = {2m,(«,—Q-—,)}/h, €q and «, being the binding energies of 
the deuteron and triton, m, the reduced mass of the neutron with respect to the 
final nucleus, and QO the Q-value of the reaction. The quantities a and 7, are 
both lengths of the order of magnitude of a nuclear radius. 

The factor G(p) is the same in both the perturbation and non-perturbation 
treatments of the problem, and depends on the internal wave functions of the 
deuteron and triton. One would expect this factor to depend sensitively on the 
form of the internal wave function for the triton and particularly on the asymptotic 
form, since one requires essentially the probability that in a triton one particle 
is a considerable distance away from the other two. However, this is found not 
to be the case as the calculations given below show. It is also found that the 
deuteron wave function has little effect, which is to be expected. The factor G(p) 
is calculated for three combinations of wave functions (we omit all multiplying 
constants). 

(a) Using Irving’s (1951) wave function for the triton, which has the correct 
asymptotic form, and an exponential wave function for the deuteron, 


xe =exp {— B(ZR; 7)? }/(2R,?)", Xa = exp {—4#Ryo}/Ryo, 
sin (4 —2n) ihe 
(A?+ B22" 

where A =1 + (4/2a/+/38)x, B= oe —x?)'l2,%4=tan+B/A. Incalculations 
we take « =0-23 x 104% cm7, and Irving’s values n=} and £8 =0-759 x 108 cm7}, 


which give the correct binding energy for the triton. 
(6) Using a gaussian wave function for the triton, and an exponential wave 


function for the deuteron, 
xe=exp {—Y2R,,"}, Xa =exp {—aRjo}/Ryp. 
(c) Using gaussian wave functions for both the deuteron and the triton, 


we obtain G(p) = 5], 


X=exp{—yZRi*}, — Xa=exp{—27Rie”}, 

we obtain CU exp cel aye ee EP ONE, (2 (d) 

In calculations we take y=0-25 x 10° cm since this gives the correct difference 
for the binding energies of 3H and *He when the coulomb energy is taken into 
account (Bruno 1948). It also gives a triton of about the same size as Irving’s 
wave function, as is seen if one calculates the mean distance and the mean square 
distance of one particle from the centre of mass of the triton. ‘The cases (4) and (c) 
are interesting since the angular dependence is the same, and it shows that the 
deuteron wave function has very little effect on the angular distribution. A plot 
of | G(p) |? is shown in fig. 1 for case (a) and cases (4) and (c), and it is seen that 
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the variation with p is very similar. Thus the form of the triton wave function 
also seems to have little effect. ‘This increases confidence in the results obtained. 
It also makes the calculations easier since the gaussian wave functions give an 
answer in terms of elementary functions. For higher values of p, however, the 
difference between the results (2a) and (2b) becomes more marked, and in this 
region one would expect Irving’s wave function to be the more satisfactory. 
These results can be applied to the angular distribution of tritons emitted 
from the °Be(d,t)®Be reaction. This has been measured by El-Bedewi (1951) 
and extended to lower scattering angles by Holt and Marsham.* 'The deuterons 
had energy 7:7 Mev in the laboratory system and the Q-valuc is 4-67 Mev. 
Figure 2 shows the experimental points together with curves drawn from the 
two theories (1a) and (14) using G(p) as given by (2a). The transition is /=1 


1-0 


+ Experimental Points 
(Holt and Marsham) 


ba 
> 


LJ 
ss 


o(6) (arbitrary units) 


02 


0 20 40 60 
6 
Fig. 1. Plot of G(p) against p. A: G(p) given Fig. 2. The °®Be(d, t)*®Be reaction. Theo- 
by (2a). B: G(p) given by (2b). G(p) retical angular distribution as given by 
has been normalized to unity at p=0. Born’s approximation (full line curve) 
and the non-perturbation method (broken 
curve). 


and the values of the radii used are, in (1a), a=6:7 x 10-¥ cm, and in (13), 
7)=5:0 x 10-43 cm. Using the selection rule, we find that if °Be has spin 3 and 
odd parity, °Be has spin 0, 1, 2, 3 and even parity. 

Holt and Marsham have also measured the angular distribution of tritons 
from the 7Li(d, t)®Li reaction, and this is also consistent with an /=1 transition, 
which means that if 7Li has spin 3 and odd parity, °Li has spin 0, 1, 2, 3 and even 
parity. ‘These results do not give any new information, but they are consistent 
with the usually accepted values of spin zero for *Be and spin 1 for ®Li. 

Comparison of the two methods is difficult since both seem to fit equally 
well. ‘The Born approximation result has the advantage of greater ease of 
numerical calculation, but the non-perturbation method is possibly the better 
founded theoretically, and has a radius which seems to be remarkably stable at 
5x10" cm. This latter method appears to fail in the (d, p) reaction, when the 
neutron enters the nucleus with energy so large that the final state is virtual. 


* T am indebted to Dr. J. R. Holt and Mr. T. N. Marsham for communicating their results 
before publication. 
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(A fuller discussion of this question of radii is given in Bhatia et al. 1952.) 
However, in this case this will not happen since the final nucleus in the (d, p) 
reaction corresponds to the initial nucleus in the present reaction, and one always 
starts with a stable nucleus. 

Finally, Wolfgang and Libby (1952) have measured the total cross section 
of the *Be(d,t)*Be reaction and have found it to be approximately constant at 
about 2 x 10-*° cm? over the range of energy from 3 to 8Mev. It is difficult to 
estimate this theoretically, but a rough estimate can be made of 10-?> cm? at 
8Mey. In the limit of high energies the cross section varies as 1/Hy. ‘These 
results are in as good agreement as can be expected with the experimental data. 

In conclusion it may be pointed out that the above analysis will apply also. 
to the (d,*He) reaction by changing all quantities relating to the triton to the 
corresponding ones for the *He nucleus. ‘The parameters in Irving’s wave 
function remain the same since these give the correct binding energy for the 
3He nucleus when the coulomb energy is taken into account. 
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Electron Capture 
I: Resonance Capture from Hydrogen Atoms by Fast Protons 
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MS. received 9th Fune 1952 


ABSTRACT. Born’s approximation is used to calculate the cross section associated with 
resonance charge exchange between protons and hydrogen atoms, it being pointed out that 
in earlier treatments of this problem unjustified simplifying assumptions were made. An 
energy range of up to 250 kev is covered. ‘The predicted energy variation of the cross 
section is in harmony with the results of Ribe’s recent laboratory work on collisions of 
protons with hydrogen molecules. 


$1) INTRODUCTION 


LECTRON capture through energetic collisions between positive ions and 
neutral atoms or molecules, besides being of interest in itself, is important 
in aurorae and in certain other natural phenomena. It is therefore 

desirable to carry out quantal calculations on the process. ‘The simplest 
possible case 
H++H(1s) > H(1s)+H*+ 
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was studied many years ago by Oppenheimer (1928) and by Brinkmann and 
Kramers (1930) who gave a closed formula for the cross section which 
could readily be generalized to more complicated collisions. Unfortunately, 
examination of their treatment shows that the interaction potential employed 
is incorrect and that as a consequence so also are both the absolute magnitude 
and the energy variation of the deduced cross sections. The present paper 1s 
devoted to a reinvestigation of the problem. Attention will be restricted to the 
collision represented by (1), but the method is applicable to other cases and 
will be extended to them later. 


§2. ANALYSIS 
2 
At high energies the cross section Q associated with the rearrangement 


collision ; 
(nucleus A+ electron e)+nucleus B+ nucleus A+(nucleus B+ electron e), 
is given by Born’s approximation (cf. Mott and Massey 1949) to be 


3 V2 +1 
Qa em pd(cos®) eee (2) 
iv -l1 


in which M is the reduced mass of the system, v; and v, are the magnitudes of 
the initial and final velocities of relative motion, 0 is the angle between unit 
vectors n, and ny parallel to these velocities, and M is the relevant transition 
integral. Denoting the position vector.of the centre of mass of (A+e) relative 
to B by pe, that of the centre of mass of (B+e) relative to A by o, and that of e 
relative to A and B by r and s respectively, 
| | =| SJ exp [e(Ain,. e — Rene. o)]b,(r)o_*(s)V(r, s)dsdo| ...... (3) 
27v; (M,+m)My k _ 2nv_e (Mg +m)M,y (4) 
h MO Motm 2 > i a, 
¢; is the initial and ¢, the final wave function of the electron, and V(r,s) is the 
interaction potential between A and (B+e) which is 
VA me ACA 
r  [r—s] 
Z,e and Zye being the nuclear charges. Previous workers have taken only the 
first term of V(r,s) into account, asserting that orthogonality would make the 
contribution from the second zero if proper wave functions were used. This 
simplification is without justification. In the standard collision theory 
treatment (cf. Mott and Massey 1949) the complete wave function ‘Y(o,s) is 
expanded in the form 


where k= 


UF ;,(6)¢,n(S) ‘ie ee (6) 
and substituted in the Schrodinger equation yielding 
2 2M - 
(V2 + ky?) F(a) = ae F(a, 5) V(r, s)oe*(s) ds 5° ee (7) 


Bm 
from which (2), with defined as in (3), follows at once on making the usual 
approximation that 
'Y(o,s) exp G2niep)odh) 2 an eee (8)t 


t+ In the case of identical nuclei, the wave function should, of course, be made antisymmetrical, 
but this is an unnecessary complication in an investigation concerned only with exchange scattering. 
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Clearly the contribution from the controversial term of V(r,s) would not vanish 
even if, in place of (8), the exact expansion (6) were used. It would instead 
give on the right-hand side of (7) the function 


822M Lpe” 

are { — eT } F(G') Fagen’ aye as oes (9) 
which is certainly not negligible, being in fact equivalent to a coulomb potential 
on the left-hand side! The corresponding function resulting from the adoption 
of (8) is very different since the ¢;(s) content of exp (zk,n,.e)¢(r) is far from 
equal to F;,(c). This is, to be sure, most grievous but it is not a reason for 
ignoring the —Z,Zpe?/|r—s| term of V(r,s) altogether. The effect of Z,e?/r, 
the remaining term of V(r,s), is also likely to be poorly represented. However, 
the similarity of the two terms suggests that the errors may partially cancel, 
especially if Z, and Zp are unity. In consequence approximation (8) may be 
not unsatisfactory. ‘The situation is analogous to that arising in the more 
familiar case of excitation by electrons for which the interaction potential is of 
essentially the same form as that under discussion. Some early investigators 
thought that owing to orthogonality one of the two terms should be omitted. 
It is now known that such a procedure is incorrect and leads to serious 
overestimates of the exchange cross sections, particularly at low. impact energies 
(cf. Bates, Fundaminsky, Leech and Massey 1950). 


peep: 


To evaluate -@ it is convenient (following Brinkmann and Kramers) to 
change the integration variables, and introduce the vectors 


a=kenetRhin,M,/(Myt+m) sevens (10) 
B= —Rn; — ken; My/(Mp +m) overs 4 one (11) 
so that (3) may be written, 


os (12) 


| M |= | | exp [i(a.r+B.s)|¢,(r)db¢*(s) a — a=} dr ds 


If each of the nuclear changes is Z, and if ¢, and ¢, are 1s wave functions 
(12) becomes 


| a|= AE* [exp f(a.r+B.s)—alr+s)]{- - ater ds 


r |r—s| 
where a@=Z/a), a) being the radius of the lowest Bohr orbit of hydrogen. 
Moreover, it may easily be shown that K, the modulus of both « and 8, and 4, 
the angle between these two vectors, are given with sufficient accuracy by the 
simple equations 

2p2 2 2 ; 
Ken (1 +008 8)+ 3, ieee 2 th (14) 
with p=27mva,/h, the subscript on v being omitted as unnecessary since the 
collision is elastic. 

Integrals similar to (13) have been discussed by Massey and Mohr (1931) in 
connection with another problem. The contribution from the first term of the 
interaction potential may be obtained by elementary methods; that from the 
second term may be obtained by choosing the z-axis parallel to either @ or 8, 
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expanding 1/]r—s| in the standard Legendre polynomial series, and using the 
formula 


e+ 
(—a)"} exp (au) P(e) du=(2m/2)"F ngapl®) =2jnl() ++ (15) 
(cf. Watson 1944). Writing t= K/a it is found that 
4nZe* . 
jel =|PZE Ly —2(—1)"P,(cosd)gu()}|p vee (16) 


where f(t) =8/(1 +27)8 


Slt) =40Z | e""7],( Kr) 
0 
x [om | Cae set Ks ase ee | ergy CRS) as| dr 
0 


=80°Z | ©, (Kr)i,(Kry" dr, 
0 


with i,(Kr) = | e-48gn425 (Ks) ds = — Koa 2 (K-02, (Kr). 
(0) 


Simple analysis shows that 
Zi 
= = DY 3 ag 
£,(o= Pa AW) B,(t) log (1+¢7)—tC, tan #], 
A,(t) and B,(t) being polynomials and C,, being a constant. ‘The first six 


members of the series (some of which have also been derived by Massey and Mohr) 
are as follows: 


A,(t)=t* +5 A,(t) =3t4+ 214-2 

Bit) =0 B,(t) =6t? -12—2t- 

Car Cy 832 

A,(t) =5t* + 852? + 24 + 12t- A,(t) = 7t4 + 2592? + 438 + 4982-2 + 1802-4 
B,(t) = 3002 +0+ 30¢-2 + 12¢-4 B,(t) = 84¢? + 252 + 672t-* + 588¢-4 + 1802-* 
Cy = 96 Corel 


A,(t) =9t* + 6152? + 2500 + 4(15710)t-? + 45202-4 + 14002-* 

B,(t) = 1802? + 1260 + 4800¢-? + 7380¢-4 + 522086 + 1400¢-§ 

Coe 320. 

A,(t) =11¢4+ 12322? + 9024 + 3024522 + 4584024 + 32550¢-§ + 88202-8 


B;(t) =33002 + 3960 + 2145022 + 50820t-4 + 613802-6 + 3696028 + 8820214 
C, =480. 


Finally, substitution of (16) in (2) yields 
2. sig n p* 2 2 
000 =[S) .MO-— DP, (sap —1) eal] ale) |naye. (17) 


‘The adoption of the infinite upper limit is justified by (14) which shows that when 
0 is zero K, and hence t, is extremely large. It may be noted that p? ,the variable in 
which Q is expressed, is m/M, times the energy of relative motion of the incident 
nucleus measured in rydbergs; thus if the incident nucleus is a proton a p? of 
unity corresponds to an energy of 24-9, kev. 
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If the g,, terms are omitted the integration in (17) can be carried out analytically 


giving 
oe") =| = : =(44Z3)-5] rag iee % elaarie cbs (18) 


which is, of course, identical with the formula obtained by Brinkmann and 
Kramers. 

It is instructive also to consider the effect of ignoring the momentum transfer 
occurring when the electron is removed from the stationary to the incident nucleus. 
If m is neglected in (10) and (11) it can be seen that in this case the vectors a and B 
are antiparallel. Hence P,,(cos8) is (— 1)", and the summation in the expression 
for Q can be replaced by 


1 4 8 
28 alt) = aire +e) 3(1 +2)? te esoai Seen, lelece. e8).0 (19) 

as may be proved by the consideration of | 
[| rpaeyexr Kn. (r—s)—a(r-+9)] de 3,078 Whee et (20) 


which can be evaluated easily (by, for example, using spheroidal coordinates). 
Since in addition the lower limit of the integral involved in it becomes zero, (17) 
thus reduces to O(p?) = {4/135p?}[1728 — 2592Z + 1067Z?|7a,? or, if Z is unity, to 


O(p?) =[6-01/p?]7a,°. Ae eA 


§3. RESULTS 

The functions f(t) and g,,(t) (n=0 to 5, t2?=0 to 100) were computed from the 
expressions given in the previous section. Great cancellation occurs when t 
is small but fortunately a useful check on the work can be obtained from (19) since 
for such ¢ the g,, series converges very rapidly. The figure shows the form of the 


Log (2 or gf) 


Log of functions f(t), go(t), 21(t), go(t), ga(t), ga(t) and g,(t) plotted against 2°. 


functions in the earlier part of the range. A logarithmic scale is used. It will be 
noted that the g,,’s decrease comparatively slowly, becoming much larger than f at 
quite moderate values of #?, and further that one of them, go, is comparable with 
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fevenat the origin. _ However, owing to the oscillatory character of the multiplying 
Legendre polynomials their influence on the magnitude of Q is not as great as 
might be supposed at first. 

Formula (17) was evaluated by numerical methods with Z unity and with p* 
between 0 and 10 (which covers energies up to some 250 kev in the case of protons). 
Account was taken of only those members of the g,, series that are illustrated. The 
omission of the higher members probably makes the calculated cross sections 
slightly smaller or larger than they should be according to whether p? is below or 
above about 2, but from some calculations on the effect of including g, and g7 
(which were estimated by extrapolation) it seems most unlikely that an error of 
as much as 10°% could be introduced. 

The cross sections finally obtained from formula (17) are displayed in the table 
together with the corresponding cross sections obtained from formulae (18) and 
(21). For convenience, the formulae numbers will be used as identifying sub- 
scripts. It can be seen that Qj, is considerably greater than Q,, at low energies 


Cross Section associated with Resonance Charge Exchange between Protons 
and Hydrogen Atoms 


seats Log of cross section (in units of 7a)”)T 

a formula (18) formula (21) formula (17)f 
(0-0) (1-7, —log p”) (0-7s—log p*) (0-73—log p*) 
0-4 1-95 Weil 0-9, 

0:8 1:4, 0-8, 0:55 

lie? 103 0-7 0:25 

1:6 0:7, 0:5, —0-0.5 

2:0 Or55 0-4, SU 

3-0 0:0. 0-35 6-7, 

4-0 —0-4 0-1, 4-1, 

5-0 —0°7; 0-03 —1-4, 

6:0 120, 0-09 —1-7; 

8-0 —1:5, —0-1, —=2:25 
10-0 —2:-0, —0:2, —2-6 


* The energy of the incident proton is (24-9,p?) kev. 
| 7a)" equals 8°8 x 10-47 cm?: 


} Formula (17) is the best approximation at present available. The results obtained 
from formulae (18) and (21) are presented purely for the sake of comparison. 


but falls off more rapidly as the energy is increased. ‘The omission of the g,,’s 
clearly provides a rather poor approximation. As would be expected Q,, equals 
Q,, when the energy is zero and falls off much more slowly so that at even moderate 
energies the difference between the two cross sections is very pronounced. The 
striking reduction of the probability of electron exchange due to the momentum 
transfer incurred is, of course, well known (cf. Mott and Massey 1949). 

‘The rehtability of Q,, itself remains to be considered. As the basic Born 
approximation is designed specifically for the high energy region it cannot be 
expected to be accurate in the low energy region. However, the range of its 
validity is not easy to predict. Results down to the zero energy limit are presented 
in the table in the hope that they may provide information on the matter when 
proper comparison data become available. At, and for some distance above, this 
limit, the cross sections given are certainly too large since in deriving them it was 
assumed that only transitions into the final state occur whereas backward transitions 


~ 
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must be important. A theoretical investigation of this is in progress. ‘The 
accuracy attained at moderate and high energies should be quite good though, as 
mentioned in § 2.1, assumption (8) is open to criticism in that what is relevant is not 
merely the general form, but a refined property, of the wave function (just as in the 
case of exchange excitation by electron impact, cf. Bates, Fundaminsky, Leech and 
Massey 1950). 

No experiments have been carried out on proton—hydrogen atom collisions, 
but several workers have studied H++H,+H’+H,+. Unlike (1) the process 
is not of the resonance type so the results obtained in the low energy region are 
of little significance in the present connection. The latest measurements in 
the moderate and high energy region are those of Ribe (1951). Over the range 
covered (for p? between about 1-5 and 6) the observed variation with energy is, 
within the stated experimental error, identical with that given by the theory. The 
absolute magnitudes of Ribe’s cross sections (per hydrogen atom) are, however, 
about 1-5 times those predicted. Part of the discrepancy may be due to the fact that 
the measurements relate to the total exchange cross sections and thus include the 
contribution from capture into excited states. Some calculations on this contri- 
bution are planned. ‘The cross sections may, of course, be influenced by the 
molecular binding but it would scarcely be anticipated that they are thereby 
increased. 
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ABSTRACT. ‘The transition curve for cosmic-ray bursts under lead at sea level is analysed 
to find the contributions from stars, extensive showers, electromagnetic interactions of 
p-mesons, and single high-energy electrons and photons. When this is done, a substantial 
fraction of the bursts under lead thicknesses less than about 2:5 cm remains unaccounted for. 
It is suggested that they are produced by ‘ narrow air showers’ consisting predominantly 
of photons of about 2-5 x 108 ev and having a sideways spread of less than 0-3 m. 


* Now at the Clarendon Laboratory, Oxford. 
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§1. INTRODUCTION 

FEATURE of the transition effect for cosmic-ray bursts in an ionization 
EN chamber appears to have escaped notice until recently (Northrop 1951). 
This is an abrupt increase in the burst rate under small thicknesses of 
lead, tentatively attributed to narrow showers, i.e. air showers having a small 
lateral extension. In the absence of shielding, bursts can be produced by nuclear 
disintegrations (stars) and by extensive showers (Bridge et al. 1948, Carmichael 
1948, Montgomery and Montgomery 1947, 1949). When the chamber is 
shielded, further burst-producing agencies become important. Bridge and 
Rossi (1949) account for the maximum of the transition curve at 10000 ft in 
terms of bursts produced by single high-energy electrons and photons. At 
larger thicknesses bursts can be produced through electromagnetic interactions 
of u-mesons (Christy and Kusaka 1941) or in nuclear interaction processes 
(Bridge and Rossi 1949, George and Trent 1951). It is the purpose of the 
present paper to point out that the above agencies are not sufficient to account 
for the transition effect of cosmic-ray bursts at small thicknesses of lead, and that 

some other agency must be presumed to be making a large contribution. 


§2. APPARATUS 

The apparatus consisted of complementary ion chamber and Geiger counter 
systems (Prescott and Mohr 1949, Campbell and Prescott 1952). The ion 
chamber was cylindrical, 57-6 cm long and 8-8 cm in diameter, filled to 
39 atmospheres with nitrogen. Three Geiger counters in coincidence were 
located immediately above the chamber; these acted as a ‘master group’. 
Counters located at 3 m and 5 m from the master group allowed bursts due to 
extensive showers to be identified. The efficiency of this arrangement for 
detecting extensive showers of the minimum density used in the analysis 
(600 particles/m?, 4:2 x 10° ion pairs) was 99%. The shielding consisted of 
hemicylindrical lead caps fitting over the chamber and extending one inch beyond 
the ends of the chamber and one inch below the centre on either side. Shielding 
thicknesses were 0-6, 1:2, 1:75, 2:3, 3-4 and 4:5 cm. The observations were 
carried out at sea level. 

§3. RESULTS 

Burst-size distributions were obtained at all lead thicknesses and, from them, 
transition curves; these were essentially the same as similar results obtained by 
other authors (e.g. Young 1937, Schmid 1941, Carmichael 1948, Chou 1948). 
The burst-size distributions for zero, 0-6 and 1-2 cm of lead show the discontinuity 
of slope characteristic of distributions made up of more than one component 
(Carmichael 1948). 

(a) A typical transition curve (for all bursts greater than 4-2 x 105 ion pairs) 
is shown in the figure, curve Z. ‘The statistical errors are within 1% for all 
points. Curve Y shows the transition effect for extensive showers as determined 
by the coincidence method described above. 

(b) ‘The rates of disintegration bursts arising in the walls of a shielded chamber 
can be found from the transition effect for stars in nuclear emulsions (Bernardini, 
Cortini and Manfredini 1949, 1950, George and Jason 1949, Schopper, Hécker 
and Kuhn 1951) and the observed rate in an unshielded chamber. These are 
equal to within a few per cent at all thicknesses of lead, since the transition effect 
is confined almost entirely to stars of three or four prongs which are relatively 
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inefficient in producing bursts in the present chamber. ‘The transition effect 
for wall disintegrations is shown as curve X. 

(c) Fast protons from disintegrations (Harding 1951) in the lead shielding 
itself might produce bursts at larger thicknesses. The transition effect for this 
supplementary contribution was calculated and is shown in curve W.* 

The main difficulty in constructing curve W is in determining the absolute 
rate. A straightforward calculation from the disintegration rates in lead gives 
only an order of magnitude (~0-1 per hour). An upper limit, which is all we 
shall require in the present case, can be found in the following way: Fahy (1951) 
finds that 68%, of the bursts under 10-5 cm of lead at sea level are due to u-mesons 
(knock-on and bremsstrahlung). Thus no more than 32% can be due to fast 
protons. Lapp (1946) gives the absolute rates of bursts of different burst sizes 
under thick lead shields. A burst of 4:2 x 10° ion pairs in the present chamber 


(per hour) 


Burst Rate 


Shielding (cm lead) 


Analysis of the transition curve in lead for bursts greater than 4-2 x 10° ion pairs 


is equivalent to a burst of 40 particles at minimum ionization, which gives, from 
Lapp’s data, a rate of four per hour for bursts produced by u-mesons and two 
per hour as an upper limit for those produced by fast protons. Curve W is 
extrapolated to two per hour at 10 cm. The true rate will actually be much less, 
since showers are also produced by nuclear collisions (George and Trent 1951, 
Fahy 1951). Indeed, Simpson, Baldwin and Uretz (1951) show that any 
increased contribution to the burst rate from stars in the lead is negligible under 
2:5 cm of lead. 

The transition curve for bursts produced by p-mesons was very kindly 
calculated for the author by Mr. I. D. Campbell and is shown in curve yp, fitted 
as to absolute rate in the above manner. ‘The important point about the curves 
W and vp is that their contribution is small for lead thicknesses less than about 
25 em. 


* This calculation is essentially the same as one outlined by Carmichael (1948). ‘The energy 
spectrum of the fast protons which enters into the calculation was constructed from Harding’s 


results (1951), 
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Fast protons from the atmosphere cannot produce bursts, since, under the 
prevailing conditions of gas pressure and bias, at least four such protons are 
required, even under the most favourable circumstances. 

(d) We now consider the contribution of single high-energy electrons or 
photons. The minimum burst-size of 4:2 x 10° ion pairs corresponds to 40 
particles at minimum ionization. The transition effect is estimated by finding 
the relative numbers of electrons that are capable of producing 40 particles or 
more under each lead thickness (Bridge and Rossi 1949, Bhabha and Chakrabarty 
1948). If we assume the integral energy spectrum of the initiating particles to 
be of the form E-* (Rossi 1949), where s is taken to be 2, we obtain curve V. ‘This 
curve includes the effect of fluctuations (assumed to be Poissonian), and the 
absolute rate is adjusted to account for as many bursts as possible. There appear 
to be no experimental observations at sea level with which this rate can be 
compared. The minimum energy involved is about 3 x10%ev. It is seen that 
single high-energy electrons or photons conspicuously fail to account for the 
bursts observed under small thicknesses of lead. A similar result follows for 
any other reasonable choice of energy spectrum. 

(e) Curve U represents the transition effect of the bursts that remain 
unaccounted for after allowing for all of the above agencies. A similar analysis 
roughly carried out on the data of Schmid (1941) and Chou (1948) leads to an 
identical conclusion, viz. that there are too many bursts at small lead thicknesses 
to be accounted for without introducing burst production by some additional 
agent. Similar general conclusions were presented in an earlier preliminary note 
(Prescott and Mohr 1949). 

§4. DISCUSSION 

It is of interest to consider what type of event could be invoked to account for 
the above result. The rapid multiplication in small amounts of lead requires 
electrons or photons as the primary agency. The maximum at about 1-2 cm of 
lead indicates an energy of about 2:5 x 10®ev and a multiplication at maximum 
of 4:5. ‘This means a minimum of, say, ten particles incident on the top of the 
shielding. If the particles are electrons this would mean an efficiency (from the 
geometry of the system) of better than 25% for the master group, i.e. a 
chamber—master group coincidence rate of about five per hour. The rate 
observed is only 0:58 + 0-07 per hour, of which roughly half can be accounted for 
by extensive showers (Campbell and Prescott 1952).* This suggests that at 
least the greater part of the particles are photons which cannot discharge the 
counters. Support for this suggestion comes from evidence for the lateral 
spread of the showers. 

The coincidence rate between two ionization chambers was observed at a 
separation of 0-32 m under shields of 0-6 and 1-2 cm of lead. No coincidences 
were recorded that could not be accounted for by extensive showers (Campbell 
and Prescott 1952), indicating that the sideways spread of the burst-producing 
showers is less than 0:32_m. Coincidences were also observed between the 
(unshielded) master group and counters at separations from 0-1 to 5m. The 
coincidence rate did not increase significantly until the separation of the counters 
was less than 0-3 m, but then increased abruptly; this type of result with Geiger 


* Of the balance, a burst can in isolated instances be identified as due to a narrow shower, 
e.g. very large bursts which fail to discharge another chamber only 0:32 m away, but the statistics are 


too poor to attempt to draw any detailed conclusions about, for instance, the proportion of ionizing 
particles present in the shower. 
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counters at lower densities is already well known (Wei and Montgomery 1949, 
Wei 1950, Alichanian and Asatiani 1945). The small spread for the shower 
particles shows either that they are produced very close above the apparatus or that, 
once produced, they are not scattered appreciably sideways. This would be the 
case if they consist mainly of photons as suggested above. The picture of the 
agency responsible for the large burst-rates under small thicknesses of lead 
consists, then, of a ‘narrow shower’ of photons of fairly low energy presumably 
accompanied by its parent particle and some slow electrons as well. If these 
“narrow showers’ are the same as those reported elsewhere (Wei and Montgomery 
1949, Wei 1950) they also contain penetrating particles. 

It does not appear possible without further experimental data to make any 
estimate of the absolute rate of narrow air showers without having to make 
arbitrary assumptions about the energy spectrum and number distribution of 
the particles in the shower and of their mode and point of origin. 
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On a Modified Definition of Riesz Potential for the Meson Case 


By L. S. KOTHARI 
Cavendish Laboratory, Cambridge 


Communicated by J. Hamilton; MS. received 3rd June 1952 


ABSTRACT. Amodified definition of the Riesz potential for the meson field is introduced. 
It is shown that the new definition is a generalization, in the w-plane, of the meson potential 
in the interaction representation. The relations satisfied by the meson potential in the 
interaction and in the Schrédinger representations are compared. 


§1. INTRODUCTION 


N a recent paper Auluck and Kothari (1951) have introduced a modified 
] definition for the Riesz potential for the electromagnetic field, and Kothari and 
Bhatnagar (1952) have shown that the new definition is a generalization of the 
Wentzel potential* in the «-plane, « being an arbitrary parameter characteristic of 
the Riesz potential. The purpose of the present note is to extend the modified 
definition for the Riesz potential to the case of neutral vector mesons. The 
extension seems to be of interest, as the new potential is a generalization of the 
meson potential in the interaction representation. ‘The definition of Fremberg 
(1946) for the meson case, on the other hand, reduces to meson potential in the 
Schrodinger representation when analytically continued to«=0. The equations 
satisfied by the meson potentials in the two representations are compared. 


$2 


We define the metric tensor g,,, a8 299 =13 211 =822 =833 = —1; 2,,=0 (uA). 
The velocity of light is taken as unity. The scalar product of two four-vectors 
A,, B,,is denoted by [AB] : 

[AB] = A,B" =A)B)— A,B,—A,B,— A,B; =A)By—(AB), 
where (AB) is the three-dimensional scalar product of the space parts of A,, and B,,. 
The (positive) length of the space part of 4,,is denoted by | A |. 


For the neutral vector meson case we define the Riesz potential 4,%(x) at 
any point x as 


A,o(3) = H(a) { (2? x?) A,, 14 Aye Nae, (1) 


where H(«) =2/T'(4«)I'(1 — 4a) =(2/7) sin 47a, k*=[RR], and x is the rest mass of 
the meson. ‘Fhe domain of integration D is the region defined by k? > y?, k, >0, 
k being the momentum four-vector of the meson. The above integral is 
convergent for«>0. ‘The coefficients A,,, are so far arbitrary. It is easily verified 
that 4,%(x) satisfies the equation 


(-4A7a=4, 7) (2) 


* It may be noted here that the Wentzel potential (Dirac 1947) and the potential for the 
electromagnetic field in the interaction representation are the same. 
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Taking f,(3) to represent the most general interaction between the meson 
and the nucleon field, we write 


1 
Anu = Tay | ful) CRS, Tenge ate ie eae (3) 


where the integral is taken throughout the space S given by 2) < x). Substituting 
this in (1), we have 


A(x) = ay eae ‘ (Ro — y2)\@-D2f (2) sin [k, x—2] d*kd*z. ...... (4). 


It is readily shown that A,**?(x), as defined by (4), vanishes when analytically 
continued to «=0. Combining this result with (2) and putting 4,°(«)=A,(x) 
we have 

Wale Ail) 0 ne er Be ont ve (5) 


Pichon 622 P = sie [ (=x 2)0-22F (2)k" cos [k, x2] dtkd'z. ....(6) 


OX, 


The right-hand side, in general, does not vanish when analytically continued 
to «=0. Thus eqns. (5) and (6) show that the Riesz potential defined by (4) 
(or by (1) and (3)) is a generalization in the «-plane of the ‘meson potential in 
the interaction representation’. In the limit «+0 the new definition reduces to 
ordinary meson potential in interaction representation. 

It is difficult to evaluate the right-hand side of (6) as « tends to zero in the 
general case. Physically also we are more interested in the interaction of a point 
particle with the meson field. We therefore proceed to evaluate (6) for this 
particular case. 

For a point particle of mesic charge g moving along a world line z=2a(r), 
7 being its proper time, the coefficients A;,,, become 


I J 1 .—ilkz’ / 
w= ay] _% Cena ee ohare (7) 
It follows from this that 
RAy,= foe kA, = fe, Te ee (8) 


Hence we obtain from (1) 


2A) OT (ue — yao sin [&, 2-2] a, 


Cx, 
SS ia oe rh ee Aine Co (9) 
where G(x) = — Rea emer sin (Ree 9 Scie (10) 


Now consider the integral 

P(x) = |, — x2)? elt dh, 
Changing the variables of integration from kp, ky, ky, kz to k, ky, ky, kz, where 
k? =k,2—|k|? (the Jacobian for the transformation is k(k? +|k|?)"), we have 


fS ae rexp {ix9(k? + | k|?)12 —2(Rx)} 
eerey 2 2__ 2)\0-2)2 b dk dT Ed pale 
ae ik bie a) (P+ [Re 
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To integrate over d3k we introduce polar coordinates, with the polar axis parallel 
to the vector (k,, ky, k;). Integration over the sa can now easily be carried 
out, giving 


T(x) = a ie (2 — x2) "22 dh : ‘ aaa eizol? + IE" feilel lel _ g—élal tel) 
2 2\1/2 
# ail, (k2 — y2)-22 adh sro | exp a x |} 
ere (11) 
The integral over a can be evaluated in terms of Bessel functions: 
inH,(ks) for x9 >| x| 
ie SE ee da= 4  inH,(ks) for |x|>x*)>—|x| 
Gis Me —inH,(ks) for —|x|>xp 
> Ae (12) 
where k =(k?)1? and s=(x9?— |x |?)¥?. 

Substituting back from (12) and (11) in (10) we get 

B(x) = H(«) [ (= gy a Ul ks\:drsiat 5 axe (13) 


Jy(ks) for xy >|x| 
where — U(ks) = 0 for |a [xp — Fen ee (14) 
—J (ks) for —|x|>%. 
The expression for G,*(x) is clearly singular on the light cone x)= |x|. The 
nature of the singularity on the light cone is of the ae 


ale oo i 7 Bl D- 8(% +] x))}= fc rap 
Thus we have 


D,x(x) = (2) | (AB — y2)-22R ‘a : al U+a()} dk. 


‘To continue the above expression to « =0 we introduce a new variable ¢, defined 
as t?=k?—y?. In terms of this 


(x) = H(«) be {7-1 im | a U{s(t? + y?)¥2} + Aa} at, Bee: (15) 
and hence 


D(x) = A(x) + Is [x] aI | U(xs) 


=A(x)= = = UGs\s/ af. ee (16) 


It will be noticed that the function Y,°(x) is the same as the D-function for a 


particle of mass x, used by Dirac (1934), Schwinger (1948) (Schwinger’s A) and 
others. 


The meson potential in the interaction representation thus satisfies the 
equations ; 
(-x)A,(x)=0, 0A, (0) /On) = — 220 a) 8” eee (17) 


These may be compared with the equations satisfied by meson potential in 
Schrodinger representation: 


(C1-x°).0,(x) =g0,, Ox0,(x)Ox,=0. oe, (18) 
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In the absence of nucleons the two sets of equations are the same. When 
nucleons are present it is the second equation that gets modified for the potential 
in interaction representation, while for the potential in the Schrédinger 
representation it is the first equation that is altered. It may also be noted that 
by putting x =0 in (17) we pass over to the Wentzel potential (Dirac 1947). 
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Alkyl Shifts in Absorption Spectra of Azulene 
and other Aromatic Molecules 


By  C-=47 COULSON* 
Wheatstone Physics Department, King’s College, London 


MS. received 18th March 1952 


ABSTRACT. An inquiry is made into the factors which cause shifts in the ultra-violet 
spectra of aromatic molecules on substitution. ‘The method of molecular orbitals is used, 
and the discussion is confined to z-electrons. It is shown that there are two chief effects : 
the first is essentially electrostatic, or inductive ; the second is essentially conjugative, or 
mesomeric. In alternant hydrocarbons the second predominates and causes a red shift. 
In other molecules (e.g. azulene, quinoline) the first effect is considerably larger and now 
predominates over the second. Different positions of substitution will lead to shifts, 
sometimes to the red and sometimes to the blue. Electron donors will cause shifts in opposite 
directions to electron acceptors. ‘The importance of these variations in characterizing the 
electronic transitions is stressed. Detailed investigation of the various methyl azulenes 
confirms the assignment of the long wavelength band with almost complete certainty. 


§1. INTRODUCTION 
T is a well-known fact that in normal aromatic molecules alkyl replacement 
] of one of the boundary hydrogen atoms causes a shift of the long-wave 
ultra-violet spectrum to the red. In azulene, however (Plattner 1941 and 
later papers), certain positions of the methyl substitution shift the band to the 
red, other positions shift it to the blue. And in heteronuclear molecules, as 
Brown (1952) has recently pointed out, the shifts vary considerably with position 
of substitution, and are often larger than the shifts with normal aromatics. 
(A ‘normal’ hydrocarbon is here defined as one with no odd-numbered rings: 
more precisely it is an alternant hydrocarbon, in the language of Coulson and 
* Now at Mathematical Institute, Oxford. 
PROC. PHYS. SOC. LXV, II—A 61 
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Longuet-Higgins (1947a).) Ever since Mulliken’s (1939) first account this 
alkyl shift has been interpreted as due to hyperconjugation. Yet, despite this, 
the phenomenon itself has by no means been satisfactorily explained; nor does 
its significance appear to have been fully recognized as a diagnostic tool whereby 
we can identify and classify electronic transitions in conjugated molecules. 
It is the object of this paper to give a fuller account than hitherto of the effect 
of alkyl and other substitution on the characteristic ultra-violet absorption of 
the parent compound. The questions to which we shall attempt to provide an 
answer are the following: 
(i) Why are normal alkyl shifts to the red? 

(ii) What factors determine the variation in shift of frequency with position 
of substitution ? 

(iii) What factors determine the absolute magnitude of the shift, and why 
are these shifts usually larger for non-alternant hydrocarbons and heteronuclear 
aromatics than for normal hydrocarbons? 

(iv) Why are these shifts for molecules other than normal hydrocarbons 
sometimes to the red, and sometimes to the blue, and why do different sub- 
stituents (e.g. —CH, and —COOH) often lead to shifts in opposite directions ? 

(v) How may an experimental determination of the directions and relative 
magnitudes of the shifts with different positions of substitution be used ‘to 
characterize the type of electronic transition involved ? 

We shall see that reasonably convincing answers can now be provided to all 
these questions. Essentially our argument rests on a distinction between two 
aspects of the total effect of substitution. We have called these the electrostatic 
effect and the delocalization effect. By this we do not mean that they can be 
entirely separated from each other, for they cannot. But the ideas behind our 
rough separation correspond to pictures which have been visualized by chemists 
for some time, and their combination certainly provides a model which accounts 
for all the main features (i)—(v) above. ) 


§2. THE MODEL 

Since the transitions with which we are concerned take place in the region 
associated with z-electrons, we deal directly only with these. Any o-electron 
effects which may exist are translated into terms applicable to z-electrons. 
As is usual in dealing with ultra-violet spectra, we find it more convenient to 
use the molecular-orbital (m.o.) description of the z-electron orbitals. In this 
each 7-electron is described by a wave function which is a linear combination 
of 2p7 atomic orbitals (a.o.) on the separate atoms. As far as possible we use 
the notation of Coulson and Longuet-Higgins (1947a). Thus, any m.o. in the 
unsubstituted molecule is of the form 


wre Cihy bins iC, Py Sule (SAY we eee (1) 


The energy € associated with this orbital is defined in terms of the coefficients c,, 
together with the Coulomb term «, of each atom and the resonance integral ,, 
of each bond. «, measures the electronegativity of the atom (i.e. the attraction 
of the atom for z-electrons), and f,, measures the strength of a double bond 
between atoms 7 and s (i.e. the attraction of the bond for z-electrons). We 
suppose that each «, and B,, is known. If necessary, as shown by Wheland and 
Mann (1949), by Altmann (1950) and by Chirgwin (unpublished), we may have 
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to determine the «, and £,, by some sort of iterative procedure. In what follows. 
we regard them all as known. ‘Then the coefficients c, for each separate m.o. 
are found from the normal secular equations. We suppose that these also are 
available. ‘They can, of course, be calculated separately for the substituted 
molecules also, as was done by Pullman, Mayot and Berthier (1950), but this, 
though it provides a-numerical value for the displacement of the band, does. 
not enable us to understand how these displacements arise, or to answer the: 
questions (i)—(v) of the introduction. 

In the simple form of the m.o. theory which we are using we calculate the: 
allowed energy levels and then allot the electrons to these levels, generally filling 
the lowest levels first. Each allocation of electrons gives rise to what we may 
call a configuration. We make the assumption that each such configuration 
by itself is a good approximation to a real state of the molecule. This means. 
that we neglect all configuration interaction of the kind discussed by Coulson, 
Craig and Jacobs (1951). It seems most unlikely that such configuration interaction 
would materially affect the small energy shifts on substitution with which we are 
here concerned. In a similar way we neglect all o-7 resonance of the type 
recently discussed by Altmann (1952). It would not be very easy to include 
this in a m.o. treatment, except through the inclusion of considerable configuration. 
interaction, and its inclusion here does not appear to be necessary in order to 
understand our phenomenon and calculate approximate values for the various 
shifts. In any complete treatment, however, it would be necessary to incorporate it. 

If we accept the foregoing assumptions, then the frequency of any band 
(strictly the 0 —O frequency, since we do not consider the vibrational pattern of 
the bands) is given very simply by the difference in energy of the two molecular 
orbitals, corresponding to an electron jump from y, ‘to yf: 


ee NE cot Oe wc ros (2) 
In what follows we shall need the formula for the energy of the m.o. . It is 
ce CV MAF Only Hoe FL CBia +++ 3) 


1+2¢,¢,Sj.+.... 
This formula follows immediately from the energy relation 


c= [utttar | | year, 


where H is the effective Hamiltonian for a single z-electron. We use the standard 
definitions of «, and f,,: 


p= | beHbedr, — Bre= | beHb.dr, 


and suppose that the overlap integrals S,, are defined by S,,=J¢,¢,d7=S,,, 
S,,=1. When we wish to distinguish the coefficients c, for the two orbitals %, 
and 7%, we shall write them c,, and c,. respectively. ‘The normalization in (3) 
is such that ¢,7+¢)2+....=1. 


sree W OLE PREC TES OF SUBSTITUTION 
Consider a methyl substituent. We know from the reactivity of toluene 
and from dipole moments, etc. (cf. Cohn, Hughes, Jones and Peeling 1952) 
that when it is substituted in an aromatic molecule the methyl group acts as 
an electron donor. We also know that it is able to conjugate with the aromatic 
nucleus of the original molecule. ‘These two effects together comprise what we 
61-2 
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call hyperconjugation. Since, however, this is a single phenomenon, we have 
no genuine right to distinguish the two partial components just listed. But, as 
stated earlier, it is interesting to do so, and helps us to understand the phenomena 
associated with substitution (cf. Matsen 1950). 

The electron donor effect is sometimes represented, in resonance theory, by 
structures such as (I), and the conjugation, or electromeric, effect by structures 
such as (II). In m.o. language (I) is very similar to an electrostatic influence 


a ee 
a. erat Zo See 
= NH at ~~ 
(1) (11) 


exerted on the aromatic nucleus, and could be called either the ‘electrostatic 
effect’ or the ‘inductive effect’, since it is induced in the original unsubstituted 
molecule by the presence of the substituent. The effect represented by (II) is 
much less dependent on this electronegative character: it is essentially an 
additional delocalization of the z-electrons of the unsubstituted molecule. It 
would exist even if the substituent were a vinyl group, where the electrostatic 
effect would be much reduced. For that reason we shall refer to these two aspects 
of substitution as the electrostatic effect and the delocalization effect. As we shall 
see, their influences on the spectral properties of the molecule are quite distinct. 
And, by means of the distinction, we are able to provide intelligible and semi- 
quantitative answers to the original questions (i)-(v). 


§4.. THE ELECTROSTATIC EFFECT OG ae 


We suppose that the substituting group X is attached at position 7 of the 
conjugated molecule. ‘The electrostatic effect arises from the fact that the bond 
C,—X is polar, so that electrons are either pushed towards C, (electron donor) 
or away from C, (electron acceptor). In either case we can say that there will 
be a small change in the effective electronegativity of atom r. Since, at this 
stage, it is chiefly o-electrons that are being displaced, this effect is essentially 
the inductive effect of Cohn et al. (1952). Now the electronegativity of atom r 
is measured by the Coulomb term «,. So with electron donors such as —CH, 
we shall have 6«,>0, and with electron acceptors such as —COOH we shall 
have 6a,<0. ‘This change in «, will affect all the z-electron m.o.’s and alter 
their energies. According to (3) the alteration in energy will be 


dc 
oe = (=) 6a,.=(c 2 Noe) oe ee (4) 
where N is written for 1 +2Xc,c)S,., and we have supposed that Sa, is sufficiently 


small for first-order perturbation theory to be applicable. We shall return to 
this condition later. 


It follows-from (4) that 


Con" Cy : 
ea) = 48 ~ Fe} Be, vine EO 


Now for the top occupied and lowest unoccupied orbitals in most conjugated 
systems 22c,cySj, is fairly small, of the order of 3. It is therefore a reasonable 
approximation to put N,=N,=1. Equation (5) then shows that 


8(Av) =(C,9? — ¢,42) x 8a, 


ao 
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This result allows us to make two important deductions: 

(a) Since d(hy) is proportional .to 5x,, we should ‘expect that the effect of 
substitution will change from bathochromic to hypsochromic, or vice versa, if 
we change the sign of 5«,. This gives us an immediate explanation of point (iv) 
in the introduction. For, as Plattner and colleagues (1948, 1951) have shown 
experimentally, carboxylic substitution affects the long wavelength band of 
azulene in the opposite direction to methyl substitution. This is exactly what 
would have been expected if the displacement represented by (6) dominated 
the phenomenon. 

(6) With the same substituent, but with different positions of substitution, 
we shall expect 6(Av) to vary as (¢,.?—c¢,,)._ Now a simple consideration of these 
coefficients in the ground and excited states shows that this expression will be 
positive for certain 7, and negative for others. ‘Thus a single substituent should 
sometimes cause hypsochromic and sometimes bathochromic shifts, in complete 
accordance with the experimental results summarized by Brown (1952). 

This completes our answer to point (iv), but it does not deal with one very 
important situation. This is the case of alternant hydrocarbons. For here, as 
Coulson and Rushbrooke (1940) have shown (see also Coulson 1948), the longest 
wavelength allowed transition is associated with molecular orbitals 4, and #, 


where, in (1), we have Cherie Vals on Seek (7) 


PR (6 ge 2 OD, ES 2 Aree (8) 


and the summation is taken over all neighbours 7s. If we call this sum A, then 
A>0, and is usually of the order of magnitude of 4. If we insert (7) in (6) we 
get a zero result. This means that we are no longer at liberty to neglect the 
overlap integral S in (5). Using (7) and (8) in (5), and supposing that S has the 
same numerical value (approximately +) for all neighbouring atoms, and that 
S'=0 for all non-neighbours, we may write 


2AS 
5(hv) = 1. 2282 Cer Od 
This is effectively 
Di Ry)=2AseexKOee nln ee (9) 


We are to use (9) rather than (6) when dealing with alternant hydrocarbons. 
But this immediately answers some of the other questions (i)-(v). For 

(c) The shifts are smaller than before. ‘This is because 2AS is of order 4; 
we should therefore anticipate shifts of the order of a quarter (or less) of the previous 
magnitudes. In this way we deal with the second part of question (ili). 

(d) The shifts associated with the same substituent, but at different positions 7, 
are now all of the same sign. For methyl substitution, where 6«,>0, we should ~ 
expect 6(hv)>0, leading to a hypsochromic shift. ‘This is the opposite of what 
is actually found, and it shows us that the delocalization effect must be considered 
as well as the electrostatic effect. Since we have seen in (c) that these latter are 
smaller in alternant hydrocarbons than in other aromatics, it is not surprising 
that delocalization effects should be relatively more important in the former 
than the latter. 

(e) In comparing a heteronuclear molecule such as quinoline with its corres- 
ponding alternant hydrocarbon, here naphthalene, we should expect corresponding 
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bands in the heteronuclear molecule to be shifted to the red as compared with 
the hydrocarbon. This is because Sz, is negative for a nitrogen atom relative to 
a-carbon. Shifts of this kind are frequently observed. Our analysis, however, 
deals only with 7—>7* transitions, and says nothing about possible transitions 
of a non-bonding nitrogen o-electron to an anti-bonding z-orbital. Recent 
analysis suggests that these may play an important part in the complete spectrum 
of the heteromolecule. They do not exist in the hydrocarbon. 


§5. METHYL AZULENES 

Before concluding this discussion of the electrostatic (or inductive) effect 
we should like to show how powerful a tool it is in the characterization of these 
electronic transitions. This is a very complex problem, and so far very few 
really definite conclusions have been drawn for polynuclear molecules other than 
benzene. It will be convenient to illustrate our argument by specific reference 
to the methyl azulenes, for which excellent data are available. Pullman, Mayot 
and Berthier (1950) have already shown that the directions of the methyl 
displacements in the five distinct positions of substitution are in accord with 


S A 
€19 
£9 
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theory. But these authors do not appear to have realized how remarkably 
clear-cut the argument is, and how definitely it characterizes the transition 
itself. 

Let us assume, then, that in the methyl azulenes the electrostatic effect is 
more important than the conjugation, or delocalization, effect. In view of our 
earlier comments and the work of Cohn et al. (1952) this seems quite probable. 
What we can now do is to discuss the detailed changes that are induced in the 
energies of the molecular orbitals as a result of the change da,. F igure 1 shows 
the numbering system for the molecule. We have also shown, in two columns, 
the energies of the orbitals for the unsubstituted compound. These are numbered 
in ascending order and fall into two groups, here called S and A for short, 
according as the m.o’.s themselves are symmetrical or antisymmetrical with 
respect to reflection in a plane perpendicular to the molecule and passing through 
its axis of symmetry 2-6. The actual coefficients c, have been calculated by 
several authors (e.g. Coulson and Longuet-Higgins 1947 b, Wheland and Mann 
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1949), using simple assumptions about the constancy of resonance integrals 
and Coulomb terms. In the present work, however, new values due to 
B. H. Chirgwin have been used. These include allowance for uneven 
distribution of charge on the various nuclei, and for differences in f,, and S,, 
due to unequal bond leagths. The writer is grateful to Mr. Chirgwin for 
providing these in advance of publication. 

If we refer to these m.o.’s in terms of their energies, there would appear 
to be five states of the molecule, comprising the ground state G and four excited 
states A-D, which might be expected to be significant in the visible and ultra- 
violet region. ‘These are shown below, together with (1) the type of transition 
assuming that absorption takes place from the ground state, (2) the direction 
of polarization, which must be either parallel or perpendicular to the symmetry 
axis, and (3) the nomenclature proposed by Mann, Platt and Klevens (1949, 
referred to as MPK) in their comparison of the observed spectra of azulene and 
naphthalene. These transitions are listed in increasing sequence of energy, 


Occupied orbitals _ Transition Polarization MPK designation 
G- €,” €,” «7 €," €;7 
A- e,7 €,7 €57 «7 €s € €5—>€, ft Ly 
Ry Gee Ge GA GaiGr €4—> €g I (bs, 
Cres ees Ne,” C, Ee, €5—>€7 TT By 


1D) Rar ear oe €3—> € A K) 


and were identified by MPK with four transitions actually occurring in the 
absorption spectrum in this order. ‘This identification is plausible and agrees 
reasonably well with the energy differences calculated. But, as a growing body 
of evidence for naphthalene is beginning to show, great caution is needed in 
interpreting transitions in conjugated molecules in any simple manner. What 
we are now going to show is that the effects of methyl substitution allow an almost 
unequivocal decision to be reached concerning the long wave transition. 
Table 1 shows the values of c,.2—c,, for the four possible transitions. An 
italic figure in this table corresponds to a hypsochromic shift on methylation 
at that position. 
Table 1. Values of c,.2—c,,? 


Band Atomr= 1 2 4 5 6 
A —0:28 0:10 0:19 —0-11 0:25 
B —0-05 —(0-21 0:16 —0-01 0-77 
Cc —0:18 0 0:09 0-17 0 
D —0:05 0:10 —0:02 —0-11 0:25 


Now the positions of hypsochromic shift with the long wave band of azulene 
are precisely 2, 4 and 6. This, as Pullman and colleagues showed, agreed 
excellently with the interpretation of this band as being A. But now we see that 
there is no other possibility. It is only if this band is really due to a transition 
of type ¢;>«, among the allowed molecular orbitals that even the correct signs 
are obtained for the shifts. It is true that we have here neglected the conjugative 
effect. But we shall show shortly that this effect can only shift energy bands 
to the red: no hypsochromic shifts are possible by this means. Hence only 
band A can be the visible band. In addition to this, we have evidence that the 
conjugative effect is of smaller magnitude than the inductive effect, since it is 
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not able to change the figures in the table by as much as 0-1. This discussion 
deals partly with point (v) at the beginning. But the agreement is really quite 
surprisingly good from another aspect—the relative methyl shifts at different 
positions. Absolute values cannot be estimated since we do not know the magnitude 
of Sa,; but if we choose Sx, so that the shift at position 1 agrees with experiment, 
then we get table 2. At the bottom of the table are shown Pullman’s calculations. 


Table 2 
Position 1 2 4 5 6 
AA expt. (mp) 41 —21 —17 18 —16 
AX calc. C.A.C. [41] —15 —26 16 —37 
AA calc. Pullman [41] —6 —6 20 —13 


It is clear that not only are the directions of the shifts correct, but so also, to good 
accuracy, are the magnitudes. The shift at position 1 is the largest, in agreement 
with experiment. By no conceivable process could an agreement even remotely 
as good as this be obtained by any other interpretation of the long wave band. 

This is also in agreement with COOH substitution. In position 5 the shift _ 
is —15. This is less than +18 and of opposite sign. The opposite sign arises 
from the change of sign in 5«,, and the smaller value is consistent with a small 
conjugating (delocalization) effect in opposition to the inductive effect in COOH 
substitution, but not in —CH, substitution. Also, in position 6, the —COOH 
shift is +63. This is of the expected sign, and should be larger than 16, since 
now inductive and conjugative effects cooperate. But it is somewhat larger 
than we should have expected, showing that the situation is really a little more 
complicated than our simple analysis suggests. Despite this, however, there 
seems little doubt but that the displacements due to methylation have served 
to identify the transition with a degree of certainty greater than that for any 
other polynuclear molecule hitherto studied. 

A slightly unexplained feature is that in all the methyl substitutions of the 
next two bands (now in the ultra-violet and not the visible region) the shifts are 
always to long waves. We should have expected that some might have been to 
shorter waves. ‘They are also much less numerically, except that a phenyl group, 
where the inductive effect can only be small, conjugates freely and gives a large 
(40 mp) shift. Since a methyl must obviously conjugate much less than a phenyl 
this confirms that normally—i.e. except in odd-numbered rings—the conjugative 
effect is of the same order of magnitude as the inductive, but in odd-numbered 
rings it will often be of a smaller order of magnitude. 


§6. THE CONJUGATION (i.e. DELOCALIZATION) EFFECT 

All that remains is to discuss the conjugation effect, and to show that in all 
normal circumstances it leads to a red-shift. A discussion of resonance energy 
and bond length along these lines has already been given by Mulliken, Rieke 
and Brown (1941). But as we are interested only in spectroscopic properties 
we are able both to simplify the analysis and make it more generally applicable. 

In order to avoid all mention of electrostatic or inductive effects we must 
suppose that the substituent is of equal or similar electronegativity to the carbon 
atom where it is attached. The simplest assumption is that instead of —C=H, 
we have a vinyl group, —CH=CH,. A similar analysis can be given for other 
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non-polar substituents. Let us therefore discuss a system such as that in the 
diagram where B stands for any aromatic residue, not necessarily benzene. It 
will be convenient to number the carbon atoms as shown, and to suppose that the 
resonance integrals are 8 and fy: 


Let us suppose that the m.o.’s for B are known and have the form 
bj; =X ,€,56,5 ; ¥=3, 4, oe n+2, fom 3s 4, eee n+2 


with energies «=e; (measured relative to a carbon z-orbital energy as zero). 
What we want to study is the way in which the «; are changed when the vinyl 
group is substituted at atom 3. ‘To do this we express the m.o.’s for the total 
system as linear combinations of ¢, and ¢,, which are the atomic orbitals of the 
vinyl carbons, and the %,, which are m.o.’s for B alone. It will be recognized 
that this way of setting up the secular equations is a mixture of the usual method, 
where we employ only a.o.’s, and the LCMO method recently developed by 
Dewar (1949), where we employ only m.o.’s._ The secular determinant is 


ah Bo ° 0 
Bo —=¢ ° C338 Cea Soto” Go te 1G C3, n+2P 
C33 Cs eG 
C348 Cy — €. 
A= =0 
0 
Conte - “Ente € 


In this section, for convenience, we have neglected all overlap integrals. In the 
lower right-hand nxn section only the leading diagonal terms of A are non- 
vanishing. Expansion of this determinant gives 
fai" 
A =(c?— By) IT(e; —€) + ef* (ec; —€)2 eile : 
p Sr 
So the new energies are the roots of the equation 


2 
(ag Satis, 53a le | en (10) 
€,—€ 
The only circumstances in which this is invalid are those in which, for symmetry 
reasons, one or more of the levels of B are unaffected by substitution. We 
exclude this. It can be dealt with simply, if necessary, when it arises. 
Solutions of (10) are most simply regarded as given by the intersections of 
the two curves re 
y=ePX Vole —Pe)s 8 haat: (11) 


€3 


€ 
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There are 1 +2 roots, as we should expect from a system containing »+2 atoms 
in the conjugating framework. Figure 2 illustrates the intersection of these 
curves for the particular case of n=4. It can be verified from this—and is easily 
proved directly from (10)—that the following relationships hold between the 
original roots <«; and the new roots «,’: 


if €j >—-Po, then €; > &j 


les] <—Pos be, (ele; pak eee (12) 


€j << Bo ej < &. 
In addition to the roots (12) there are the ‘vinyl roots’ +>. These may be 
shifted either upward or downward, but the shifts will not be expected to be 
large. 
Now in the ground state of these molecules all the levels for which «<0 are 
doubly occupied, all the others being empty. Thus in the unsubstituted B of 


Fig. 2. Curves to illustrate the displacement of term values due to vinyl-substitution. The arrows 
along the axis of « denote the energy displacements. 


fig. 2 the lowest transition is of the type «,>«;. Quite clearly, since the result 
of vinyl substitution is to shift both «, and «, closer towards «=0, the energy 
difference «,’—e,’ is less than «;—e,. ‘This argument will hold for all these 
molecules provided only that |e,| and |e;| are less than fy. This is satisfied 
in all cases of interest. It follows that the delocalization effect is always in the 
direction of a shift to long waves. ‘This is the basic explanation of point (i) in 
the introduction. We can also see that this conclusion is unaffected if, for 
symmetry reasons, referred to above, one of the levels «4, ¢; is unchanged: for 
the other level will be changed, and in such a direction as to cause a bathochromic 
shift. ; 

There is a simple pictorial description of the phenomenon we have just been 
describing. We could think of the vinyl substitution as a conjugation of B and 
—CH=CH,. If we place the orbital energies <; of the original molecule B and 
the energies +, of the vinyl group on the same diagram the resulting energies 
can be thought of as derived from these unperturbed energies by supposing 
that levels from different ‘parts’ repel each other. It will be recognized that 
the rule (12) is merely a transcription of this. 
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Both experience and calculation show that the delocalization shifts are not 
large. We may therefore solve eqn. (10) by approximate methods. These show 
that, to first order, e,Btes? 

t j 34 


Gy se SS eel we | wi IS 
Soh fuk Ss ete) cae Pa Sat See Se) 
UF9 


A rough verification of this formula may be obtained, as follows: Since, for 
transitions in the same region of the spectrum, «,6? will not vary by a factor 
greater than about 2: 1, we should anticipate that a strong correlation should 
exist between e;’—«,; and c;,?._. Now if there are n carbon atoms in the unsubstituted 
molecule BH, the normalization condition shows that c,,? is of order 1/n. So 
we should expect that with alternant hydrocarbons the delocalization shift of 
frequency should vary approximately as 1/n. Individual c? values may differ 
considerably from the value 1/n, and so, where possible, we should consider 
the mean alkylation shift, when the methyl group is substituted, in turn, at all 
positions in the molecule. This we have done for a series of molecules, with the 
results shown in table 3. It is quite extraordinary that although the number of 
atoms in the molecule and the mean Av change by a factor of nine or more, the 
total spread of nAv values is less than 30°4 from the mean. This somewhat 


Table 3 
Molecule n Ute fein Mise Dy nAv Comments 
of spectrum (A) (cm) 

Ethylene 2 2000 3800 7600 (1) 
Benzene 6 2600 1100 6600 
Naphthalene 10 2800 530 5300 (2) 
Anthracene 14 3800 500 7000 (3) 

320 5800 (4) 
Benzanthracene 18 3600 { 50 4500 (5) 


(1) Comparison of ethylene and propylene ; (2) data in American Petroleum Institute 
Research Programme 44, serial numbers 88, 89, 130-2, ultra-violet absorption 
spectrograms ; (3) see Jones 1945; bands labelled F, G, H. ; (4) long wave band H ; 
(5) average of bands F, G, H, I in data of Jones 1940, and of Badger, Pearce and Pettit 
1952. 


unexpected result provides confirmation of the general description of these 
alkylation shifts that we have been investigating. We must, however, be careful 
not to conclude from table 3 that the conjugation effect necessarily far outweighs 
the inductive effect. ‘The fact, as represented by eqns. (9) and (13), seems to be 
that both effects have a magnitude which depends on c¢,”. It is this situation 
which is confirmed by the values calculated in the table. 


§7. CONCLUSION 
Our conclusions are best put in the form of answers to the five questions 
which were posed in §1. ‘They are: 

(i) Normal alkyl shifts in the long wavelength band are to the red as a result 
of greater delocalization of the z-electrons; this effect is more important than 
the electrostatic, or inductive, effect for alternant hydrocarbons. 

(ii) The variation in frequency shift with position of substitution depends 


essentially on the square of the coefficient (c,”) in the LCAO molecular orbital 


for the top occupied level. 
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(iii) In alternant hydrocarbons the delocalization effect appears to be larger 
than the electrostatic effect: and in fact this latter by itself would give rise to 
blue shifts. But in other molecules, including both non-alternant hydrocarbons — 
such as azulene, and heteronuclear molecules such as quinoline, the electrostatic 
effect is shown to be several times larger than before, and dominates the phenomenon. 

(iv) In non-alternant hydrocarbons and heteromolecules the electrostatic 
effect depends on the product of an induced change of coulomb term da, at the 
position of substitution, and the difference (¢,.”—¢,;") of the squares of the LCAO 
coefficients in the two orbitals involved. The variation with 6a, shows that the 
electron-donating groups (e.g. methyl) and electron-accepting groups (e.g. 
carboxyl) should give shifts in opposite directions. The variation with c,.?—¢,,? 
shows that even with one substituent the shifts will be in different directions 
with different positions of substitution. 

(v) A careful study of the predicted directions and relative magnitudes of 
the shifts due to methylation at all possible positions on the perimeter of the 
molecule can sometimes be used to characterize the electronic transition in an 
unambiguous way. A detailed analysis of the methyl azulenes gives a clear 
verification of the interpretation of the long wave band put forward by Mann, - 
Platt and Klevens on other grounds. 

Almost ‘all of the above is in accord with experiment. ‘There are a few 
exceptions, not properly understood, and there seem to be difficulties in applying 
these arguments to other than the longest wavelength allowed transitions. But — 
certain new features have been suggested, notably the variation of the mean 
shift for all possible positions in a given molecule with the increase of the number 
of carbon atoms. ‘This applies, naturally enough, only to alternant hydrocarbons. 

One final comment needs to be made.. We have where necessary assumed 
that first-order perturbation theory could be used. This is equivalent to the 
assumption that the shifts due to two or more substituents are additive. This 
appears quite satisfactory for alternants. Thus the successive shifts in the series 
benzene > toluene p-xylene are 78A and 75A respectively. In the methyl 
azulenes a slightly less strict additivity obtains; but it is certainly good enough 
for our purposes provided that we exclude—as obviously we must—cases of 
steric interference. 


Note added later 


Since this paper was communicated a paper has appeared by Longuet-Higgins- 
and Sowden (1952) in which part of the theory given here has been discussed from 
a somewhat different point of view. The same distinction is made between the 
inductive and conjugative effects, but in discussing the inductive effect the 
overlap integral is omitted so that the replacement of eqn. (6) by eqn. (9) for 
alternant hydrocarbons does not occur. Instead it is assumed that the inductive 
effect is actually zero. However, we have shown that the first-order effect with 
methyl substitution is necessarily hypsochromic. A further difference occurs 
with the hyperconjugative effect, where our argument is quite general, and that of 
Longuet- Higgins and Sowden is based on second-order perturbation theory. 


It is gratifying that in so far as the two accounts overlap they vee to sundlar 
conclusions. 


CC ee 
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ABSTRACT. ‘The half-lives, mass assignments and radiation characteristics of three 
neutron-excess gold isotopes have been determined as follows : 


200A 4 201A) 203A y 
Half-life 48 min 26 min 55 sec 
Max. f-energy (Mev) 2:2 1:5 1:9 
y-energy (Mev) 0-39; 1-13 0-55 0:69 


§1. INTRODUCTION 


HE only stable gold isotope is Au. 1*Au and !°Au are well known, 
| and have half-lives of 2:69 and 3-3 days respectively (National Bureau of 
Standards 1950). Sherr, Bainbridge and Anderson (1941) produced a 
48-minute isotope by a (n, p) reaction on mercury, and Maurer and Ramm (1942) 
found the same by a (n,«) reaction on thallium. Its mass number is therefore 
200 or 202. Butement (1951) found a new 27-minute isotope produced by a 
(y, p) reaction on mercury, and therefore of mass number 201 or 203. Its radiation 
characteristics were not determined. This paper describes work directed towards 
clarifying the properties of all the gold isotopes with mass numbers 200 to 204 
inclusive. 
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§2. PRODUCTION OF GOLD ISOTOPES 


The isotopes were produced by the appropriate irradiations of mercury metal — 
or thallium nitrate. The stable isotopes of these elements are as follows (National 
Bureau of Standards 1950): 


Mercury 196 198 199 200 201 202 204 
Abundance (%) 0:16 10-0 16°9 23°1 (NS ie 29-7 6°8 
Thallium 203 205 

Abundance (%) 29°5 70°5 


Fast neutrons were produced in the cyclotron by proton bombardment of a 
beryllium target. 20Mev protons gave neutrons, with energies up to 18 Mev, 
which were used to produce (n,p) reactions in mercury and (n,«) reactions in 
thallium. By increasing the proton energy up to 50 Mev, neutrons with energies 
up to 48 Mev were produced, and used to initiate also (n, pn) and (n, «n) reactions. 
X-rays of maximum energy 28 Mev, obtained from an electron synchrotron, were 
used to produce (y, p) reactions in mercury. 
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Fig. 1. Beta and gamma decay of ?°%Au. Fig. 2. Decay of gold isotopes from (y, p) 
reactions on enriched mercury isotopes. 


Except as described below, separation and purification of the gold isotopes 
from other activities were carried out by conventional radiochemical methods. 
The following half-lives were found for the gold isotopes by analysis of the 
decay curves. 

Mercury (n,p). 26 minutes, 48 minutes, and approximately 3 dae: the 
latter explicable as a mixture of 18Au and 1° Au. There was no period less than 
26 minutes except the tail end of a period of about 25 seconds which could just be 
observed in samples isolated very rapidly. This activity was not studied in 
detail, but its probable assignment is discussed at the end of this section. 

Mercury (n, pn) plus (n,p). ‘The above periods and also a 55-second activity 
(fig. 1). ‘This new short period first appeared when the maximum neutron 
energy was raised to 23 Mev, and the yield, compared with that of the 26-minute 
activity, increased as the neutron energy was raised successively to 28, 38 and 
48 Mev. 

Mercury (y,p). 55 seconds, 26 minutes, 48 minutes and 3 days. A full 
radiochemical identification of the 55-second activity as gold was not possible 
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in the time available. This activity was prepared by irradiating very pure 
mercury which was then poured on to a few milligrams of gold foil, distilled off, 
and the residue heated to redness, thus leaving only goldand platinum. _Consider- 
ation of the yield, which was comparable to that of the other radioactive gold 
isotopes showed that the 55-sec activity must arise from a (y, p) reaction, since in 
this region of atomic number (y,«) reactions leading to platinum, although just 
possible energetically, have a yield which is normally less than one hundredth 
that of (y, p) reactions. 

Thallium (n,v). 48 minutes. 

Thallium (n,n) plus (n,«). 26 minutes, 48 rainutes and 3 days, the latter 
showing the radiation characteristics of }®Au. 

Each of the above periods showed both beta and gamma activity. From the 
modes of production the following mass assignments can be deduced: 7° Au 
48 minutes, 2°%Au 26 minutes, ?°3Au 55 seconds. 

These assignments were confirmed by studying the gold isotopes produced 
by (y, p) reactions on enriched mercury isotopes. The concentrations (per cent) 
of the important mercury isotopes were as follows: 


Mass number 200 201 202 204 
Enriched 7°4Hg 11:8 Dye Ome 2-0) 
Enriched ?°°Hg 3-8 ibe Whibes ie7 
Enriched *°*Hg 4:8 AS) Wey eS} 


The yields for the three short period activities were computed after correction 
for decay during irradiation and processing, and compared with the results 
for natural mercury. 

Enriched °'Hg gave mainly a 48-minute gold (7°Au, fig. 2), the ratio of the 
yields of the 48-minute to the 26-minute activities being six times greater than the 
ratio from natural mercury, an increase in agreement with the change in the ratio 
of the concentrations of 7°'Hg and #°Hg. Enriched *°*Hg gave a nearly pure 
26-minute gold (?°1Au, fig. 2), the 48-minute activity being too small for accurate 
estimation, as might be expected, since the concentration of 7°'Hg was only 
1-8%. Enriched ?°*Hg gave mainly the 55-second activity (?°%Au), the 26-minute 
activity being too weak for very accurate measurement. 

A possible explanation of the absence of any period for ?°*Au, which must 
be produced in large yield both by Hg (n, p) and 'Tl (n,«), would be that it has 
a half-life almost equal to that of 48-minute ?°Au. ‘This explanation was shown 
to be untenable by irradiating equal weights of enriched °°Hg and natural mercury, 
for 20 minutes with 18 Mev neutrons under the same experimental conditions. 
The cyclotron beam current was constant during each irradiation, and a small 
correction was made for a change in the integrated neutron flux, as measured by 
the activities induced in standard copper monitor foils. ‘The gold activities were 
quantitatively isolated and the 48-minute activity from enriched mercury was 
only 0-18 of that from natural mercury. ‘This was to be expected if this activity 
derived only from a (n, p) reaction on the ?°°Hg present, the predicted ratio in 
this case being 3-8/23-1=0-16. If a 48-minute activity were also produced 
from 2°Hg with a similar reaction cross section, the total yield from enriched 
mercury would have been greater than that from natural mercury by a factor of 
about (91-5 + 3-8)/(29-7 + 23-1) =1°8. 

The decay energies of ?°?Au and Au calculated from the semi-empirical 
Fermi formula tabulated by Metropolis and Reitwiesner (1950) are 2-6 and 3-3 Mev 
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respectively. Using the graphs of Feenberg and Trigg (1950) relating logf with 
beta-decay energy, and a possible log ft value of 5-0, the calculated half-lives of 
these two isotopes could be 20 and 10 seconds respectively. The 25-second 
activitv is probably due to one of these isotopes. 


§3. RADIATION 
The gamma-radiation was studied with the conventional type of scintillation 
spectrometer (Nal-TI crystal), the pulses being displayed on an oscillograph and 
photographed. The sources were prepared by the neutron or X-ray irradiation of 
mercury, and photographs were taken at regular intervals in order to follow the 
fading out of individual components according to the half-lives of the isotopes 


responsible. The results for the gamma-ray energies were: ?°Au 1-13 and 0-39, 
201Au 0-55 and 73Au 0-69 mev. 
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Fig. 3. Microphotometer trace taken Fig. 4. Absorption of beta-radiation of 
across negative photograph of gold isotopes. 


oscillograph screen of gamma 
scintillation spectrometer. 


Figure 3 shows the record obtained by scanning with a microphotometer across 
a photograph taken just after the ?8Au had decayed away. There are peaks 
corresponding to the radiations from ?°°Au and ?*1Au, a broad peak due to the 
combined effect of the 0-156, 0-207 and 0-23 Mev gamma-rays of 1°°Au (Beach, 
Peacock and Wilkinson 1949) and a peak at about 75kev attributable to the 
combined effect of the 70 kev gamma-ray and the mercury K x-rays from !Au. 
The relative intensity of the low energy radiations is much exaggerated by the 
increased efficiency of the counter in this region. ‘The pulses in the region 0-6 to 
1-0 Mev are, in part at least, due to the Compton electrons from the 1-13 Mev 
gamma-ray, but there may also be superimposed the effect of one or more weak 
gamma-rays from 7°Au. . 

‘The beta-particle energy was determined by the usual technique of absorption 
in aluminium. For each isotope calculated corrections were made for the effect 
of the longerlived isotopes present. For ?°Au a separate sample was used to 
obtain each point on the graph (fig. 4). Because of necessary corrections only the 
initial part of each absorption curve could be accurately determined. From a 
knowledge of the efficiencies of the beta and gamma counters, it was estimated that 
only about 20%, of the beta-particle emission of 2°°Au is accompanied by gamma- 
rays, assuming the two gamma-rays are not in cascade. The corresponding value 
for ?'Au was 5%, and for 7%Au 10%. A reasonable approximation to the 
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maximum beta-particle energy could therefore be obtained by using the nomogram 
of Behrens (1950) relating the absorption in aluminium with the energy of a simple 
beta-emitter. 

The results for the beta-particle energies were 2Au 2-2, Au 1-5, 2°Au 
1-9mey. The value of Sherr et al. (1941) for 2Au was 2°5 Mev. 
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HET PERS Oe beh DECLOR 
Spin—Orbit Coupling in the D-D Reaction 


It is well known that the observed angular distribution of the outgoing particles in the 
‘D-D reaction cannot be accounted for adequately without the introduction of strong 
‘spin-orbit coupling (Beiduk et al. 1950, to be referred to as BPK). This has resulted in 
the suggestion that, because of such coupling, the outgoing particles may be polarized 
(Blin-Stoyle 1951, Cini 1951). Recently this possibility has been confirmed both for 
outgoing neutrons and protons (Bishop et al. 1952, Langley et al. 1952*); in particular at 
deuteron energies of 300 kev the proton polarization appears to be about 30%, 

A theoretical expression for the polarization at energies up to about 300 kev (when only 
incoming S, P and D waves need be considered) can be written 


A(E) sin 6 cos 0 

1+ A(E) cos? @ ’ 
where § defines the direction of motion of the outgoing particles in the centre of mass 
system and n is a unit vector normal to the plane of the reaction. A(E) is the energy 


dependent asymmetry factor occurring in the angular distribution as the coefficient of cos? 6. 
C can be written in the notation of BPK as 


P(#)=nC 


1 1 
C=—tiF [ (rr! Pl) +(710°)* (V1) + V6 (By2)* (39119 —2ya. t+ aD By) 


5 
= vi a oaths / ER ea enk (2) 


BPK attempted to calculate values for the various matrix elements occurring in C with 
- -the assumption that tensor forces alone were responsible for the spin-orbit coupling. 
‘Under these circumstances (2) reduces to 


C=254/6| Bill Yall sitadi(A:tS4o:), kaye anne (3) 
where ¢ is the phase difference between fj, and y,,’. 


* Prof. P. Huber has very kindly informed us that similar results have been obtained in Basel. 
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Assuming their calculated values for |fj,| and |y1,°| it is found that at 300 kev the 
maximum possible value of the polarization is about 7%, which is considerably less than the 
experimentally determined value. 

This fact can be interpreted in two ways : either (a) tensor forces alone are insufficient 
to account for the spin-orbit coupling, or (b) the theory of BPK is too crude to give © 
reasonable values for the matrix elements. 

In case (a) internucleonic forces having a spin dependence of the form L;;(o;+-0;) 
might be acting (such forces have been used fairly successfully in explaining nucle 
nucleon collision phenomena (Case and Pais 1950)). Using the same method as BPK for 
such an interaction yields the following expression for C 


5 
C=—-61.9 | 200") Ja ots | [(4r+4400. tee (4) 


It is then easy to show that values can be given to the various matrix elements which, 
whilst being consistent with the observed angular distribution, also easily account for the 
magnitude of the polarization. 

In case (0) the expression (3) might still be retained since it results mainly from symmetry 
arguments. An attempt can then be made to choose values for |/y2|, |y11°| and @ which 
will account for the form of the angular distribution and the sign and magnitude of the 
polarization. In this connection the behaviour of K, (see BPK), which reflects the effect 
of spin-orbit coupling on the angular distribution and total cross section, is of considerable 
interest. 

If only central and tensor forces are effective K, is given by 

K,=6| Yar |P+ | B12|2— 24/6 | B42 | Yu°| cos d; tee (5) 
and, according to BPK, should be energy independent. Its value can be determined from 
the measured partial and total cross sections. In practice, slightly better agreement with- 
the experimental results for outgoing protons is obtained if K, is assumed to vary with 
energy. This energy dependence can be explained if the fact that the outgoing F-waves 
(implied by £1.) have to penetrate a potential barrier is taken into account. We therefore 
write Bj.=P,1/"B,.’ where the penetrability P; is calculated by standard methods. 


Thus 
K,=6| yu°|?+ | Bie’ |?Ps—24/6| Bre" | | Yu | VP; cos d- tees (6) 
_ Ignoring the variation of P;!/2 with energy (which is very slow) and taking a mean value 
P,12, we can write 


KP aes Ff a Siegen (7) 


where X=6| y11°|?—24/6] Bre’| | yur°| Pst? cosh; Y= |fy2’|? and are taken to be energy 
independent. 

Using all ‘ thin target’ published results it is found that a good fit with experiment is 
obtained if X=0-024 and Y=0-087. 

Further, to account for the sign and magnitude of the proton polarization C=0-75. 

Such values of X, Y and C can be obtained by choosing |y1,°|=0-185, |B1.2| =0-295, 
and d=12° 35’. 

These values should not be taken too seriously since there are large errors on the values. 
of X and Y. The results do indicate, however, that tensor forces alone may well be 
sufficient to account for the spin-orbit coupling necessary to produce the observed 
polarization. If, subsequently, more quantitative theories show that this is not the case we 
have seen that a non-central force of the L . § variety might satisfy the requirements. 

The author would like to express his thanks to Mrs. M. Owen who investigated the 
fitting of K, with experiment and to the Pressed Steel Co. Ltd. for a research fellowship. 
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7th August 1952. 


Berpuk, F. M., Pruert, J. R., and KonopinskI, E. J., 1950, Phys. Rev., 77, 622. 

BisHop, G. R. , WESTHEAD, Ns M. , PRESTON, G., and Hasan, H., 1952, Nature, Lond., 170, it Fot 
BLIN- STOYLE, ‘Ri Je, 105 weeberoce "Phys. Soc. A, "64, 700. 

Case, K. M. , and Pais, A., 1950, Phys. Rev., 80, 203. 

Cintr, M. , 1951, Nuovo Cim. soap 

LANGLEY, El.) vererie eRe 'N,, and Stys, J. M., 1952, Phys. Rev., 86, 419. 


Letters to the Editor g5I 


The Nuclear Spins of °Ru and ‘Ru 


The paramagnetic resonance spectrum from a crystal of Ru(NH3),Cl, diluted with the 
diamagnetic salt Co(NH3)sCl; has been measured at a wavelength of 1:2 cm and a temper- 
ature of 20°K. The hyperfine structure indicates that the nuclear spins of 1°4Ru and ®*Ru 
are both 5/2, and that the ratio of their magnetic moments (101 :99)'1s1-09+ 0°03. 

Trivalent ruthenium has five electrons in the 4d shell and is analogous to trivalent iron 
in the 3d transition group. The spectrum and the magnetic behaviour of the Ru(NHs3)¢ 
complex will not be discussed in detail here, but briefly the behaviour seems to be consistent 
with that expected from a very strongly bound complex (Van Vleck 1935), the ground state 
being a Kramers doublet and the effective spin 4, which is a similar situation to that found 
in a ferricyanide. 

The hyperfine structure (see figure) has been resolved by diluting Ru(NH;).Cl; with 
the diamagnetic salt Co(NH;)sCl; in the ratio Ru: Co ~1: 200. The overall line width 
at 20° Kis then about 25 gauss. Since a width approaching 20 gauss is to be expected from 


® 
‘i 
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the field of the eighteen protons in the Ru(NH3), cluster, any appreciable improvement 
in the resolution could probably only be obtained by using deuterium rather than hydrogen 
inthe ammonia. In the figure the strong line is attributed to the even isotopes of ruthenium 
(Ru 96, 98, 100, 102, 104, total abundance 70-21%), and the satellites to **Ru (12°81%) 
~ and !°!Ru (16-98%). These abundances are given by Engelkemeir (1944). ‘The satellites 
consist of the outside pairs of lines AA’, G’G which can just be resolved, suggesting that 
these are from the two odd isotopes. The inner lines B and F are of the same total intensity 
(area under curve) as the outside pairs, within the experimental error, suggesting that these 
are also from both isotopes, but are unresolved. There is no sign of any hyperfine structure 
outside A and G. It is therefore concluded that the spins of both odd isotopes are 
I=5/2, [(2I-+1) lines from each]. Also the overall hyperfine separation AG of the most 
abundant isotope (1°'Ru) is seen to be slightly greater than that of **Ru (A’G’), giving the 
ratio of the nuclear magnetic moments [449;/499= 1:09. 
62-2 
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The usual shell model scheme suggests that the ground states of ®°Ru and 1°'Ru should 
be 4d;). and 5gy/2, respectively. However, as pointed out by Klinkenberg (1952), there is 
some evidence that the 5g,/. shell for odd neutrons fills in pairs, and he proposes that 
99°Ru with 55 neutrons may have 50 neutrons in closed shells with the remaining five in the 
configuration 5g7/2(4), 4ds5/.(1), while #°*Ru has its remaining seven in the configuration 
587/2(6), 4d5jo(1). ‘The present measurements are consistent with this proposal. 

Paramagnetic resonance has also been observed in K;[Mo(CN)s]; K;[MoCl,]; 
K,[MoCl;H,0]; [Ru(NHs)¢]Cls, 3HgCl,; [Ru(NHs3).]Cl, (concentrated); and in a 
mixed crystal of (NH,).[IrCl,] and (NH,).[PtCl,]. Investigations are being continued 
in these and other compounds of the Pt and Pd group elements. 

It was not easy to find a suitable ruthenium compound for this experiment, and it is a 
pleasure to thank Dr. Powell and Mr. Lever of Messrs. Johnson, Matthey & Co. Ltd. both 
for their valuable advice on the chemical side, and for supplying the ruthenium compound. 


Clarendon Laboratory, J. H. E. GrirFirus. 
Oxford. J. OweEN. 
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Nuclear Spin of Vanadium 50 


Vanadium contains an isotope of mass 50, whose natural abundance is about 0:2%. 
By observation of the paramagnetic resonance spectrum at 1:3 cm wavelength of a sample 
enhanced to about 20°%, the nuclear spin has been found to be six. This is the highest 
nuclear spin yet established. 

Single crystals were grown of K,V(CN),.3H.,O diluted with the isomorphous 
diamagnetic ferrous salt. ‘The crystals are monoclinic, with two ions in unit cell, one 
derived from the other by reflection in the crystallographic ac-plane. "The vanadous ion 
is in a 3d? state, and in the solid the orbital motion is almost completely quenched, leaving an 
effective spin S=3/2. The resonance spectrum of the abundant isotope, *!V (for which 
I=7/2), can be fitted to a Hamiltonian 


H =afH.S+D{S2—49(S+1)} + E(S2—Sy2}+ 48.1, 


where g is isotropic and equals 1-992, D= —0:0264 cm, E=0-0072 cm, and the hyperfine 
structure constant A is isotropic and equals —0:0056 cm“. ‘The relative signs are deter- 
mined from second-order effects in the spectrum, and the absolute sign of D from the 
intensity variation at low temperatures due to the difference in Boltzmann factors of the 
uppermost and lowermost states. 

With the sample enhanced in °°V, lines of about one-seventh the intensity of those due 
to >1V were observed for the M=} <> — } transition, whose relative separation corresponded 
to the known nuclear gyromagnetic ratio of °°V/*'V. From the frequency ratios 
(°V/?D) and (5'V/1H) given respectively by Walchli, Leyshon and Scheitlin (1952) and 
Knight and Cohen (1949) the ratio (°°V/°V) is 0-3792, and we find the ratio of the hyperfine 
structure constants A,;»/A,, to be 0:380+0-001. More accurate measurements of this ratio 
would be of interest to determine whether there is a nuclear volume effect, as, according 
to the theory of Abragam (1951), the hyperfine structure is due entirely to s-electrons, through 
admixture of the 3s3p%3d*4s state. 

The figure shows two weak lines, corresponding to the nuclear magnetic quantum 
numbers 5 and 6, of °°V, situated between the strong lines m=3/2, 5/2 of ®1V. Since the 
hyperfine structure lines are very nearly equally spaced, it is clear that there is no line 
corresponding to m=8 for ®°V, and the principal question was whether a line corresponding 
to m=7 would be resolved from the strong m=5/2 line of °'V. The position which such 
a line would occupy is indicated by the arrow in the figure, and it seems certainly to be 
absent. ‘lo check this point, experiments will be carried out at 3 cm wavelength, when 
a °°V, m=7 line would be further displaced from the °4V, m=5/2 line. This arises from 
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second-order effects which displace the lines downwards in field by (A2/2hy) {I([+1)—m?}. 
The displacement is different for the two isotopes, and at 3 cm wavelength an m=7 line 
would be as well resolved as the m=6 line shown in the figure at 1:3 cm. From other 
second-order effects it should be possible to check the positive sign for the nuclear moment 
of °°V obtained by Walchli, Leyshon and Scheitlin. 

For °°V, (Z=23, N=27) the odd proton and odd neutron should be both in f,), states. 
‘The observed spin J=6 illustrates a tendency for the spins of odd—odd nuclei not to take 


A fi 


i 


Spectrum of vanadium. Strong lines: m=5/2, 3/2 for *1V. The weak lines m=6, 5, 4 
(right to left) for ®°°V are indicated, and the position which an m=7 line would occupy is. 
shown by the arrow. 


the extreme values which the states of the odd particles would permit. Previous examples 
are ®Li, ??Na, ?4Na, °°Cl, 4°K, while (excluding cases where both odd particles are in $ states) 
only 1°B, ®°Rb have extreme values (cf. table II, Klinkenberg 1952). 

The authors wish to record their gratitude to the Director and the Electromagnetic 
Separator Group of the Atomic Energy Research Establishment, Harwell, who supplied 
the enhanced vanadium, and to Dr. P. F. D. Shaw, who converted it to the desired chemical 
form. 
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Oxford. B. BLEANEY. 
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Paramagnetic Resonance in Potassium Ferricyanide 


Potassium ferricyanide, K,Fe(CN),, has long been of considerable interest because of 
its anomalous magnetic properties. Unlike most ferric compounds, which show an 
isotropic susceptibility characteristic of the 3d°, °S state, the cyanide has a much smaller 
and very anisotropic susceptibility. Pauling (1939) has explained the small susceptibility 
on the hypothesis of covalent bonding, the set of octahedral bonds d’sp* to the six cyanide 
groups using up two of the 3d orbitals. Thus only three d-orbitals are available for the 
electrons which normally go in the 3d shell. ‘These are completely filled by six electrons, 
making potassium cobalticyanide and potassium ferrocyanide diamagnetic. With only 
five electrons, the single hole remaining in the 3d shell gives a susceptibility equal to that 
for S=4, with an orbital contribution. Van Vleck (1935) has shown that a similar result 
can be obtained by crystal field theory, and Howard (1935) has attempted to fit the 
anisotropic susceptibility using a rhombic component of the crystalline field comparable 
in magnitude with the spin-orbit coupling. 
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Potassium ferricyanide grows in monoclinic, pseudorhombic crystals, and resonance 
was detected at 20°K. In the concentrated salt a single resonance line is found, corres- 
ponding to a spin S=4, which reaches extreme positions along the rhombic axes. The 
g-values are: a-axis, 2:30; b-axis, 2:18; c-axis, 0-94. Here the accuracy (+ 0-03) is not 
very high because of the abnormally large line width. The half-width at half intensity 
due to magnetic dipole interaction with neighbouring ferric ions should be less than 0:2 ke, 
but the smallest observed width (along the a-axis) was 0:8 kc, while along the b- and- 
c-axes it was over 1:5kc. 'The line-width did not change on pumping down to 12°x, 
showing that the anomalous line width is not due to spin-lattice relaxation, since then the 
line width would be inversely proportional to the relaxation time 7, which varies rapidly 
with temperature. Such relaxation effects appear at higher temperatures, and 7 is estimated 
as about 210-19 sec at 62°K. This is very much shorter than that for normal ferric 
salts, as would be expected from the presence of low-lying excited levels. ‘The behaviour 
of the susceptibility shows that these are only a few hundred cm™ above the ground state. 

Preliminary measurements (Bowers, unpublished) with crystals diluted with K;Co(CN), 
show that there are two similar but differently oriented ions in unit cell, whose principal 
g-values (2:37, 2:11 and approx. 0°85) are not greatly different from those for the crystalline 
axes. The line width falls steadily with dilution, showing again that it must be due, not 
to spin—lattice relaxation, but to interaction with neighbouring paramagnetic ions. ‘The 
minimum so far reached ( = 30 gauss) is still much above that found by Bowers (1952) 
for diluted K,;Cr(CN),. 

The observed g-values are in good agreement with susceptibility measurements (Jackson 
1938, Guha 1951), extrapolated to 0°kK. Howard’s theory has been revised by Bersohn 
(1949), who finds that if the g-values for the a- and b-axes are fitted, that for the c-axis 
comes to 1:5. We find that the closest fit to all the principal g-values that can be obtained 
is the set (2:55, 2:27, 0:98). Fair agreement can then be obtained with the susceptibility 
measurements if the spin—orbit coupling coefficient is taken as about 390 cm™, which is 
only 10°% lower than the estimate of Kotani (1949). The fact that the theoretical g-values 
are slightly too high is due to the approximations in the theory. 

Preliminary measurements by Baker (unpublished) show that the g-values of the 
isoelectronic ion Mn?+ in K,Mn(CN),- 3H.O, diluted with the isomorphous diamagnetic 
ferrous salt, show a similar anisotropy. The line width also seems abnormally great along 
the direction of the lowest g-value. On the other hand, diluted crystals of K,V(CN),. 
3H,O show a normal spectrum with very narrow lines. ‘Thus the anomalous behaviour 
of the (3d)? ions (Mn?+, Fe?+) is in marked contrast with the normal behaviour of the 
(3d)? ions (V?+, Cr?+). The magnetic properties of the latter depend little on whether they 
are ionically or covalently bound, while those of the former are strongly affected. 

Our thanks are due to Mr. J. M. Baker and Mr. K. D. Bowers for allowing us to quote 
unpublished work. 
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The ‘Extra’ Levels in Rare Earth Salts 


Interpretations of the optical absorption spectra of rare earth crystals have sometimes 
shown that there are more low-lying levels than can be accounted for by assuming splittings 
by a crystal field (see Yost, Russell and Garner 1947 for a survey). It has been suggested 
that these may be due to a coupling between the electronic states and thermal vibrations. 
We would like to put forward evidence which points to another explanation. 
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We suggest that the levels may be created by the electrostatic fields in the crystals, but 
that the fields are not quite the same for all the ions. Thus the absorption spectra may be 
a superposition of the individual, slightly different, spectra. Magnetic resonance results 
on samarium ethyl sulphate would seem to require an explanation of this sort. It has been 
observed that a number of crystals, both concentrated and diluted with lanthanum ethyl 
sulphate, show two resonance lines with very similar g-values. The weaker of the two 
lines varies considerably in intensity from one crystal to another, and in fact it is possible 
to find crystals which show only one line. (A crystal of this last sort was used in obtaining 
the previously reported measurements by Bogle and Scovil (1952).) When the crystal was 
rotated in the magnetic field, the two lines, when present, did not cross as would be expected 
if the crystal was twinned. Also the absence of twinning was confirmed by the spectra of 
gadolinium and neodymium impurities, which were not doubled. There were thus two 
samarium ions with very similar ground states. A similar phenomenon has been reported 
in cobalt fluosilicate (Bleaney and Ingram 1951). 

Most of the optical absorption measurements seem to have been made on conglomerates. 
It would therefore be of considerable interest to make measurements on a number of single 
crystals to see if there are still too many levels. This work might then be followed by 
resonance measurements on the same crystals. 


Clarendon Laboratory, He Ee Ds Scoviu: 
Oxford. K. W. H. STEVENS. 


25th August 1952. 


BLEANEY, B., and INcrRaM, D. J. E., 1951, Proc. Roy. Soc. A, 208, 143. 

Boctez, G. S., and Scovit, H. E. D., 1952, Proc. Phys. Soc. A, 65, 368. 

Yost, D. M., Russeti, H., and Garner, C. S., 1947, The Rare Earth Elements and Their Compounds 
(New York: Wiley). 


The Mean Free Path of Conduction Electrons in Bismuth 


It has recently been pointed out by Sondheimer (1952) that in metals containing 
a very small number of conduction electrons, such as bismuth, the electronic mean free 
path 7 should be abnormally long, of the order of 1 u at room temperature, because the 
requirements of momentum and energy conservation forbid collisions between the electrons 
and the higher-energy phonons. ‘The experimental results of Eucken and Fo6rster (1934) 
on the apparent resistivity of thin bismuth wires in fact indicate a mean free path at room 
temperature of about 7 yu, while those of Justi, Kohler and Lautz (1951) on the thermo- 
electric power of thin bismuth films give a value of about 0:2 u. The discrepancy between 
these values may be due partly to the intrinsic difficulty, in thin-conductor measurements, 
of ensuring that the properties of the material are the same as those of the bulk metal; 
it should also be noted that the value deduced from the experiments of Justi et al. depends 
critically on the way in which / is assumed to vary with the energy of the electron. In view 
of the discrepancy between these values, further data seem desirable. 

It has been found that measurements on the ‘ anomalous skin effect’ in metals will 
yield reliable estimates of the temperature-independent quantity o//, where o is the d.c. 
conductivity, and hence of / at any temperature. ‘The quality factor Q of a resonant 
specimen is plotted against o/® (o being varied by varying the temperature), and the 
resulting curve compared with the theoretical curve (Reuter and Sondheimer 1948) 
of AX plotted against a/*, where AX ox Q and al§=(37w/109)"%(a/1) -V8oN?; here w is 
the angular frequency and o is measured in ohm~! cm~!. Thus o// can be found from 
the parameters fitting the experimental to the theoretical curve. 

Measurements have been made on bismuth at two frequencies, 1250 and 9060 Mc/s, 
and at three temperatures, 290°K, 80°K and 4°k, using the techniques described elsewhere 
(Pippard 1947, 1950, Chambers 1950, 1952). The 1250 Mc/s specimen was in the form 
of a single-crystal ‘ hairpin’ of wire about 6 cm long (cast in a silica tube), forming a 
shorted quarter-wave resonator ; a section from one limb of this specimen was used as 
a half-wave resonator at 9060 Mc/s. ‘The triad axis of the crystal was perpendicular to 
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the specimen axis. The residual resistivity at 4°k was about 0:07 times the room- 
temperature value. Because of the high resistivity of bismuth the resonators had rather 
low values of Q, which restricted the accuracy of the measurements. The figure shows. 
the observations fitted to the theoretical curve of Reuter.and Sondheimer ; it is assumed 
that on colliding with the surface of the metal the electrons are reflected diffusely rather 
than specularly (cf. Chambers 1950, 1952). The same value of o// has been used in fitting 
the results at both frequencies, namely 2:74 x 10? ohm~! cm~?, corresponding to a value 
of / at room temperature of 3-4. These figures differ considerably from those for normal 
metals ; for example, for copper the value of o// (found by the same technique) is. 
1-54. 101! ohm! cm~?, corresponding to a room-temperature value for / of 0-038 p. 


05;- © 1250 Me/s 


A 9060 Mc/s 


ak 


Experimental results fitted to theoretical (Reuter and Sondheimer) curve. 


The estimates for o/] and IJ are subject to the following uncertainties: (i) There is 
a possible error of + 20% due simply to the low accuracy of the measured values of Q. 
(ii) If the surface of the specimen is imperfect, the value of o// will be an underestimate 
and therefore that of 7 an overestimate. Under the microscope the specimen showed. 
a good surface finish ; experience with other metals (Chambers and Pippard 1952, 
Chambers 1952) leads to the belief that the true value of o// is highly unlikely to be greater 
than 4x10’ ohm cm~ ; correspondingly the true value of 7 at room temperature is. 
unlikely to be less than 2:3. (iii) These figures have been deduced by comparing the 
experimental results with the isotropic free-electron theory of Reuter and Sondheimer. 
For a single crystal, the values of o// and / obtained in this way depend on orientation and 
probably represent average values for electrons travelling at small angles to the surface 
(Pippard 1950, Chambers 1952): the values given above are therefore probably weighted. 
in favour of electrons travelling perpendicular to the triad axis. 

It is concluded that, in agreement with Sondheimer’s suggestion, the mean free path 
of the conduction electrons in bismuth is much longer than in normal metals, and that 
(for electrons travelling perpendicular to the triad axis) it probably lies between 2 and 4 uw 
at room temperature. 
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Zener’s Treatment of Ferromagnetism 


Zener (1951 I)* has suggested an interpretation of the effects of electron exchange 
interaction in solids which, as developed by Zener (II), (III) and (IV) and Carr (1952), 
appears to account for a wide range of structural and magnetic properties. 'The object 
of the present note is to consider critically the validity of this treatment. 

The basic postulates introduced by Zener (I), subsequently referred to as Z (i) to (iii), 
are : (i) ina single atom the electrons in the incomplete d shell are in the state of maximum 
multiplicity; (ii) the exchange interaction between electrons in adjacent d shells favours 
anti-parallel setting of their spins (i.e. antiferromagnetic); (iii) exchange interaction 
between d shells and the conduction electrons tends to a parallel setting (i.e. ferromagnetic). 

For a pure metal Zener expresses the interaction energy in the form 

B= toSq°—pSaSen soos es (1) 
where Sq, S, are the mean resultant spins of the d shells and conduction electrons in Bohr 
magnetons per atom and aq, f, y are appropriate coupling coefficients, all positive. The 
magnetic moment is then w=Sq+ So. 

Sa Ree .. ae Ol, WOE ae 
If 6?—ay<0, the minimum energy equilibrium condition is a5. 7 aSa =0, giving 
Cc 
S,=Sq=0, i.e. the resultant moment is zero. If however, f?—ay>0, the lowest energy 
occurs when either Sq or Sg, or both, attain their greatest possible value : 


for Bly < So max/Sa max : Sd= Sd max) So=(6/y)Sa max and w= Sq max(1 +B/y); ke) 
Bla < Sa max/Se max * Sa=(B/«)Se max» Se= Se max and w=S, max inp) &) 3 SGD) 
and «/B <.S¢ max/Sa max <6/y : Sa=Sa max» Se=Se max and 4=S¢maxt+Sa max: - - (4) 


E 
By writing the equilibrium condition as =(), Zener (I) obtains only case (2), but the 
c 


necessary condition for this case to apply may not be satisfied generally. ‘The validity of 
Zener’s (III) estimate of f/y is very doubtful, but this question will not be considered in 
this note. His estimate gives B/y ~. Since Sqmax may be as large as 5, case (2) will 
not apply if S_ max is less than about 2. 

For Fe, Co, Ni, assuming one conduction electron per atom, Sgmax=3, 2, 1. Then 
with £/y as calculated by Zener (III), (2) gives =3-45, 2:30, 1:17, as compared with the 
corresponding experimental values of 2:22, 1:71, 0-61. Zener (III), however, assumes 
the number of conduction electrons per atom to be 0:22, 0:7, 0:6. These numbers are 
consistent with the band scheme, but it is not clear how they arise in Zener’s scheme. 
If all the ions in the pure metal have similar d shell configurations, non-integral numbers of 
conduction electrons cannot arise. It is then necessary to assume that several types of 
d shell configuration occur (as Zener (IV) does for Fe), in which case Carr’s approach 
(see below) is more appropriate, and the calculations in Zener (III) have no significance. 

Carr (1952) considers the behaviour of a large range of ferromagnetic binary alloys. 
The sign of the interaction between d shells is taken to depend on whether the sum of 
their radii is greater or less than an arbitrarily chosen value. The values of the radii used are 
quite inaccurate (cf. Slater 1939, p. 349) and the errors may well be several times as large 
as the differences between adjacent tabulated values. Further, Carr assumes the sign of 
the interaction to be independent of the separation of the ions concerned, which seems 
very improbable for the range of alloys considered. Errors arising in these two ways will 
be particularly important just on either side of the critical value. ‘Thus the basis for this 
table is very weak, and there is no theoretical reason for assigning the signs to the interactions. 
in this particular way. It may be noted that the adoption of this table is contrary to Z (ii), 
and is therefore inconsistent with Zener’s general prediction of antiferromagnetism in 
non-conductors. 

Carr’s treatment of the equilibrium conditions is mathematically incorrect. A rigorous 
solution of this problem shows that the conditions given by Carr for the applicability of 
his solutions are not sufficient, and that additional solutions (i.e. values of Sg, Sg etc. for 


* To be referred to as Zener (I) and his other papers as Zener (II), Zener (III) and Zener (IV). 
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equilibrium) may be obtained. Since the values of the parameters involved are not known, 
however, a quantitative comparison is not possible. It may also be shown that an equally 
adequate interpretation is provided using the table of interaction coefficients but without 
recourse to Z (iii). 

Neither Zener nor Carr consider the implication of their treatment on temperature 
dependent properties. Discarding Z (iii), Carr’s model is similar to Néel’s (1948) model 
of a ferrimagnetic, for which the inverse susceptibility 1/y depends hyperbolically on T. 
As indicated above, a detailed development of the model for Fe, Co, Ni may lead to a 
similar conclusion. Experimentally, Fe, Co and the alloys FeNi, NiCo, FeV show an 
almost linear relation between 1/y and T, and the slight curvature appearing in the curves 
for Ni and some NiCr, NiCo alloys receives a natural interpretation on the band scheme 
(cf. Wohlfarth 1949). 2 . 

Thus without further modification Zener’s treatment, quite apart from the lack of a 
firm theoretical basis, does not lead to agreement with experiment, and Carr’s results cannot 
be considered as support for Zener’s postulates. 

I wish to express my thanks to Professor E. C. Stoner for guidance in this and other 
work, and to the Department of Scientific and Industrial Research for a grant. 


Physics Department, A. 'TEVIOTDALE. 
University of Leeds. 
25th August 1952. 
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The Thermal Neutron Capture Cross Section of *Sr 


When *Sr is irradiated in a nuclear pile, two radioactive isotopes are produced by 
neutron capture, ®*°Sr and a metastable state ®°™Sr. From the work of Sunyar et al. (1952) 
it follows that 86%, of **™Sr (T=70 min) decays by y-emission to the ground state of 
8°Sr (T=65 d), which in turn decays by K-capture and y-emission (0°514 Mev) to stable 
*>Rb. Thus if natural strontium is irradiated in a pile, measurements of the disintegration 
rate of the y-activity due to *°Sr and of the disintegration rate of the f-activity due to 
89Sr (T=52-5d) will enable the thermal neutron cross sections o of §*Sr and ®8Sr to be 
approximately compared. 

We have recently been using pile-irradiated strontium in biological experiments in 
which the longer life y-radiation from *°Sr has been most useful in verifying the Sr content 
of various specimens. In the course of these experiments the two disintegration rates 
mentioned above were measured. A rough comparison of the isotopic cross sections did 
not appear to support the value of 100 barns* for o®4Sr obtained by Pomerance and Arnette 
(1951). We decided therefore to compare these two cross sections more exactly. 

‘Specpure ’} strontium carbonate was irradiated in the Harwell Pile and, after allowing 
two days for the shorter-life y-activities of **™Sr and 8’™Sr (T=2-8 h) to decay, the 
f-activity of a sample of the activated salt (~ 2 mg) was compared with that of a standard 
source of **P using a GM4 counter with a fixed geometry. It was confirmed that the y-ray 
sensitivity of the counter was negligible, so that the y-rays from the decay of ®°Sr did not 
affect the counting ratio which was determined to within 1%. Since the variation of 
sensitivity of the counter for B-particles of 1:5 and 1-7 Mev is very small, the disintegration 
rate of the ®*Sr in the counting sample could be calculated. Absorption tests were used to 


verify that there was no appreciable contamination of the B-particle emission from each 
source. 


* Difficulty of estimating the purity of the samples of ®'Sr used produced large errors in this 
value. More recent estimation-would indicate o®Sr to be <3 barns (Pomerance (1952), private 
communication). 


+ High purity ‘Specpure’ salts were obtained from Messrs. Johnson, Matthey and Co. Ltd. 
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The same sample was then placed in a fixed position in a scintillation counter using a 
thallium-activated sodium iodide crystal. The y-ray activity of the sample was compared 
successively with that from a standard source of “!I with a mixed y-ray spectrum but 
predominantly of energy 0:36 Mev and !’Cs (y=0-663 Mev) supplied by the Isotopes 
Division at Harwell. The efficiency of the counter was thus interpolated for a y-ray energy 
of 0-51 Mev, and so the disintegration rate of *°Sr was obtained. 

Accepting the value of 0-0050 barn for the isotopic capture cross section of *8Sr (Seren 
et al. 1947), the value obtained for the isotopic capture cross section for ®Sr was 
0-32 barn+ 20%. As a test of the purity of the irradiated strontium salt, the B- and y-ray 
activities were determined at various times over a month. The half-life of the B-activity 
was 52°5 days (**Sr) and that of the y-activity 67 days. 

If the present value of o is assigned to **Sr of relative abundance 0:55°% the value 
of o for the stable isotope §’Sr can be calculated. Using the value of Colmer and Littler 
(1950) for the total capture cross section of natural Sr for thermal neutron (c=1-35 barn) 
together with the value of 1-3 barn obtained for 87™Sr by Seren, Friedlander and Turkel 
(1947), the value of o for §7Sr is about 17 barns. 


M.R.C. Radiobiological Research Unit, G. E. Harrison. 
Atomic Energy Research Establishment, F. D. SEYMoUR. 
Harwell, Berks. 
17th September 1952. 
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The Origins of Background Tracks in Shielded 
Electron-Sensitive Emulsions 


In a previous communication (Fremlin and Walters 1950) it was shown that a 200 u 
Ilford G5 emulsion accumulated electron tracks at the rate of 31+3 mm~? per day when 
screened within 2-inch lead walls at a depth of 600 metres below ground. By the use 
of cleaner shielding it has been possible to reduce this background accumulation rate to 
2:6+0-1 mm~? per day. This is considerably larger than the rate to be expected from 
cosmic radiation, which is estimated from the measurements of Miyazaki (1949) to be of 
the order of 10-* mm~? per day, and experiments have now been made to discover what 
proportion of the background arises from radioactive contamination of the plates themselves. 


Observed rate of formation 


Origin of tracks B (mm? per day) 


Thorium present in glass backing 0:57+ 0-05 
40K present in glass backing 1:00 + 0-05 
14C present in gelatine 0-83 + 0-06 
Shielding (entering plate at upper surface) 0-14+ 0-02 
Total 2:54+ 0-10 

Total track density (separate count) 2:61 + 0-09 


As previously pointed out, the thorium and potassium contamination of the plate 
backing and the C present in the emulsion layer may make an appreciable contribution 
to the background when external radiation has been reduced to a minimum. ‘This has 
now been confirmed by an analysis of the spatial distribution in the emulsion layer of 
tracks formed during a two-month exposure. ‘The analysis showed that a considerable 
proportion of the observed background arises from the glass backing, since the number of 
electron tracks intersecting the lower emulsion surface was found to be a factor of 2:4 greater 
than the number crossing the upper surface. The results of the analysis are shown in the 
table. ; 
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The thorium activity of the glass was estimated from the results of a count of «-particle: 
tracks observed at the emulsion-glass interface, the difference between this and the total 
glass activity being ascribed to *“K. The potassium content of the backing glass was. 
found by analysis* to be 0:15°%, while the activity quoted in the table is estimated to. 
correspond to a concentration of 0:09°%4. C present in the gelatine will produce tracks 
at the rate of 0-76+ 0-07 mm~? per day if present in the concentration of 0-95 x 10-™ g/g 
of carbon (Anderson and Libby 1947). In the table all low-energy tracks starting within 
the emulsion layer have been ascribed to the decay of !C. 

The analysis indicates that the minimum background so far attained cannot be reduced 
by more than 10% with improved shielding if the normal G5 plate is used. A further’ 
reduction can, however, be achieved by coating the emulsion on to a support free from 
potassium and thorium. This has been verified using +,-inch Perspex sheet as a backing. 
During a four-month exposure these plates accumulated tracks at the rate of only 
0:86+ 0-11 mm~ per day in a 200 w emulsion. This background may be attributed almost 
entirely to the 4C present in the gelatine. 

It is considered that the use of such plates for the detection of rare events such as. 
double B-decay (Fremlin and Walters 1952) will enable the limiting half-life to be increased 
by as much as a factor of 10, since the comparatively short range of the tracks resulting 


from the decay of !C will, in many cases, permit them to be distinguished from the events. 
under investigation. 


The Physics Department, M. C. WALTERS. 
University of Birmingham. 
11th August 1952. 


ANDERSON, E. C., and Lissy, W. F., 1947, Phys. Rev., 72, 931. 


FREMLIN, J. H., and Watters, M. C., 1950, Proc. Phys. Soc. A, 63, 1178; 1952, Ibid., 65, 911. 
Mrvyazakl, Y., 1949, Phys. Rev., 76, 1733. 


* We are much indebted to Messrs. Chance for making this analysis. 


REVIEWS OF BOOKS 


Les machines a penser, by Louis CourFiGNaL. Pp. 153. (Paris: Les Editions. 
de Minuit, 1952.) No price given. 


The author defines a thinking machine as a machine which can replace a human being” 
in the performance of operations belonging to some class, however limited. Calculating” 
machines clearly come into this category and their development is described beginning with 
the work of Pascal and ending with modern large-scale digital and analogue machines. 
Some comparisons are drawn between calculating machines and the cerebral cortex, both 
in relation to action and structure. This leads to a discussion of symbolic logic and ways of 
mechanizing it, and finally to a chapter on the problem of thought—but the net result is, 
as the author admits in a postscript, the statement of a series of questions to be answered by 
future research rather than a conclusion. 

The author’s approach to the problem of thought is to consider the brain as a machine 
which can derive systematically the consequences of propositions expressed in symbolic. 
form and held in its memory, and which can store these consequences and use them as new 
data for a further series of deductions. He considers that the logical rules followed may not- 
be the same as those of Aristotelian logic, and that light may be shed on them by a statistical 
examination of texts in a wide range of subjects in order to see what concepts are normally 
linked together. A diagram is given showing how certain logical processes can be carried 
out in a machine composed of elementary logical and storage elements. A comparison—- 
which may perhaps strike some readers as a little naive—is made between the connections in 


this machine and the arrangement of nerve fibres in the brain. The philosophical aspects. 
of the subject are not embarked on. 
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The author appears to have been less influenced by a consideration of the properties of 
large universal calculating machines than have other writers who have stressed the ability 


-of such machines to discriminate between different kinds of stimuli received through their 


input mechanisms and to modify their own programmes of instructions accordingly. In 
fact the amount of space devoted to universal calculating machines is sufficient for little more 
than a general description. 

It would have been helpful if more references had been given to the sources on which the 


arguments are based, and if an index had been provided. Apart from this, however, the 
book is well written. M. V. WILKES. 


An Introduction to Mathematical Physics, by R. A. Houstoun. Pp. x+262. 
(London: Blackie, 1952.) 25s. 


The familiar pre-war textbook of the same title which was first published in 1912 has 
been completely revised and extended in scope. In addition to chapters on attractions, 
hydrodynamics, heat conduction, wave motion, electromagnetic theory and thermodynamics, 
there are two new ones on quantum theory and relativity. As is to be expected in a book 
of this size, only the elementary parts of these subjects are included. 

Where necessary, each topic is introduced by a short statement of the basic physical 
principles on which it rests, and, though these statements are rather over-simplified in places, 
they are in general sufficient introduction to the mathematical development which follows. 

The reader is expected to have a fair knowledge of calculus, but the analysis required 
is derived at appropriate places in the text. Because of its extremely clear and straight- 
forward presentation, the book will be of great assistance to physics students (for whom it 
was mainly written). ‘ 

There are over 100 examples. Vector notation is not used. A. N. GORDON. 


Fourtersynthese von Kristallen und thre Anwendung in der Chemie, by WERNER 
Nowacki. (Lehrbiicher und Monographien aus dem Gebiete der exakten 
Wissenschaften. Chemische Reihe, Band VI.) Pp. 237. (Basel: Verlag 
Birkhauser, 1952.) 34.30 fr. 


Professor Nowacki’s book inevitably invites comparison with Dr. Booth’s Fourier 
Technique in X-ray Organic Structure Analysis. Cursory examination of the tables of 
contents shows that the topics considered are almost identical, but there is a considerable 
difference in the approach. The difference does not quite reach that implied in the old 
story of the competition for the best book on the elephant, in which the English entry was 
entitled Elephants I have Shot and the German Der Elefant als Grundbegriff in Philosophie 
und Metaphysik, but the tendency is in that direction. Without intention of disparaging 
either book, it might be said that Dr. Booth is concerned with what he has himself found 
useful and Professor Nowacki with what might be useful to someone. 

The first section (99 pages) is concerned with electron-density synthesis. Sixteen 
‘pages are devoted to explaining the use of Professor Lonsdale’s structure-factor and 


-electron-density formulae, with special reference to the space groups Pna (four pages) 


P2,2,2, (four pages) and P2,/c (two and one-half pages). Projections of the whole cell or 
parts of it are treated. There is considerable discussion of the effect of errors, both experi- 
‘mental and termination-of-series, and the various suggestions for avoiding them are 
described. The general tendency is to discourage the use of an artificial temperature 


factor as introducing greater errors and uncertainties than the ‘ diffraction’ effects it is 


intended to remove. This conclusion might be anticipated from the known mathematical 
theorem that, of all finite trigonometrical series of given degree, the terminated Fourier 
series is the best mean-square approximation to the function to be represented. Other 
methods, such as continuing the series with calculated F’s, are not subject to the same 
criticism. Methods of refinement (differential syntheses, least squares, steepest descents) 
are fully discussed, as well as methods of structure determination (trial and error, iso- 
morphous replacement, heavy atom, algebraic methods, inequalities) and related questions 
such as the determination of absolute scale. 
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Patterson syntheses, including Harker and Buerger syntheses, occupy 78 pages. 
Dr. Booth devotes less than six pages to these topics, as “‘ only on rare occasions have the 
Patterson or Patterson—Harker syntheses been of service”’. The difference is partly due to 
Professor Nowacki’s full discussion of space-group determination, based on Buerger’s 
diffraction symbols. (Table 19, a more detailed version of these, is supplied as a separate 
booklet in a pocket of the cover. One wonders whether the convenience of this is out- 
weighed by the risk of complete loss.) If enantimorphous pairs are counted as single space 
groups, x-ray (Laue) symmetry and absences alone are capable of resolving only 68 groups out 
of the total 219. Crystal class in addition would resolve all but eighteen groups, occurring 
in pairs. A detailed study of the heights and positions of the peaks in Patterson (or Harker 
or Buerger) syntheses would in principle resolve even these, but overlapping of maxima, 
which so often reduces the usefulness of Patterson syntheses, might make the resolution 
difficult in practice. The statistical methods of space-group resolution were developed too 
late for mention in Professor Nowacki’s book; when applicable they would involve much 
less labour than the Patterson-synthesis methods. 

The third and last section (46 pages) describes methods of numerical calculation, 
including several types of strips and analogue and digital computers. ‘The use of strips 
appears to be given in sufficient detail to make reference to the papers of their originators 
unnecessary, but most of the computers are described only very briefly. Optical methods 
(‘ x-ray microscope ’ and ‘ fly’s eye’) are described in fair detail. 

A few misprints were noted : a nicht seems to have been omitted in the definition of 
gemischte Indizes on p. 25, the bar over h is omitted (or has failed to print) in some places on 
p. 29, and z is printed for y on p. 144. The reviewer found a number of abbreviations 
(u.a., u.a.m., u.U., u.W., z.T.) puzzling, but they would presumably be clear to readers with 
a better knowledge of German. Still another meaning, this time |F|*, is given to the 
expression ‘ structure factor’. 

The book is highly recommended as a full, clear and detailed account of the matter within 
its field. A. J. C. WILSON. 


Nuclear Data, Supplement 3 (January 1951 to July 1951) to National Bureau of 
Standards Circular 499. Pp. 66. (Washington, D.C.: U.S. Department 
of Commerce, 1952.) Price $4.25 includes Circular 499 and three 
supplements. 


A Guide to Tables of the Normal Probability Integral, National Bureau of 
Standards Applied Mathematics Series 21. Pp.iv+16. (Washington, 
D.C.: U.S; Department of Commerce,.1952,). Ive: 


Mass Spectrometry, Report of a Conference organized by The Institute of 
Petroleum. Pp. vi+205. (London: The Institute of Petroleum, 1952.) 
30s. 


CORRIGENDUM 


A Theoretical Attempt to Predict the Excited States of Nuclei in the Neighbourhood ~ 
of Pb, by M. H. L. Pryce (Proc. Phys. Soc. A, 1952, 65, 773). 


The statement (§ 6, p. 784) that the interaction energy of a hole and a particle has the 
_same sign as for two particles is incorrect. 'The proper statement is that the interaction 
energy is given by the same formula as for two particles, provided «, and ¢, are replaced by 
— («t+ €s) and 3(€,—3e,) respectively. The conclusions regarding the doublet configuration 
in *°*T'] remain unaffected, but the order of levels in the quartet should be reversed. The 
agreement with the intensities of the alpha groups of 7!*Bi is thereby spoilt. 


. 
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ABSTRACTS FOR SECTION B 
X-Ray Analysis of Proteins, by W. L. Brace. 


ABSTRACT. ‘The attempt to discover the atomic arrangement in the protein molecule 
seems very ambitious. Ever since Bernal first showed that crystals of protein gave x-ray 
diffraction pictures, it has been clear that a protein molecule of a given type is a structure 
with a definite individual form; the x-ray diffraction spots are very sharp and reproducible 
and extend to regions corresponding to interatomic distances. "she molecules are, however, 
of great complexity. It has been a triumph of x-ray analysis to pass from simple substances 
like rock salt to such molecules as strychnine or penicillin with about one hundred atoms. 
We are now trying to analyse a molecule such as haemoglobin, which contains ten 
thousand atoms. ‘The reward, if an analysis were completed, would be great, because the 
determination of any one protein would undoubtedly cast a flood of light on the character 
of these bodies, which Nature has selected as the basis of living matter. 

Work on protein is going on in a number of centres. I wish to speak here in particular 
about the work on haemoglobin which has been going on in Cambridge for the last twelve 
years under the direction of Dr. Perutz. More x-ray data for haemoglobin are available 
than for any other protein. ‘Though still a long way short of the goal, considerable 
advances have been made recently, so great as to encourage us to hope that the problem 
is soluble and even that a solution lies not far ahead. 


The Dynamics of Small Transient Cavities, by W. D. CHESTERMAN. 


ABSTRACT. Experimental methods have been developed for creating and observing 
single spherical cavities in water and certain organic liquids. Rapid photography permits 
motion analysis to be made. 

The cavities show several cycles of expansion and collapse during their lifetime of a few 
milliseconds. At each volume minimum a sharp pressure maximum occurs in the liquid. 

It appears that these cavities form initially on minute gas nuclei in the liquids, and the 
existence of these nuclei is essential for the type of cavitation studied. ‘The relation of 
cavity life to maximum diameter is discussed. 
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The Temperature Variation of the Magnetization of Nickel in Low and Moderate | 
Fields, by R. S. TEBBLE, J. E. Woop and J. J. FLORENTIN. | 


ABSTRACT. An account is given of measurements on the reversible changes in magneti- 
zation accompanying change in temperature in annealed nickel, and the results are discussed 
in relation to the work of Bates and others on the temperature changes produced by adiabatic 
changes in magnetization (the magneto-caloric effect). It is found that at relatively high 
fields, where the magneto-caloric effect is mainly reversible, the results are in good agreement 
with those obtained from direct magneto-caloric measurements, while for low field values a _ 
quantitative separation of the reversible and irreversible components of the magneto- 
caloric effect has been carried out. The result of an analysis of the reversible heating effect 
is compared with the estimate given in the theoretical investigation of Stoner and Rhodes, 
and it is suggested that, over the low field range considered, rotational processes predominate, 


The Effect of Temperature of Deformation on the Electrical Resistivity of Cold- 
worked Metals and Alloys, by 'T. BRoom. 


ABSTRACT. A special apparatus has been used to draw wires of various metals at 
temperatures between —183°c and 100°c, and to measure their resistances at the tempera- 
ture of drawing. The change of resistivity with deformation was found to depend on the 
temperature of deformation relative to the recrystallization temperature. The results 
explain qualitatively the apparently anomalous resistivity increases previously reported 
for tungsten and molybdenum. It is suggested that a unified theory of the effect of deform- 
ation on the resistivity of both pure metals and alloys can be based on stacking faults. f 


Yield Points in Bending Experiments on Zinc Crystals, by M. J. DUMBLETON 
and B. W. How ert. 


ABSTRACT. It is shown that sharp yielding can be observed in previously unstrained 
zinc crystals by means of bending experiments devised to avoid premature overstraining 
of the specimen. ‘This resolves a difficulty in the interpretation of previous tensile experi- 
ments on crystals which show the yield phenomenon. 


Some Observations on the Yield Point in Zinc, by H. L. Warn. 


ABSTRACT. _ Evidence is given confirming that the yield point in zinc is due to nitrogen 
in solution. ‘The theoretical prediction that a yield point should appear in a susceptible 
crystal without the necessity for any prior deformation is verified experimentally, and an 
explanation is given of the observation that a yield point may appear in a specimen even 
though free dislocations have been introduced subsequent to ageing. 

Under certain conditions the yield point values in zinc decrease as identical strain-ageing 
treatments are repeated, and a smooth stress-strain curve is ultimately obtained. This is 
attributed to an accumulation of work hardening which causes the dislocation density to 
increase to such an extent that the dislocations can no longer be effectively anchored by the 
available nitrogen atoms. ‘The difference between this behaviour and that of iron on 
repeated strain ageing is discussed. 


A Critical Study of the Asymmetrical Temperature Gradient Thermoelectric Effect 
in Copper and Platinum, by N. FuscHiLuo, 


ABSTRACT. A method is described by which metallic wires can be subjected to easily 
measurable but very pronounced asymmetrical temperature gradients without straining the 
wire. ‘The electromotive forces developed by the application of symmetrical and asym- 
metrical gradients to annealed ‘ thermopure ’ platinum wire and specially prepared cables of 
copper wire, at both high and low temperatures, have been measured at many points along 
considerable lengths of wire. Precautions to minimize or eliminate spurious effects 
enabled measurements significant to 0-01 pzv to be taken. 

The results provide no justification for the supposition that an asymmetrical temperature 
gradient thermoelectric effect exists, under the conditions of investigation, which cannot be 
explained in terms of inhomogeneity effects. The nature of these effects is discussed and 
some evidence is provided which indicates that the variation in the Thomson thermoelectric 
potential with strain and impurity is significant. 


